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ABSTRACT 


This  paper  documents  the  development  of  a  non-para— 
metric  statistics  package  for  the  TRS-80  microcomputer. 

The  package  is  comprised  of  ten  programs  with  the  major 
emphasis  on  non-parametri c  hypothesis  testing. 

The  programming  language  is  TRS-80  Level  II  Disk 
Basic.  The  package  is  compatible  with  any  disk  based 
TRS— 80  Model  I/III  compatible  microcomputer .  With 
modification  of  the  screen  display  and  the  disk  access 
commands,  the  package  is  transportable  to  microcomputers 
produced  by  other  manuf acturers. 

The  statistical  analysis  capability  implemented  on  a 
relatively  inexpensive  system  provides  a  useful  tool  to  the 
student  or  the  trained  analyst. 
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I.  INTRODUCTION 


Statistical  Analysis  of  data  is  utilized  as  a  decision 
aid  in  most  every  field.  Personnel  in  technical  fields 
have  been  exposed  to  and  require  the  use  of  many 
statistical  analysis  procedures  during  the  course  of  their 
career.  In  a  career  in  which  statistical  analysis  of  data 
is  required  ,  but  not  on  a  daily  basis,  the  knowledge  of 
the  proper  statistical  procedure  is  often  retained  where 
the  mechanics  of  implementing  the  procedure  requires 
extensive  review.  A  software  package  which  eliminates  this 
review  and  handles  the  mechanics  of  implementing  the 
procedure  is  most  valuable. 

The  personal  computer  has  evolved  to  the  degree  that 
it  is  now  well  suited  for  the  implementation  of  statistical 
analysis  procedures  for  small  to  moderately  sized  data 
bases.  The  statistical  analysis  software  that  is 
commercially  available  is  generally  unsophisticated  and 
unsuited  for  the  serious  user. 

An  excellent  package  which  models  parametric 
distributions  and  implements  parametric  hypothesis  testing 
was  developed  by  Mr.  R.  P.  Isbell  and  was  the  subject  of 
his  thesis  work  at  the  Naval  Postgraduate  School.  The 
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purpose  of  this  thesis  is  to  provide  a  software  package 
which  implements  the  more  commonly  used  non— parametric 
procedures,  accurately  approximates  their  distributions, 
and  provides  a  convenient  data  base  manager  useful  for  both 
this  package  and  that  developed  by  Nr.  Isbell. 


B 


II. 


OVERVIEW 


A.  PACKAGE  CONTENTS 


The  package  may  be  logically  divided  into  three 
partitions;  the  data  base  manager,  non-par ametric 
distributions,  and  non— parametric  hypothesis  tests.  The 
data  base  management  programs  permit  data  entry  and  storage 
to  be  in  a  columnar  -format  or  in  a  tabular  -format.  The 
programs  provide  a  convenient  -file  editing  capability  as 
well  as  a  flexible  file  printing  capability. 


The  more  commonly  used  non-par ametric  distributions 
and  their  inverse  distributions  are  approximated.  The 
distributions  are: 

1.  Wilcox son  signed  rank  distribution 

2.  Mann-Whitney  distribution 

3.  Smirnov  distribution 

4.  Kolmogorov  distribution 

5.  Lilliefors  modification  to  the  K-S  test  for 
Normal  Distribution 

6.  Lilliefors  modification  to  the  K-S  test  for 
Exponential  Distribution 

7.  Chi— Square  distribution 


The  package  contains  procedures  for  implementing  the 
following  hypothesis  tests: 


1. 

r> 


o  . 

4. 


5. 


Wilcoxson  signed  rank  test 
Mann-Whitney  test 
Smirnov  test 

Test  for  Normal  Distribution 

a.  Mean  and  variance  unknown 

b.  Mean  known,  variance  unknown 

c.  Mean  and  variance  known 
Test  for  Exponent i al  Distribution 
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a.  Parameter  unknown 

b.  Parameter  known 

£>.  Contingency  Table  Chi-Square  test 

B.  MICROCOMPUTER  SYSTEM 

The  package  was  developed  on  an  LNW  Model  I 
microcomputer  system  with  4SK  bytes  of  memory  and  one  5-1/4 
inch  floppy  disk  drive.  The  package  is  totally  compatible 
with  any  TRS— SO  Model  I  or  Model  III  disk  based 
microcomputer  system  and  is  not  dependent  on  any  unique 
features  of  the  advanced  TRS-BO  Model  I/I I I  operating 
systems. 

The  programming  language  is  TRS-80  Disk  Basic 
implemented  by  Microsoft.  Compatibility  with  all  TRS-80 
Model  I/I  1 1  systems  is  maintained  by  using  only  TRSDOS 
compatible  commands.  The  package  is  not  directly 
transportable  to  non— TRS— 80  compatible  computers;  however, 
the  Microsoft  interpreter  is  commonly  used  by  other 
manuf acturers.  The  display  commands  would  have  to  be 
altered  to  be  compatible  with  mi crocomputers  produced  by 
other  manuf acturers. 

C.  LIMITATIONS 

With  any  program  development  on  a  microcomputer ,  the 
programmer  must  be  cognisant  of  its  speed,  memory,  and 
precisian  limitations.  To  achieve  maximum  utilization  of 
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the  computer  memory,  the  -floppy  disk  is  used  as  a  virtual 
memory  device  pulling  in  routines  on  an  as  needed  basis. 
This  method  permits  a  minimum  o-f  3000  data  observations  to 
be  resident  in  the  computer  memory. 

The  speed  o-f  execution  is  a  major  concern  when  using 
an  interpreted  language.  To  achieve  an  acceptable 
execution  speed,  asymptotic  distributions  were  used  in  lieu 
o-f  exact  distributions  and  faster  but  more  memory  intensive 
sort  routines  were  routines  were  used. 

Precisian  is  certainly  a  concern  in  this  package; 
however,  single  precision  variables  proved  to  be  adequate 
in  most  cases.  Double  precision  variables  were  used  only 
in  the  approximations  to  the  Wilcox son  and  the  Mann— Whi tnev 
distributions. 
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III.  DISCUSSION  Or"  PROGRAMS  AND  ALGORITHMS 

A.  STRUCTURE 

The  package  consists  of  ten  programs,  i.e.  a  control 
module,  two  data  I/O  programs,  a  program  calculating 
distributions  and  six  programs  performing  hypothesis  tests. 
The  choice  of  the  program  to  be  executed  is  made  through 
the  menu  driven  control  module .  The  user  is  free  to 
execute  any  program  in  the  package  through  the  control 
module  with  all  variables  being  initialized  as  a  new 
program  is  executed. 

All  data  I/O  is  conducted  through  the  two  data  I/O 
programs.  The  data  is  entered  to  the  hypothesis  test 
programs  through  a  menu  of  established  data  files  which 
were  created  by  the  data  I/O  programs. 


B.  PROGRAM  1,  MENU/BAS 

MENU/BAS  is  the  control  module  from  which  all  programs 
are  executed.  The  program  consists  of  the  master  menu  and 
two  sub— menus.  The  functional  flow  of  the  program  is  shown 
in  figure  1. 
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BLOCK  DIAGRAM  Oh  tlENU/BAS 


FIGURE  1. 
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C.  PROGRAM  2,  DBMGR/BAS 


DBMGR/BAS  provides  the  data  I/O  -for  the  hypothesis 
test  programs.  The  program  is  menu  driven  and  offers  the 
following  selections: 

1.  Create  New  Data  Files 

2.  Edit  Data  Files 

3.  Purge  Data  Files 

4.  Print  Data  Files 

5.  Return  to  Master  Menu 

A  special  purpose  data  file  named  EDITDATA/DIR  is  used 
by  this  program  to  contain  the  filenames  of  the  data  files 
created.  It  is  used  as  a  data  file  directory  and  is  called 
by  the  hypothesis  test  programs  to  display  the  filenames  of 
all  current  data  files  and  to  permit  a  convenient  menu 
selection  of  the  data  files. 

C.l.  NEW  FILE  CREATION  ROUTINE 

The  new  file  creation  routine  prompts  the  user  for  the 
number  of  files  to  be  created.  The  user  is  then  prompted 
for  the  disk  drive  number,  i.e  0  —  3,  to  contain  the  files. 
The  data  files  on  the  selected  drive  are  displayed  to 
prevent  any  undesired  duplicate  file  naming.  The  user  is 
then  prompted  for  the  new  filenames.  If  a  previously  used 
filename  is  chosen,  the  user  is  alerted  that  a  duplicate 
filename  has  been  used  and  the  user  is  given  the  option  of 
renami ng  the  file  or  overwriting  the  existing  file. 
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After  the  files  are  named,  the  input  screen  is 
displayed.  The  screen  will  display  four  files  (columns) 
and  five  observations  (rows).  The  rows  and  columns  may  be 
scrolled  to  display  other  observations  as  needed. 

The  input  area  is  located  in  the  top-left  section  of 
the  screen.  Depressing  ENTER  transfers  the  input  area  to 
the  bracketed  destination  area. 


There  are  several  special  function  keys  which  must  be 
defined: 


ENTER  -  transfers  contents  of  the  input  area 

to  the  destination  area.  If  the  input 
area  is  blank,  the  destination  area  is 
moved  to  the  right. 

Arrows  -  moves  the  destination  area 


Shift  arrow 


moves  the  destination  area  to  its 
most  extreme  position 


D 


deletes  the  current  row 


I 


inserts  a  blank  row  immediately  before 
the  current  row 


CLEAR 


Returns  to  menu 
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writes  the  files  to  disk  and  returns  to 
the  menu 


The  data  is  written  to  disk  in  a  sequential  file  with 
the  first  field  containing  the  number  of  observations  and 
each  field  thereafter  containing  the  observations.  The 
routine  was  designed  for  maximum  user  convenience;  however. 
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these  convenience  -features  will  permit  undesired  "zero" 
valued  observations  if  the  user  is  not  aware  of  the  design 
of  the  program.  Some  precautions  and  guidelines  to  follow 


are: 


1.  All  displayed  files  will  be  assigned  a  common 

value  for  the  number  of  observations.  Therefore,  only 
files  which  will  have  the  same  number  of  observations 
should  be  entered  simutaneously. 


It  is  sometimes  tempting  to  "play”  with  the  arrow  keys 
to  move  the  destination  area  around  the  screen.  The 
number  of  observations  is  determined  by  the  maximum  row 
reached  by  the  destination  area;  therefore,  "playing" 
with  the  arrow  keys  may  result  in  an  unwanted  value  for 
the  number  of  observations. 


C.2.  EDIT  EXISTING  FILES  ROUTINE 


The  file  editing  routine  operates  similiarly  to  the 
new  file  creation  routine  with  only  minor  variation.  The 
user  is  prompted  for  the  disk  drive  containing  the  files 
and  then  prompted  for  the  selection  of  the  files  to  be 
edited.  The  files  are  displayed  on  the  screen  and  editing 
is  performed  exactly  as  described  in  the  new  file  creation 
routine.  Prior  to  writing  the  files  to  disk,  the  user  is 
given  the  opportunity  to  rename  the  file. 


C.3.  PURGE  FILES 


The  data  files  are  displayed  one  by  one  and  the  user 
is  asked  if  it  is  desired  to  erase  the  file  from  the  disk. 


6 


All  erasing  of  data  -files  should  be  performed  -from  the 
purge  routine  and  not  -from  the  operating  system  routines  as 
the  purge  routine  also  removes  the  -filename  from  the  data 
file  directory,  EDITDATA/DIR. 

C.4.  PRINT  DATA  FILES 

This  routine  permits  formatting  the  data  files  for 
printing.  The  user  is  prompted  for  the  files  to  be 
printed,  the  printer  line  width,  column  headers,  and  given 
the  option  of  numbering  the  observations.  The  maximum 
number  of  spaces  between  columns  is  calculated  and 
displayed.  The  user  is  prompted  for  the  desired  number  of 
spaces  between  columns.  Finally  the  user  is  asked  for  the 
report  title  and  the  files  are  printed. 


D.  PROGRAM  3,  CTABLE/BAS 


CTABLE/BAS  is  used  to  create  tables  used  in  the 
contingency  table  test  program.  The  program  operates 
similiarly  to  the  DBMGR/BAS  program  with  the  differences 
listed  below: 

1.  The  data  file  names  are  stored  in  CTABLE/DIR  versus 
EDITDATA/DIR. 

2.  In  the  new  file  creation  routine,  the  user  is 
prompted  for  the  number  of  rows  and  columns  in  the 
table. 

3.  The  ”Q"  and  "I"  keys  have  no  function. 
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E.  PROGRAM  4,  WILCOX/BAS 


WILCOX/BAS  performs  the  Wilcoxson  Signed  Rank  test 
using  data  stored  in  data  files  created  by  DBMGR/BAS.  The 
user  is  prompted  for  the  two  data  -files  -far  the  test.  The 
test  statistic  is  computed  and  the  user  is  prompted  for  the 
hypothesis  testing  problem,  i.e.  one  sided  or  two  sided. 

The  test  statistic,  sample  size  and  p-value  are  displayed 
and  the  user  is  prompted  for  the  significance  level  of  the 
test.  The  significance  level  is  compared  with  the  p-value 
and  acceptance  or  rejection  of  the  null  hypothesis  is 
determined. 

E.l.  WILCOXSON  SIGNED  RANK  TEST 

The  data  consists  of  two  paired  samples  of  n 
observations,  <  x  i ,  Y2*  ?  (;<n,  yn  >  - 

The  difference  dj  =  Xj— y^  is  computed  for  each  of 
the  n  pairs.  Pairs  with  a  difference  of  zero  are  omitted. 
The  absolute  differences,  id^t,  of  the  remaining  pairs 
are  ranked  with  rank  1  assigned  to  the  pair 
with  the  smallest  absolute  difference.  If  one  or  more 
pairs  have  the  same  absolute  difference,  the  average  rank 
is  assigned  to  the  pairs.  The  following  test  statistics 
are  computed: 

T+  =  Sum  of  ranks  with  positive  d^’s 

T~  =  Sum  of  ranks  with  negative  dx’s 


IS 


T  =  MIN  <T+, T  ) 

The  assumptions  -for  the  test  are: 

1.  The  distribution  o-f  each  is  symmetric 

2.  The  pairs  constitute  a  bivariate  random 

sample. 

The  hypothesis  may  be  stated  in  one  o-f  three  ways: 


A.  Hq  :  E  (X)  *  E  (Y) 

Hj  :  E  <  X )  s*  E  ( Y) 

B.  Hq  :  E  <  X )  i  E(Y) 

:  E  <X)  >  E(Y> 

C.  Hq  :  E  (X)  2.  E(Y) 

HD  :  E ( X )  <  E(Y) 

The  decision  rule  for  the  three  hypotheses  is  to  accept 
Hq  i*  : 

A.  P(Test  Statistic  <  T)  >  Significance  Level /2 

B.  PtTest  Statistic  <  T+)  Significance  Level 

C.  P(Test  Statistic  <  T_>  >  Significance  Level 


E.2.  DISTRIBUTION  THEORY  AND  ALGORITHM 


The  distribution  of  the  Wilcoxson  statistic,  T,  may  be 
found  by  randomizing  the  signed  ranks.  Since  there  are  two 
possible  signs,  +  and  — ,  and  n  ranks,  the  number  of 
combinations  in  the  randomization  is  2n.  The  enact 
cumulative  distribution  may  be  computed  by  examining  the 
2n  combinations  and  counting  those  with  the  sum  of  ranks 
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less  than  or  equal  to  the  test  statistic  and  dividing  by 
2n.  Such  a  method  is  obviously  slow  -for  all  but  very 
small  sample  sizes. 

A  number  of  approximations  to  the  exact  distribution 
have  been  proposed.  The  most  accurate  of  the  proposed 
approximations  is  the  Edgeworth  approximation  which  is 
documented  in  References  1,  2,  and  3.  This  approximation 
is  based  upon  the  fact  that  the  distribution  of  T  is 
asymptotically  normal  with  mean  and  variance: 

Uy  =  n (n+1 ) /4 

O^2  =  n(n+l)  (2n+l)/24 

The  normalized  Wilcoxson  Statistic,  z,  is  defined: 
z  =  (T  -  Uy )  /  <r^- 

The  Edgeworth  expansion  to  terms  of  size  l/n~  is 
P(T  <  t)  =  Fvz)  +  L4f(3)(z)/4!  + 

L6f  (5)  <z)  /4>'+35<L4)2f  (7)  /S! 

where  F(z)  is  the  cumulative  standard  normal  distribution 
and  : 

f*1^  is  the  i ’ th  derivative  of  the  normal  density 
function 

L4  =  -12C3n2+3n— 1) /5n<n+l) <2n+l> 

L6  »  576  (3n4+6n3-3n+l  > /7<n  (n-*-l  >  <2n+l)  >2 
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The  algorithm  implemented  in  this  package  uses  a  speed 
optimized  exact  distribution  when  the  execution  time  is 
small  and  the  Edgeworth  approximation  otherwise.  The 
decision  rule  is  : 


Approximation 


if  T  >  27  or  (N  >  9  and  T  >  N) 


Exact 


if  otherwise 


The  algorithm  must  use  the  normal  cumulative 
distribution.  The  algorithm  used  in  this  package  -for  the 
normal  cumulative  distribution  was  taken  from  the  thesis 
work  of  Mr.  R.P.  Isbell  and  is  documented  in  Reference  4. 


F.  PROGRAM  S,  MANN/BAS 


MANN/BAS  performs  the  Mann-Whit.ney  test  using  data 
stored  in  data  files  created  by  DBMGR/BAS.  The  user  is 
prompted  for  the  two  data  files  for  the  test.  The  test 
statistic  is  computed  and  the  user  is  prompted  for  the 
hypothesis  testing  problem,  i.e.  one-sided  or  two-sided. 
The  test  statistic,  sample  sizes  and  p-value  are  displayed 
and  the  user  is  prompted  for  the  significance  level  of  the 
test.  The  significance  level  is  compared  with  the  p-value 


and  acceptance  or  rejection  of  the  null  hypothesis  is 


determined. 


F.  1. 


MANN-WHITNEY  TEST 


The  data  consists  of  two  random  samples, 

- »*n  and  yi»Y2» *  * ‘ 'Vm*  The  twD 

samples  are  combined  and  ranked  -from  1  to  n+m  with  the 
smallest  observation  being  assigned  the  rank  o-f  1.  The 
test  statistics  are  : 

Ty  =  nm  +  n(n+l)/2  -  (Sum  of  ranks  for  sample  ;< ) 
Ty  =  nm  +  m(m+l)/2  -  (Sum  of  ranks  for  sample  y) 
T  =  MIN(Tx,Ty) 

The  assumptions  for  the  test  are  : 

1.  Both  samples  are  random  samples. 

2.  The  two  samples  are  independent. 

3.  If  there  is  a  difference  between  the  two 
population  distributions,  that  difference  is 
in  location  only,  i.e.  F,.  (u)  =  Fy(u-*-c> 
where  c  is  a  constant. 

The  hypothesis  may  be  stated  in  one  of  three  ways  : 


A. 

Ho  ! 

:  E  (  X  ) 

= 

E(Y) 

Hi  : 

:  E  (X) 

* 

E  (Y) 

B. 

Ho  5 

:  E  ( X ) 

E  (Y) 

Hi  : 

:  E  ( X ) 

E  (Y) 

C. 

Hq  : 

:  E  (  X  ) 

E  ( Y) 

Hi  : 

;  E  (  X  > 

< 

E  (Y) 

The  decision  rule  for  the  three  hypotheses  is  to  accept 


Ha  if  : 

A. 

P (Test 

Statistic  <  T) 

Significance 

Level /2 

B. 

P (Test 

Statistic  <  T., ) 

>  Significance 

Level 

C. 

P (Test 

Statistic  <  T, ) 

>  Significance 

Level 

no 


F.2. 


DISTRIBUTION  THEORY  AND  ALGORITHM 


The  distribution  of  the  Mann-Whitney  statistic  may  be 
found  by  randomizing  the  ranks  of  the  combined  data.  There 
are  (n+m)  ? / <n  • )  (m ! )  ways  in  which  the  ranks  may  be  grouped. 
The  exact  distribution  may  be  computed  by  examining  all 
possible  combinations  of  ranks  and  counting  those  with  the 
rank  sum  less  than  or  equal  to  nm+n(n+l)/2  minus  the  test 
statistic  and  dividing  by  (n+m) ! / (n ! ) (m ! > .  Such  a  method 
is  slow  for  all  but  very  small  sample  sizes. 

A  number  of  approximations  to  the  exact  distribution 
have  been  proposed.  The  most  accurate  of  those  proposed  is 
the  Edgeworth  approximation  documented  in  Reference  5.  The 
approximation  is  based  on  the  fact  that  the  test  statistic 
is  asymtot i cal 1 v  normal  with  mean  and  variance  : 

Uj.  =  nm/2 

=  nin(n+m+l)  /12 

The  normalized  Mann— Whitney  Statistic,  z,  is  defined  : 

z  =  C  T  -  U£  )  /  9^ 

The  Edgeworth  expansion  is  : 

P(t<t)  =  F(z)  +L4f (3) <z> /4!  +  L6f <5)  <z) /6\  + 

35(L4)2f (7) <z) /8! 

where  F<z)  is  the  standard  normal  cumulative  distribution 


is  the  i  ’  th  derivative  o-f  the  standard  normal 
densi ty 

=  -  (m“+n^+n+ffl++n(n)  /  < 20nm  (n+m+1 )  ) 

LS  =  NUM  /  DENOM 
where  : 

NUM  —  2 (n^+m^) +4nm <n“+m2) +6n2m2 

+4  (n +7nm  (n+m)  +n*‘+m“+2nm-n-m 
DENOM  =  210n2m2<n+m+l)2 

The  algorithm  that  is  implemented  in  this  package  uses 
a  speed  optimized  exact  distribution  -for  small  sample  sizes 
and  the  Edgeworth  approximation  -for  1  arger  sample  sizes. 

The  decision  rule  is  : 

appro:-:  i  mat  i  an  i -f  m  >  9  or  T  >  27  -  n(n+l)/2 

exact  if  otherwise 

where  n  and  m  are  the  smaller  and 
larger  sample  sizes  respectively 

G.  PROGRAM  6,  SMIRNOV /BAS 

SMIRM0V/BA3  performs  the  Smirnov  test  on  data  stored 
in  data  files  by  DBMGR/BAS.  The  user  is  prompted  for  the 
two  data  files  for  the  test.  The  test  statistic  is 
computed  and  the  user  is  prompted  for  the  hypothesis 
testing  problem,  i.e.  one-sided  or  two-sided.  The  test 
statistic,  sample  sizes,  and  p-value  are  displayed.  The 


user  is  prompted  for  the  significance  level  of  the  test. 


The  signif icance  level  of  the  test  is  compared  to  the 
p— value  and  acceptance  or  rejection  of  the  null  hypothesis 
is  determined. 

S. 1.  SMIRNOV  TEST 

The  data  consists  of  two  independent  random  samples 
with  unknown  cumulative  distribution  functions  Fiji)  and 
G<>:).  S^(::)  and  S-,(x)  are  their  empirical  cumulative 

distribution  functions.  The  test  statistics  are  : 

D*  =  MA X  ( S ^  ( :•: )  -  S2  (:•;>> 

D“  =  MAX<S2<::>  ~ 

D  =  MAX(D+,D~) 

The  assumptions  for  the  test  are  : 

1.  The  samples  are  random  samples 

2.  The  two  samples  are  independent 

3.  The  random  samples  are  continuous 
The  hypotheses  may  be  stated  : 


A . 

u 

■  3 

: 

F  ( :: ) 

= 

5  < ::  ) 

Mi 

• 

F  ) 

A 

G  < > 

B- 

Ho 

: 

F  < ) 

G  <:: ) 

H1 

■ 

F  ) 

GCi) 

C. 

Ha 

: 

F  (:; ) 

G  ■  ■: ) 

Hi 

: 

F  < ::  5 

G  i ::  ) 

The  decision  rule  for  the  three  hypothese  is  to  accept 
Hq  if  : 

A.  PCTest  Statistic  L  D)  Significance  Level/2 

3.  F'lTest  Statistic  v  D+) 


Significance  Level 


•4 


C.  PtTest  Statistic  <  D  )  >  Significance  Level 

G. 2.  DISTRIBUTION  THEORY 

The  exact  distribution  o-f  the  test  statistic  is 
difficult  to  analyze  for  the  general  case.  E.F.  Drion,  in 
Reference  6,  derived  the  distribution  for  the  special  case 
of  equal  sample  sizes.  It  is  known  that  the  test  statistic 
has  the  asymptotic  distribution  : 

G(S)  =  1  -  ^  (-1 ) i-1exp <-2i2S2> 

i 

where  GCS)  is  known  as  the  Smirnov  distribution 

S  =  D  SQR(  (n+m)  /nnt)  where  n  and  m  are  sample  sizes 

It  is  generally  agreed  that  sample  sizes  of  100  or 
more  are  required  cor  the  overestimating  effects  of  the 
asymptotic  distribution  to  become  negligible.  The  problem 
is  then  to  improve  upon  the  Smirnov  distribution  for  sample 
sizes  of  less  than  100.  P.J.  Kim,  in  Reference  7, 
presented  an  improvement  which  reduced  the  over est i mat l ng 
effects  for  the  smaller  sample  sizes.  The  improvements 
presented  by  P.J.  Kim  and  implemented  in  this  package  are  : 

3  =  D  •  (SDR  (rid/  (n+m)  )  )  +2/  (3  •  (SQR  (m)  )  )  for  m=kn,  k=l,2, . 

S  =  < D- 1 / 2m ) •  (  SQR ( n ) )  f  or  n / m  <  0 .  1 

3  =  D  •  (SQR  (nm/  (n+in)  !  )  ■+•  2/  (5  ♦  (SQR  (n )  )  )  otherwise 
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H.  PROGRAM  7,  NL ILL /BAS 


NLILL/BAS  performs  a  test  for  normal  distribution  on  a 
data  set  stored  in  a  data  file  created  by  DBMGR/BAS.  The 
test  may  be  performed  for  the  three  cases  : 

1.  mean  and  variance  unknown 

2.  mean  known  and  variance  unknown 

3.  mean  and  variance  known 

The  user  is  prompted  for  the  data  file  to  be  tested 
and  prompted  for  the  specific  case  to  be  tested.  In  the 
cases  where  a  parameter  is  known,  the  user  is  prompted  for 
the  known  parameter.  The  test  statistic,  sample  sice  and 
p— value  are  displayed.  The  user  is  prompted  for  the 
significance  level  of  the  test.  The  significance  level  is 
compared  with  the  p-value  and  acceptance  or  rejection  of 
the  null  hypothesis  is  determined. 

H.l.  TEST  FOR  NORMAL  DISTRIBUTION 

The  data  consists  of  a  random  sample  x  ;<n 

from  an  unknown  cumulative  distribution  Fix).  G(::)  is  the 
normal  cumulative  distribution  with  parameters  either 
specified  or  estimated  from  the  sample,  depending  upon  the 
specific  case.  The  test  statistics  are  : 

D*  =  MAX (i /n  -  )  ) 

D~  =  MAX  ( G  t )  -  <i-l)/n>> 

D  =  MAX  CD*, D~) 


The  assumption  -for  the  test  is  that  the  sample  is  a 


random  sample.  The  hypotheses  are  : 


Case  1.  Hq  :  F(x)  is  normal  with  unspecified 
mean  and  variance 
Hj  :  F<x)  is  not  normal 

Case  2.  HQ  :  F ( >; )  is  normal  with  specified  mean 
:  F(x)  is  not  normal  with  specified 
mean 


Case  3.  HQ  :  F(>:)  is  normal  with  specified  mean 
and  variance 

Hj  :  F<x)  is  not  normal  with  specified 
mean  and  variance 


The  decision  rule  is  accept  HQ  if  : 

P(Test  Statistic  <  D)  >  Significance  level 


H.2.  DISTRIBUTION  THEORY  AND  ALGORITHM 


The  distribution  of  the  test  statistic  differs  for 
each  of  the  three  cases.  For  Case  I,  where  the  mean  and 
variance  are  unknown,  the  distribution  was  approximated 
using  Monte  Carlo  techniques  by  H.W.  Lilliefors  and 
published  in  Reference  8.  The  distribution  is  presentd  in 
tabular  form  and  is  not  described  in  an  analytical  sense. 

In  Reference  8,  Lilliefors  recognized  that,  for  a 
particular  significance  level,  the  distribution  of  the  test 
statistic  differed  by  a  multiplicative  constant  from  that 
of  the  Kolmogorov  distribution.  This  package  capitalizes 
on  this  concept;  however,  the  multiplicative  value  is  a 
function  of  the  value  of  the  cumulative  distribution.  The 
solution  must  therefore  be  solved  iteratively.  The 
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algorithm  -follows: 

T  < -  Test  Statistic 

Denominator  < -  1 

Do  Until  ABS (Tdenom  -  Denominator)  <  .0001 
D  < -  T/Denominator 

S  < -  D*(-.01  +  SQR(N)  +  0. S5/SQR (N) ) 

Tdenom  < -  Denominator 

CDF  < -  1  -  <.135)S 

If  CDF  <  .85 

Then  Denominator  < -  1  —  . 24xCDF 

I-f  .85  <  CDF  <  .90 

Then  Denominator  < -  1.337  —  . 66xCDF 

I-f  .90  <  CDF  <  .95 

Then  Denominator  < -  1.339  —  . 64xCDF 

If  .95  <  CDF  <  .975 

Then  Denominator  < -  1.775  —  1.12xCDF 

If  CDF  >  .975 

Then  Denominator  < -  2.133  —  1.467xCDF 

End  Do 


For  Case  2,  where  the  mean  is  known  and  the  variance 
is  unknown,  the  distribution  was  approximated  by  M.A. 
Stephens  and  published  in  Reference  9.  A  similiar  solution 
was  found  to  exist  for  this  case  as  existed  for  Case  1, 
i.e.  the  cumulative  distribution  must  be  calculated 
iteratively.  The  algorithm  for  the  iterative  solution  is  : 

T< - Test  Stastistic 

Denominator  < -  1 

Do  Until  ABS  ( Tdenom  -  Denominator)  <  .0001 
D  < -  T/Denominator 

S  < -  D*(.35  +  SQR(N)  +  . 283/SQR(N)) 

Tdenom  < -  Denominator 

CDF  < -  1  -  (.135) 5 

If  CDF  <  .9 

Then  Denominator  < -  1.33  -  -QxCDF 

If  .9  <  CDF  <  .95 

Then  Denominator  < -  1.47  -  -4xCDF 

If  .95  <  CDF  <  .975 

Then  Denominator  < -  1.85  —  .8xCDF 

If  CDF  >  .975 

Then  Denominator  < -  1.07 


End  Do 


For  Case  3,  where  the  mean  and  variance  are  completely 
specified,  the  distribution  is  the  Kolmogorov  distribution. 
The  Kolmogorov  distribution  is  a  special  case  of  the 
Smirnov  distribution  with  one  sample  size  approaching 
infinity  and  has  the  asymptotic  distribution  : 

G (S)  =1-2  <-l> 1-1EXP(-2i2S2> 

! 

with  S  =  D  *  SDR  <N> 

M.A.  Stephens  proposed  improvements  to  the  asymptotic 
distribution  in  Reference  9  which  improve  the  accuracy  of 
the  asymptotic  distribution  for  small  sample  sizes.  This 
package  implements  the  improvements  proposed  by  Stephens  in 
Reference  9,  i.e.  : 

3  =  D  •  ( .  12  +  SQR(n)  +  .U/SQR(n)> 

I.  PROGRAM  8,  EL ILL /BAS 

EL  ILL/BAS  performs  a  test  for  exponential  distribution 
on  data  stored  in  data  files  by  DBMGR/BAS.  The  test  may  be 
performed  for  the  cases  where  the  parameter  is  known  and 
where  the  parameter  is  unknown. 

The  user  is  prompted  for  the  data  file  to  be  tested 
and  for  the  specific  case,  i.e.  parameter  known  or  unknown. 
In  the  case  where  the  parameter  is  known,  the  user  is 
prompted  for  the  parameter  value.  The  test  statistic, 


sample  size  and  p- value  are  displayed.  The  user  is 
prompted  -for  the  significance  level  of  the  test.  The 
significance  level  is  compared  with  the  p-value  and 
acceptance  or  rejection  of  the  null  hypothesis  is 
determined. 

1.1.  TEST  FOR  EXPONENTIAL  DISTRIBUTION 

The  data  consists  of  a  random  sample  x  >:n 

from  some  unknown  cumulative  distribution  F<x>.  G(x>  is 
the  exponential  cumulative  distribution  function  with  the 
parameter  either  estimated  from  the  data  or  specified, 
depending  upon  the  specifiic  case.  The  test  statistics 
are: 

D+  =  MAX (i/n  -  G(x) ) 

D~  =  MAX (G (x )  -  (i-l)/n) 

D  =  MAX(D+,D“> 

The  assumption  for  the  test  is  that  the  sample  is  a 
random  sample.  The  hypotheses  are: 


Case  1 . 

Ho  s 

F  <  x  ) 

i  s 

exponential  —  parameter 

specified 

Hi  : 

F  <  x ) 

i  s 

not  exponential 

Case  2. 

Ho  : 

F  ( x ) 

i  s 

exponential  -  parameter 

unspeci f i ed 

»l  * 

F  <x  ) 

i  s 

not  exponential 

The  decision  rule  is  to  accept  H0  if  : 

P(Test  Statistic  <  D)  >  Significance  Level 
1.2.  DISTRIBUTION  THEORY  AND  ALGORITHM 

The  distribution  of  the  test  statistic  differs  for  the 
two  cases.  For  the  case  where  the  parameter  is  unknown. 
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the  distribution  was  approximated  using  Monte  Carlo 
techniques  by  H.W.  Lilliefors  and  published  in  Reference 
10.  M.A.  Stephens,  in  Reference  9,  expanded  the  work  of 
Liffiefors  by  providing  improvements  for  the  small  sample 
cases.  The  algorithm  is  an  iterative  method  which  assumes 
a  solution  and  computes  the  test  statistic  which  would 
yield  the  assumed  value  for  the  cumulative  distribution 
function  (CDF).  The  assumed  CDF  is  updated  until  the 
difference  between  the  calculated  test  statistic  and  the 
actual  test  statistic  is  sufficiently  small.  The  algorithm 
i  5  : 


T  < -  TEST  STATISTIC 

DENOMINATOR  < -  .26  +  SQR(n)  +  .5/SQR(N) 

FI  < - 0 

F2  < - 1 

DO  UNTIL  ERROR  <  .0001 


CDF  < - 

■  (F2 

- 

FI)  /2 

IF 

CDF  > 

.95 

THEN 

TS  <- 

— 

5. 1 25x  CDF 

- 

3. 7S08 

IF 

.90  < 

CDF 

.95 

THEN 

TS  <- 

— 

2. 02xCDF 

8367 

IF 

.80  < 

CDF 

< 

.90 

THEN 

TS  <- 

— 

1. 124xCDF 

.0249 

IF 

.70  < 

CDF 

.80 

THEN 

TS  <- 

— 

0. 743xCDF 

-f 

.  2799 

IF 

.50  < 

CDF 

.70 

THEN 

TS  <- 

— 

0. 546:tCDF 

■f 

.4187 

IF 

•  jO  ( 

CDF 

.  50 

THEN 

TS  <- 

— 

0. 464xCDF 

.4588 

IF 

.20  < 

CDF 

< 

V 

■  oO 

THEN 

TS  <- 

— 

0. 435x CDF 

.4525 

IF 

.  10  < 

CDF 

< 

.20 

THEN 

TS  <- 

— 

0.  587 :<  CDF 

.  4321 

IF 

CDF  < 

.  10 

THEN 

TS  <- 

— 

0.  346:;  CDF 

.4062 

TS  < -  TS /DENOMINATOR  +  . 2/N 

ERROR  < -  T-TS 

IF  ERROR  >  0  THEN  FI  < - 

END  DO 


CDF  ELSE  F2 


CDF 


For  Case  2,  with  the  parameter  known,  the  distribution 
is  the  Kolmogorov  di stribution.  The  distribution  and  the 
algorithm  used  in  this  package  is  described  in  the  -fully 
specified  case  of  the  test  for  normal  distribution. 

J.  PROGRAM  9,  CHISQU/BAS 

CHISQU/BAS  performs  the  r  ;<  c  contingency  table  test 
using  data  stored  in  data  files  created  by  CTAELE/BAS.  The 
user  is  prompted  for  the  data  file  to  be  tested.  The  test 
statistic,  degrees  of  freedom,  fraction  of  cells  with 
expected  frequency  less  than  5,  number  of  cells  with 
expected  frequency  less  than  1,  and  the  p— value  are 
printed.  The  user  is  prompted  for  the  significance  level 
of  the  test.  The  significance  level  is  compared  with  the 
p-value  and  acceptance  or  rejection  of  the  null  hypothesis 
is  determined. 

J.l.  CONTINGENCY  TABLE  TEST 

The  contingency  table  test  may  be  used  to  present  a 
tabulation  of  data  contained  in  several  samples  to  test  the 
hypothesis  that  the  probab i 1 i ti es  do  not  differ  from  sample 
to  sample.  Another  use  of  the  r  ::  c  contingency  table  is 
with  a  single  sample  where  each  element  in  the  sample  may 
be  classified  into  one  of  r  different  categories  according 
to  one  category  and  into  one  of  c  different  categories 


according  to  a  second  category. 


The  two  applications  are 


treated  the  same  in  the  statistical  analysis,  but 
basic  diff erences  between  the  two  applications  justi-fy 
separate  discussions. 

J.1.1  TEST  FOR  DIFFERENCES  IN  CELL  PROBABILITIES 

The  data  consists  of  random  samples  drawn  -from  r 
populations.  Each  observation  in  each  sample  may  be 
classified  into  one  o-f  c  different  categories.  Let  0^  ^ 
be  the  number  of  observations  from  the  i ' th  sample  that 
fall  into  category  j.  Let  n^  be  the  number  of 
observations  from  the  i  ’  th  population,  c^  be  the  number 
of  observations  in  the  j ' th  category,  and  N  be  the  total 
number  of  observations.  The  test  statistic  is  : 

T  —  21  2(0^  -  E.  j)2/Eij 
where  Ei j  =  n^c^/N 

The  assumptions  are  that  each  sample  is  a  random 
sample,  the  samples  are  mutually  independent,  and  each 
observation  may  be  categorised  into  exactly  one  of  c 
categories.  The  hypotheses  may  be  stated  : 

Hq  :  All  of  the  probabilities  in  the  same  column 
are  equal 

Hj  :  At  least  two  of  the  probabi 1 i ti es  in  the  same 
column  are  not  equal 

The  decision  rule  is  to  accept  the  null  hypothesis  if  : 

P(Test  Statistic  <  T)  >  significance  level 
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J.1.2  TEST  FOR  INDEPENDENCE 


The  data  consists  of  one  random  sample  o-f  size  N.  The 
data  may  be  classified  by  one  o-f  two  criteria.  Using  the 
first  criteria,  each  observation  is  associated  with  one  of 
r  rows  and  using  the  second  criteria,  each  observatin  is 
associatd  with  one  of  the  c  columns.  Let  0^  j  be  the 
number  of  observations  in  row  i  and  column  j,  n^  be  the 
number  of  observations  in  the  i ’ th  row,  and  Cj  be  the 

number  of  observations  in  the  j’th  column.  The  test 

statistic  is  as  defined  in  the  test  for  differences  in  cell 
probabi 1 i ties. 

The  hypothesis  may  be  stated  : 

H0  :  P(row  i,  column  j>  =  P(row  i)  x  P (column  j ) 

:  P(row  i,  column  j)  =  P(row  i)  x  P  (column  j) 

The  decision  rule  is  to  accept  the  null  hypothesis  if  : 

PCTest  Statistic  <  T)  >  Significance  Level 
J.2  DISTRIBUTION  THEORY  AND  ALGORITHM 

The  exact  distribution  of  T  is  difficult  to  calculate 
for  all  but  the  2  ::  2  case.  It  is  well  known  that  T  has 
asymptotically  a  Chi-Square  distribution.  The  Chi-Square 
approximation  is  valid  for  large  exec ted  cell 
probabilities,  i.e.  E^ , 


The  approximation  is 


considered  to  be  satisfactory  i-f  no  Ei  j  is  less  than  1 
and  not  more  than  2.0'/.  of  the  j '  s  are  less  than  5. 

The  algorithm  used  in  this  package  for  the  Chi— Square 
distribution  was  extracted  from  the  thesis  of  Mr.  R.P. 
Isbell.  The  details  of  his  approximation  to  the 
distribution  are  detailed  in  Reference  1.  For  the  2x2 
case,  a  correction  factor  proposed  by  F.  Yates  in  Reference 
11  is  used  to  improve  the  Chi-Square  approximation  to  the 
exact  distribution. 

K.  PROGRAM  10,  NPSTAT/BAS 

NPSTAT/BAS  contains  seven  non-parametr i c  probability 
distributions,  the  normal  distribution,  and  their  inverse 
probability  distributions.  The  distributions  approximated 
are  : 

1.  Normal  Distribution 

2.  Wilcoxson  Signed  Rank  Distribution 

3.  Mann— Whitney  Di str i but i on 

4.  Smirnov  Distribution 

5.  Kolmogorov  Distribution 

6.  Lilliefors  Test  for  Normal  Distribution 

7.  Lilliefors  Test  for  Exponential  Distribution 

The  user  is  provided  a  menu  to  select  the  distribution 
desired.  The  user  is  prompted  for  the  required  inputs  and 
the  cumulative  distribution  is  displayed. 

The  approx i mati ons  to  the  distributions  are  described 
in  detail  in  the  hypothesis  test  programs. 
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TABLE  2 


MANN-WHITNEY  DISTRIBUTION  ACCURACY  COMPARISON 


Samp  1 e 
Size 

Sample 

Size 

Test 

Statistic 

Exact 

Approx 

6 

h 

7 

.  0465 

.  U-463 

6 

6 

10 

.  1201 

.  1200 

6 

6 

12 

.  1970 

in 

4} 

■*-< 

6 

3 

1 1 

.  053? 

.  0542 

6 

3 

14 

.1142 

.1143 

6 

a 

16 

.  1725 

.  1724 

3 

g 

16 

.  0524 

.  0524 

3 

a 

20 

.1172 

1  173 

a 

a 

.  1641 

.  1641 

B 

10 

21 

.  0506 

.  0506 

5 

1 0 

25 

.  1 0  1 5 

0 

f'i 

a 

10 

23 

.  1577 

'  t=,~T'7 
a  X.  / 

10 

io 

23 

.  0526 

a  0526 

i  o 

10 

.  1038 

.  i  :>ae 

10 

10 

36 

•  1575 

.  1575 

TABLE  2 


MANN-WHITNEY  DISTRIBUTION  ACCURACY  COMPARISON 


Sample 

Size 

Test 

Statistic 

E;:  act 

Appro;: 

b 

7 

.  0465 

.  0)463 

6 

10 

.  1201 

.  1  200 

a 

12 

„  1970 

.  1965 

3 

u 

.  0539 

.  0542 

3 

14 

.  1 142 

.1143 

a 

1 6 

.  1725 

.  1724 

s 

16 

.  0524 

.  0524 

a 

20 

.1172 

.1173 

a 

.  1641 

1 

.  1641 

10 

21 

.  0506 

.  0506 

1  0 

25 

.1015 

.  1 0  1 6 

10 

23 

.  1577 

.  1 577 

10 

23 

.  0526 

.  0526 

10 

o 

.  1088 

.  i  ose 

1 0 

36 

.  1 575 

.  1575 

TABLE  3 


SMIRNOV  DISTRIBUTION  ACCURACY  COMPARISON 


Probablity  =  0.95 


m  = 

25 

m  =  50 

m  =  lOO 

n/m 

Exact 

Approx 

Ex  act 

Approx 

Ex  act 

Approx 

.  1 

.5209 

.5185 

.3797 

.3737 

.2 

.5321 

.  5339 

.  3903 

.3910 

.2830 

.2833 

•  3 

.  3393 

.  3434 

.2443 

.2463 

.4 

.4188 

.4277 

.3044 

■  3086 

.2192 

.2213 

.5 

.2769 

.2765 

.2003 

.2003 

.6 

.3704 

.  3733 

.2680 

.2694 

.  1925 

.  1932 

.7 

.2579 

.2571 

.  1848 

.  1844 

.3 

.3428 

.3429 

.2474 

.2475 

.  1774 

.  1775 

.9 

.2409 

.2397 

.  1726 

.  1719 

1.0 

.  3087 

.3082 

.2256 

.2258 

.  1634 

.  1635 

Probablity  =  0.975 


m  = 

25 

m  = 

50 

m  = 

lOO 

n/m 

Exact 

Approx 

Ex  act 

Approx 

Exact 

Approx 

.  1 

.  5775 

.4224 

.4142 

.5910 

.5996 

.4348 

.3152 

■  ^IdI 

m  Ct 

.3777 

.  2719 

.  2742 

.4 

.  4764 

.4779 

•  3390 

.3441 

.2440 

.2464 

.5 

.  3096 

.  3094 

.  2236 

.2235 

,  6 

.  4126 

.4172 

.2985 

.3004 

.2143 

.2151 

.7 

.2871 

.  2867 

.2056 

.  2053 

.8 

.3818 

.  3832 

.  2754 

.2759 

19744 

.  1976 

.9 

.  2682 

.2673 

.  1920 

.  1914 

1.0 

.3464 

.3462 

.  2524 

.  2526 

.  1822 

.  1825 

TABLE  4 


KOLMOGOROV  DISTRIBUTION  ACCURACY  COMPARISON 


Samp 1 e 

1! 

a. 

0.90 

P  = 

0.95 

P  = 

0.975 

Size 

Exact 

Approx 

Exact 

Approx 

Exact 

Approx 

5 

.447 

.446 

.509 

.509 

.  563 

.564 

10 

.  323 

.323 

.369 

.369 

.409 

.409 

15 

.266 

.267 

.304 

.304 

.  338 

.  338 

20 

.232 

.232 

.265 

.265 

.294 

.294 

25 

.208 

.209 

.238 

.238 

.264 

.264 

30 

.  194 

.  191 

.218 

.218 

•  242 

.242 

35 

.  177 

.  177 

.202 

.202 

.  224 

.  224 

40 

.  165 

.166 

_ 

.  189 

.  189 

.210 

.210 

I 
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TABLE  5 


LILLIEFORS  NORMALITY  TEST  DISTRIBUTION 
ACCURACY  COMPARISON 


Sample 

Size 

P  =  < 

).  90 

.95 

P  =  0.99 

Exact 

Approx 

Exact 

Approx 

Exact 

Approx 

5 

.  315 

.314 

.337 

.343 

.405 

.396 

10 

.  239 

.239 

.  258 

.261 

.294 

.  302 

15 

.201 

.  200 

.  220 

.219 

.257 

20 

.  174 

.  176 

.  190 

.  192 

.231 

25 

.  15S 

.  159 

.  1 73 

.  200 

.  200 

30 

.  144 

.  146 

.  159 

.  187 

.  184 

TABLE  6 


LILLIEFORS  EXPONENTI ALITY  TEST  DISTRIBUTION 
ACCURACY  COMPARISON 


Samp 1 e 
Size 

o 

tl 

a. 

.80 

■8 

.90 

P  =  0 

.95 

Exact 

Approx 

Exact 

Approx 

Exact 

Approx 

5 

.  3603 

.3615 

.  4045 

.  4028 

.4420 

.4400 

10 

.2626 

.2642 

.  2955 

.2956 

.  3244 

.  3239 

15 

.  2174 

.2185 

.2448 

.2448 

.  2690 

.  2636 

20 

.  1393 

.  1 905 

.2345 

.2346 

25 

.  1703 

.1711 

.  1918 

.2110 

.  21 10 

30 

.  1539 

.  1567 

.  1756 

.  1760 

.  1932 

.  1933 

35 

.  1447 

.  1454 

.  1630 

.  1 63w 

.  1793 

.  1795 

40 

.  1356 

.  1362 

.  1528 

.  1531 

.  1681 

.  1683 

45 

.1281 

.  1286 

.  1443 

.  1 445 

.  1583 

.  1589 

50 

.  1217 

.1221 

.1371 

a  1  vj>7  O 

.  1509 

.  1510 

60 

.1118 

.1117 

.  1255 

.  1256 

.  1381 

.  1381 

70 

.  1033 

a  1  Ow5 

.1164 

.1164 

.  1281 

.  1231 

30 

.  0968 

.  0969 

.  1090 

.  1 09 1 

.  1 200 

.  1200 

90 

.  09 1 4 

.0914 

-  1029 

.  1029 

.1132 

.  1 132 

100 

.  0868 

.  0868 

.  0977 

.  0977 

.  1  075 

m 

PROGRAM  1 


MENU/BAS 


10  CLEAR  350 

20  DEFINT  I-K  s  DEFSTR  Z 

30  LB*=CHR*(123)  s  RB*=CHR* ( 1 25 )  :  MB*=CHR* ( 176)  : 

BL*=CHR* (252) 

40  CLS: PRINT  TAB (10)  "NON-PAR AMETR I C  STATISTICAL  ANALYSIS 
SOFTWARE" 

50  PRINT  TAB (31)  "by" 

60  PRINT  TAB (21)  "Robert  Lee  Zangmei ster " 

70  PRINT  s  PRINT  TAB(13)  "Submitted  in  partial  fulfillment 
of  the"  :  PRINT  TAB(17)  "requirements  for  the  degree  of" 

SO  PRINT  :  PRINT  TAB (19)  "MASTER  OF  APPLIED  SCIENCE" 

90  PRINT  :  PRINT  TAB (28)  "from  the” 

100  PRINT  :  PRINT  TAB (19)  "NAVAL  POSTGRADUATE  SCHOOL" 

110  PRINT  TAB (25)  "December  1982" 

120  PRINTsPRINT  TAB(17)  "Press  SPACE  BAR  to  continue" 

130  ZI=INKEY*  :  IF  ZI  =  "“  THEN  130 

140  CLS  :  PRINT  TAB (20) ; "PROGRAM  SELECTION  MENU" 

150  PRINT  3128, "(1)  Data  Base  Manager" 

160  PRINT  "(2)  Probability  Distributions" 

170  PRINT  ” (3)  Hypothesis  Tests" 

180  KP=896  :  KS=6  :  GOSUB  500 

190  ON  10  GOTO  200,250,1000 

200  CLS: PRINT  "DATA  BASE  OPTION" : PRINT 

210  PRINT  "(1)  Column  Format  Data  Entry" 

220  PRINT  "(2)  Table  Format  Data  Entry" 

230  KP=896:KS=2: GOSUB  500 

240  IF  IQ  =  1  THEN  RUN  "DBMGR/BAS"  ELSE  RUN  "CT ABLE/ BAS" 
250  RUN  " NPST AT / BAS " 

500  REM  *  *  *  KEYBOARD  ENTRY  ROUTINE  *  * 

510  PRINT  3KP+5, "OPTION  DESIRED  ” ;  LB*;  MB*;  RB*,- 
520  ZI=INKEY*  :  IF  ZI=“"  THEN  520ELSE  IO=VALCZI) 

530  PRINT  3KP+21,  ZIj:  FOR  1=1  TO  100  :  NEXT 

540  IF  I0<1  OR  I  OKS  THEN  PRINT  3KP,  BL*;  :  PRINT  3KP,  "** 

NOT  A  VALID  OPTION  ***";:  FOR  1=1  TO  lOOO  :  NEXT  :  PRINT 

3KP,  BL*;  :  GOTO  510 

550  RETURN 

1000  CLS: PRINT  TAB (20) ; "HYPOTHESIS  TEST  MENU " : PR I NT 
1010  PRINT  3128, "(1)  Wilcoxsan  Signed  Ranks  Test" 

1020  PRINT  "(2)  Mann-Whitney  Test" 

1030  PRINT  "(3)  Chi-Square  Contingency  Table  Test" 

1040  PRINT  " (4>  Smirnov  Test" 

1050  PRINT  "(5)  Test  for  Normality" 

1060  PRINT  "(6)  Test  for  Exponential i ty" 

1070  KP=896:KS=6: GOSUB  500 

1080  ON  10  GOTO  1100,1110,1120,1130,1140,1150 
1100  RUN  "WILCOX/BAS" 

1110  RUN  "MANN/BAS" 

1120  RUN  "CHISQU/BAS" 

1130  RUN  "SMIRNOV /BAS" 

1140  RUN  "NL  ILL /BAS" 

1150  RUN  "ELILL/BAS" 
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PROGRAM  2 


DBMGR/BAS 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFSTR  2 

40  DIM  F*(95) ,F1*(11) ,M*(11) 

50  LB*=CHR*(123) : RB*=CHR* ( 125) :MB*»CHR* < 176) : BL*=CHR* (252) 
60  CLS 

70  PRINT  TAB (20) ; "DATA  FILE  MANAGEMENT  PROGRAM" 

80  PRINT! PRINT: PRINT  "MENU:" 

90  PRINT  "(1)  Create  New  Data  Files" 

10O  PRINT  "(2)  Edit  Existing  Files" 

HO  PRINT  "(3)  Purge  Existing  Files” 

120  PRINT  "(4)  Print  Data  Files” 

130  PRINT  "(5)  Return  to  Master  Menu" 

140  KP=896 :  KS=5 :  GOSUB  180 

150  ON  10  GOSUB  260,310,330,340,170 

160  GOTO  60 

170  RUN  "MENU/BAS” 

180  REM  *  *  KEYBOARD  INPUT  SUBROUTINE  * 

190  PRINT  3KP+5, "OPTION  DESIRED  “ ; LB*; MB*; RB*; 

200  21= INKEY*: IF  21=”"  THEN  200  ELSE  ID=VAL(2I) 

210  PRINT  3>KP+21 ,21;:  FOR  1  =  1  TO  100:  NEXT 

220  IF  IOC1  OR  IO>KS  THEN  PRINT  3)KP ,  BL* ;  :  PR  I  NT  5>KP,"** 

NOT  A  VALID  OPTION  ***";:  FOR  1=1  TO  1000:  NEXT:  PRINT 
3KP,  BL* ; : GOTO  190 
230  RETURN 

240  DN*= INKEY*: IF  DN*=""  THEN  240  ELSE  I=VAL(DN*) 

250  IF  I<0  OR  I >3  THEN  240  ELSE  RETURN 

260  *  SUB  TO  PERFORM  NEW  DATA  FILE  CREATION  ROUTINES 

270  CLS  :  INPUT  "How  Many  Variables  are  to  be  Used 

(1-11)  " ;  NF :  CT =NF :  RT  =  1 

280  IF  NF< 1  OR  NF>11  THEN  270 

290  GOSUB  14000  :  GOSUB  12000  :  N=1  :  CR=1  :  GOSUB  6020 

300  CLS: GOSUB  15050: RUN 

310  GOSUB  350  :  CR=1  :  GOSUB  6020 

320  GOSUB  15000  :  RUN 

330  GOSUB  5000  :  RUN 

340  GOSUB  350  :  GOSUB  4000  :  RUN 

350  REM  *  *  *  READ  MULTIPLE  ARRAYS  *  * 

360  K=0  :  DL*=" / " 

370  CLS:PRINT  "Disk  Drive  Containing  Files  (0—3)?’* 

380  GOSUB  240 

390  SL= 1715: GOSUB  1200: IF  SL=0  THEN  GOTO  370 
400  GOSUB  1400  :  PRINT  "" 

410  PRINT: IF  10  =  2  PRINT  "Enter  the  Files  To  Be  Edited 
(e.g.  1/3/4)";  ELSE  PRINT  "Enter  the  Files  To  Be  Printed"; 
420  INPUT  B*  : P1=0 
430  P2= INSTR ( 1 , B* , DL*  > 
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440  IF  P2=0  THEN  NV=VAL<B«) :  Pl=l  :  GOTO  460 

450  NV=VAL(LEFT*(B*,P2-1> )  s  B*=RIGHT* <B*, LEN (B*) -P2) 

460  IF  NV>FO  THEN  PRINT: PRINT"  *  *  *  INVALID  FILE  #  * 

* " : GOTO  410 

470  K=K+1  :  IF  K=12  THEN  NF=11  :  GOTO  530 

480  F1*(K)=F*<NV) 

490  IF  P1=0  THEN  430 

500  B*=" PRINT : PRINT  “Do  You  Wish  To  Read  Files  From 
Another  Disk (Y/N) ? 

510  ZI=INKEY*:IF  ZI  =  “"  GOTO  510  ELSE  IF  ZI="Y“  PRINT 
Z I : GOTO  370 

520  NF=K:CT=K: PRINT  "N" 

530  MM=— 1 

540  FOR  1=1  TO  NF 

550  OPEN  “I  ,1,F1*(I)  :  INPUT  41, M  :  CLOSE  1 
560  IF  M  >  MM  THEN  MM=M 
570  NEXT  I 
580  RT=MM 

590  GOSUB  12000  :  PRINT 
600  FOR  J=1  TO  NF 

610  PRINT  "READING  IN  FILE  — >  "jFl*<J> 

620  OPEN  "I", 1,F1*'J>  :  INPUT  #1,N 

630  FOR  1=1  TO  N  :  INPUT  #1,X(I,J>  :  NEXT  I 

640  CLOSE  1 

650  NEXT  J 

660  N=MM  :  RETURN 

lOOO  REM  DIRECTORY  CHECK  *  * 

1010  IF  F0=0  THEN  RETURN 
1020  P2=0 

1030  FOR  Pl=l  TO  FO 

1040  OPEN  " I " , 1 , F* (P1-P2) 

1050  CLOSE  1 

1060  NEXT  PI 

1070  IF  P2=0  THEN  RETURN  ELSE  GOSUB  3400  :  RETURN 
llOO  REM  *  *  *  CORRECT  DIRECTORY  *  * 

1110  IF  PI— P2=F0  THEN  1130 

1120  FOR  I =P 1 -P2  TO  FO— 1  :  F* ( I ) =F* (1+1 >  :  NEXT  I 
1130  FO=FO— 1  :  P2=P2+1  :  RETURN 
1200  REM  *  *  *  READ  DIRECTORY  *  * 

1210  CLS 

1220  F*="EDITDATA/DIR: "+DN* 

1230  OPEN  " I ” , 1 , F * 

1240  INPUT  #1 , FO 
1250  IF  F0=0  THEN  1290 
1260  FOR  1=1  TO  FO 
1270  INPUT  #1,F*(I) 

1280  NEXT  I 

1290  CLOSE  1  :  GOSUB  1000 
1300  CLOSE  1 
1310  RETURN 

1400  REM  t  *  *  DISPLAY  DIRECTORY  *  * 
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1410  CLS 

1420  IF  F0=0  THEN  1510 

1430  FOR  1*1  TO  FO  STEP  4 

1440  FOR  J=1  TO  4 

1450  IF  I+J— 1 >F0  THEN  1510 

1460  IF  J*1  AND  LEN <F* < I+J— 1 ) ) < 12  AND  I>1  THEN 

PRINT 

1470  PRINT  TAB(16*(J-1)  );”  < " ; STR< < I+J-l >  ;  " > 

" ; F< ( I+J-l ) ; 

1480  NEXT  J 

1490  NEXT  I 
1500  PRINT 
1510  RETURN 

2300  REM  *  *  *  CREATE  DIRECTORY  *  * 

2305  CLOSE  1 

2320  F*="EDITDATA/DIR:  "+DN* 

2340  OPEN  "0",1,F* 

2345  PRINT  #1,0 
2360  CLOSE  1 
2370  RETURN 

3400  REM  *  *  *  ENTER  NEW  FILES  *  * 

3410  F*="EDITDATA/DIR: "+DN* 

3420  OPEN  “0",1,F* 

3430  PRINT  #l,FO 
3440  FOR  1*1  TO  FO 
3450  PRINT  #1,F<<I) 

3460  NEXT  I 
3470  CLOSE  1 
3480  RETURN 

3500  REM  *  *  *  FILE  NAME  CHECK  *  * 

3510  EX=0 

3520  FOR  1=1  TO  FO 

3530  IF  F*=F*<I)  THEN  EX=1  :  RETURN 

3540  NEXT  I 
3550  RETURN 

3600  REM  *  *  *  DIRECTORY  DISPLAY  *  * 

3610  CLS 

3620  INPUT  "DISK  DRIVE  NUMBER  <0, 1,2,3)  TO  FIND  DIRECTORY 
OF " ; DN< 

3630  IF  VAL (DNS) <0  OR  VAL(DN*)>3  THEN  3610 

3640  GOSUB  1200  :  GOSUB  1400  :  PRINT  ” "  :  INPUT  "PRESS 

ENTER  TO  CONTINUE"; J1 

3650  RETURN 

4000  REM  *  *  *  PRINT  DATA  FILES  *  * 

4005  CLS  :  INPUT  "Input  Printer  Width";PW 
4010  PF=0  s  ML=- 1  :  DL*=”  :  TL=-1 

4020  FOR  1=1  TO  N 
4030  FOR  J=1  TO  NF 

4040  IF  LEN<STR<<X (I, J) ) )  >  ML  THEN 

ML=LEN ( STR*  <  X  < I , J ) ) ) 

4050  NEXT  J 
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4060  NEXT  I 

4061  ML=ML+1 

4070  CLSs PRINT: PRINT  "Enter  Column  Headers.  Separating 
lines  with" 

4080  PRINTMthe  delimiter  to  produce  a  two  column 

header " 

4100  FDR  1*1  TO  NF 

4110  PRINT:PRINT  "Enter  the  Column  Header  -for  Column";  I 
4120  INPUT  MS<I) 

4130  NEXT  I 
4140  MS=-1 
4150  FOR  1=1  TO  NF 
4160  P1=0 

4170  P2=INSTR(P1+1,MS(I> ,  DLS) 

4180  L=P2— PI— 1 

4190  IF  P2=0  THEN  L=LEN (MS ( I > > -PI 

4200  IF  L>MS  THEN  MS=L 

4210  IF  P2=0  THEN  4240 

4220  P1=P2 

4230  GOTO  4170 

4240  NEXT  I 

4250  IF  ML>MS  THEN  MX=ML  ELSE  MX=MS 

4260  FOR  1=1  TO  NF 

4270  P1=0  :  NL=0  :  TS="‘* 

4290  P2=INSTR(P1+1,MS(I) ,  DLS) 

4300  L=P2— PI — 1 

4310  IF  P2=0  THEN  L=LEN (MS < I > > -PI 

4320  SF=INT ( (MX— L) /2)  :  SB=MX-L-SF 

4330  TS=TS+STRINGS <SF, " 

")+MIDS(MS(I> , Pl+1 , L) +STRINGS (SB, "  ") 

4335  NL=NL+1  :  IF  NL>TL  THEN  TL=NL 

4340  IF  P2=0  THEN  MS ( I ) =TS  :  GOTO  4370 

4350  P1=P2 

4360  GOTO  4290 

4370  NEXT  I 

4380  PRINT:PRINT  "Would  you  like  to  number  the  observations 
(Y/N)?  "; 

4385  AS= INKEYS: IF  AS=""  THEN  4385  ELSE  PRINT  AS 
4390  IF  AS="Y"  THEN  PI =5  ELSE  PI=0 
4400  MS=INT ( (PW— MXtNF— PI ) /NF) 

4410  PRINT: PRINT"The  maximum  #  o-f  spaces  that  can  be  used 
between  columns  is";MS;"." 

4420  PRINT:  INPUT"Enter  the  number  o-f  spaces  between 
col  unins"  ;SP 

4445  IF  SP<0  THEN  INPUT  "TO  RENAME  COLUMNS  TYPE  Y  ELSE 

N " ; TS  :  IF  TS="Y"  THEN  4010  ELSE  RETURN 

4450  IF  SP<0  OR  SP>MS  THEN  RETURN 

4455  PRINT: INPUT  "Enter  Report  Title";TS 

4472  P1=(NF-1) * (MX+SP) +PI+MX-1 

4474  P1=INT (Pl/2)  :  P2=PW-P1  :  IF  LEN(TS)>PW  THEN 

T*=LEFTS(TS,PW> 
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4476  IF  LEN(T*)>P2  THEN  P1=INT( (PW-LEN(T*) > /2> 

4480  LPRINT  TAB(P1);T* 

4490  LPRINT  TAB(Pl);  :  FOR  1=1  TO  LEN(T$)  :  LPRINT 
jNEXT  I  s  LPRINT 
4500  LPRINT  ““  :  LPRINT 
4505  FOR  J=1  TO  TL 
4510  FOR  1=1  TO  NF 

4520  LPRINT 

TAB ( ( I— 1 ) * (MX+SP) +PI ) ;MID*(MS(I> , (J-l) X+1,MX) ; 

4530  NEXT  I 

4540  LPRINT  H” 

4550  NEXT  J 

4560  A*= " "  :  FOR  1  =  1  TO  MX  :  A*=A*+,,-n  :  NEXT  I 
4570  FOR  1=1  TO  NF 

4580  LPRINT  TAB ( ( 1-1 ) * (MX+SP) +PI ) ; A*; 

4590  NEXT  I 

4600  LPRINT  i  K=PI+MX-1  — INT  < (MX— ML) /2> 

4610  FOR  1=1  TO  N 

4620  IF  PI  =5  THEN  LPRINT  STR* ( I ) ; " : " ; 

4630  FOR  J=1  TO  NF 

4640  LPRINT 

TAB ( (J-l) * (MX+SP) -LEN(STR*(X (I, J) ) >+K) ; X (I, J) ; 

4650  NEXT  J 

4660  LPRINT  •••• 

4670  NEXT  I 
4680  RETURN 

5000  REM  *  *  *  PURGE  DATA  FILES  *  * 

5010  CLS 

5020  PRINT  “Disk  Drive  Containing  Files  to  Be  Purged 
(0-3)?” 

5030  GOSUB  240 

5040  SL=50 10  :  GOSUB  1200  :  IF  SL=0  THEN  GOTO  5010 

5045  K=0 

5046  PRINT  TAB (20);"  *«***«  PURGE  DATA  FILES  t**t" 

5050  FOR  1=1  TO  FO 

5060  A*= " " : PR I NT  F* < I-K) ; "  KILL  IT?  (Y/N/Q)  " ; 

5061  ZI=INKEY*:IF  ZI=""  THEN  5061 

5062  IF  Z I = " Q "  THEN  PRINT  Z I : GOTO  5095 

5064  IF  Z I  = " Y "  THEN  PRINT  Z I :  GOSUB  5100  ELSE  PRINT  **N” 
5090  NEXT  I 
5095  GOSUB  3400 

5099  RETURN 

5100  KILL  F*(I-K) 

5110  IF  I— K=FO  THEN  5130 

5120  FOR  J=I— K  TO  FO-1  :  F* < J ) =F* ( J+l )  :  NEXT  J: 
5130  FO=FO— 1  :  K=K+1 
5140  RETURN 

6000  REM  ***  EXAMINE  DATA  ** 

6005  INPUT  CR 

6010  IF  CR  <  1  THEN  CR=1 

6015  IF  CR  >  N  THEN  CR=N 
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PI =0:  CLS 


6020  E*=CHR*<31)  :  N*=,,M 
6025  C1=1:R1=1:B1=1 
6030  CR~R1:CC=C1 
6035  GOSUB  9500 
6040  GOSUB  6300  :  GOSUB  6350 
6045  P1=0 

6050  A*= INKEY*  :  IF  A*=,,“  THEN  6050 

6055  IF  A*=CHR*(91>  AND  P1=0  THEN  N*=,,,,r  GOSUB  6500  :  GOTO 
6050 

6060  IF  A*=CHR*<27)  AND  P1=0  THEN  N*= CR= 1 : GOSUB 
6520: GOTO  6050 

6065  IF  A*=CHR* (10)  AND  P1=0  THEN  N*=“":  GOSUB  7000  :  GOTO 
6050 

6070  IF  A*=CHR*  (26)  AND  P1=0  THEN  N*="  **:  CR=RT:  GOSUB 
7020: GOTO  6050 

6075  IF  A*=CHR*<8>  AND  P1=0  THEN  N*=”":  GOSUB  7500  :  GOTO 
6050 

6080  IF  A*=CHR* (9)  AND  P1=0  THEN  N*=" “ :  GOSUB  8000  :  GOTO 
6050 

6085  IF  A*=CHR*<24)  AND  P1=0  THEN  N*=” “: CC=1 : GOSUB 
7530: GOTO  6050 

6090  IF  A*=CHR*<25>  AND  P1=0  THEN  N*=“  '* :  CC=CT :  GOSUB 
8030: GOTO  6050 

6095  IF  ( A*  >  '*  /  “  AND  A*<  ”  :  ** )  OR  A*=,,.“  OR  A*="-M  OR  A**'^" 
OR  A*= “E"  THEN  Pl=l  :  GOSUB  9000  :  GOTO  6050 
6100  IF  A*=CHR*<13>  THEN  P1=0:  GOSUB  8500  :  GOTO  6050 
6105  IF  A*=CHR*<8)  AND  Pl=l  THEN  GOSUB  9020  :  GOTO  6050 
6110  IF  A*="D"  AND  P1=0  THEN  N*=" GOSUB  6320:  GOSUB 
lOOOO; GOTO  6050 

6115  IF  A*=” I “  AND  P1=0  THEN  N*= GOSUB  6320: GOSUB 
10500: GOTO  6050 

6120  IF  A*=CHR*(31)  AND  P1=0  THEN  RUN 
6125  IF  A*=''3''  AND  P1=0  THEN  N=RT  :  RETURN 
6130  GOTO  6050 

6300  REM  *  *  *  CLEAR  ENTRY  AREA  *  * 

6310  CC=1 
6320  PRINT  5)2,  " 

6330  PRINT  5)  2,CHR*(95) 

6340  RETURN 

6350  REM  *  *  *  BRACKET  ENTRY  t  * 

6360  K=324+ 1 28* (CR-R1 ) +15* (CC-C1 ) 

6370  TY=15360+K 

6380  POKE  TY, 183  :  POKE  TY+14,187 
6390  RETURN 

6400  POKE  TY, 128  :  POKE  TY+14,128 
6410  RETURN 

6500  REM  *  *  *  PROCESS  UP  ARROW  *  * 

6510  CR=CR— 1 : IF  CR< 1  THEN  CR=1:G0T0  6540 
6520  IF  CR<R1  THEN  R1=CR: RL=R 1+4: GOSUB  9500 
6530  GOSUB  6400: GOSUB  6350 
6540  RETURN 
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7000  REM  *  *  *  PROCESS  DOWN  ARROW  *  * 

7010  CR=CR+1 

7020  IF  CR >RL  THEN  RL=CR: Rl=RL-4; GOSUB  6400: GOSUB  9500 
7030  IF  CR  >  RT  THEN  RT=RT+1 
7040  GOSUB  6400: GOSUB  6350 
7050  RETURN 

7500  REM  *  *  *  PROCESS  LEFT  ARROW  *  * 

7510  CC=CC-1 

7520  IF  CC<1  THEN  CC=1:G0T0  7550 

7530  IF  CC<C1  THEN  C1=CC: CL=C 1+3: GOSUB  9500 

7540  GOSUB  6400: GOSUB  6350 

7550  RETURN 

8000  REM  *  *  *  PROCESS  RIGHT  ARROW  *  * 

8010  CC=CC+1 

8020  IF  CC  >  CT  THEN  GOSUB  9200 

8030  IF  CC>CL  THEN  CL=CC: Cl=CL-3: GOSUB  9500 

8040  GOSUB  6400: GOSUB  6350 

8050  RETURN 

8500  REM  *  *  *  CHANGE  VALUE  IN  FIELD  *  * 

8510  GOSUB  6320 

8520  IF  N*=*"*‘  THEN  GOSUB  8000  :  RETURN 
8530  IF  CR >RT  THEN  RT=CR 
8540  X <CR,CC)=VAL<N*) 

8550  N*= " " 

8560  K=325+12B* (CR-R1 ) +15* (CC-C1 ) 

8570  PRINT  5>K+1,X(CR,CC>; 

8580  RETURN 

8590  K-K-2: PRINT  5)K,  CHR*  ( 128)  ; 

8600  PRINT  USING  ZZ;X(CR,CC>; 

3610  RETURN 

9000  REM  *  *  *  PROCESS  NUMBER  ENTRY  *  * 

9010  N*=N*+A*  :  PRINT  5)2,"  ”  :  PRINT  5)  2,N*;:PRINT 

CHR* (95)  :  RETURN 

9020  IF  LEN(N*)>0  THEN  N*=LEFT*  (NS,  LEN  (N*) -1  >  :  PRINT  5)  2," 

PRINT  5)  2,N* 

9030  RETURN 

9200  REM  *  *  *  LAST  CDLUMN  PROCESSING  *  * 

9210  CC=1  :  Cl  =  CC  :  CL=Cl+3 
9220  CR=CR+1 

9230  IF  CR>RL  THEN  RL=CR  :  Rl=RL-4 
9240  GOSUB  6400 

9250  IF  CT  >  4  OR  CR  =  RL  THEN  GOSUB  9500 
9260  IF  CRJRT  THEN  RT=RT+1 
9270  RETURN 

9500  REM  *  *  *  SCREEN  DISPLAY  *  * 

9510  PRINT  5)  1 96 ,  CHR*  ( 1 43 )  :  IF  CT<Cl+4  THEN  CL=CT  ELSE 

CL=C 1 +3 

9520  RL=Rl+4 

9530  FOR  I=C1-1  TO  CL-1 

9540  PRINT  3197+15* ( I +1 -Cl ) , STRING*  <6, CHR* ( 143) ) ; : PRINT 
CHR*<  128)  ; '.PRINT  CHR*  (65+ 1 )  ;  :  PRINT  CHR*  (  128)  ;:  PRINT 


STRING*  (6, CHR*  ( 143) ) 

9550  NEXT  I 

9560  PRINT  3195, CHR* (170) : PRINT  3259, CHR* ( 1 70) 

9570  FOR  I=R1— 1  TO  RL-1 

9580  IF  I  >  8  THEN  K1*0  ELSE  Kl=l 

9590  PRINT  3320+128* ( I+l-Rl > +K1, STR* < 1  +  1 ); : PRINT  CHR* (170) 
9600  FOR  J=C1-1  TO  CL-1 

9610  PRINT  3326+15* ( J+l— Cl ) +128* ( 1+1— R1 ),X(I+1,J+1) 

9620  NEXT  J 

9630  PRINT  3387+128* ( I+l-Rl ), CHR* (170) 

9640  NEXT  I 
9650  RETURN 

lOOOO  REM  *  *  *  DELETE  A  LINE  *  * 

10010  FOR  I=CR  TO  RT-1:F0R  J=1  TO  CT: X ( I , J) =X ( 1+1 , J)  :  NEXT 
J  :  NEXT  I 

10020  FOR  J=1  TO  CT:  X(RT,J)=0  :  NEXT  J 

10030  RT=RT-1  :  IF  R1  >  1  AND  RL  >  RT  THEN  R1=R1-1  : 

CR=CR— 1 

10040  GOSUB  9500  :  GOSUB  6350 
10050  RETURN 

10500  REM  *  *  *  INSERT  A  LINE  *  * 

10510  FOR  I=RT  TO  CT  STEP  -1  :  FOR  J=1  TO  CT  : 

X(I+1,J)=X(I,J)  :  NEXT  J  :  NEXT  I 

10520  FOR  J=1  TO  CT  :  X(CR,J)=0  :  NEXT  J 
10530  RT=RT+1  :  GOSUB  9500  :  GOSUB  6350 
10540  RETURN 

11000  REM  *  *  *  VERIFY  FILENAME  *  * 

11010  FW=0 

11020  P1=INSTR(1,S*, "/"> 

11030  IF  PI  <>  0  THEN  IF  LEN (S*> -PI >3  THEN  FW=1  :  RETURN 

11040  IF  PI  >8  THEN  FW=1  :  RETURN 

11050  IF  P1=0  THEN  IF  LEN(S*)>8  THEN  FW=1  :  RETURN 
11060  IF  PI  <>  O  THEN  IF  ASC (MID* (S*, Pl+1 , 1 ) ) <65  OR 
ASC (MID* (S*, Pl+1 , 1 ) ) >90  THEN  FW=1  :  RETURN 
11070  IF  ASC ( LEFT * ( S* ,  1 )  ) < 65  OR  ASC (LEFT* (S*, 1 ) ) >90  THEN 
FW=*1  :  RETURN 

11080  IF  P1=0  THEN  P2=LEN(S*)  ELSE  P2=P1-1 
11090  P3=l 

11100  FOR  I =P3  TO  P 2 
11110  A*=MID* (S*, 1,1) 

11120  IF  (ASC (A*) >47  AND  ASCCA*)<58)  OR  (ASC (A*) >64 

AND  ASC (A*) <91  )  THEN  11140 
11130  FW= 1  :  RETURN 

11140  NEXT  I 

11150  IF  P3=l  AND  PI  <>  O  THEN  P3=P1+1  :  P2=LEN(S*)  :  GOTO 
11100 

11160  RETURN 

12000  REM  *  *  *  DIMENSION  ARRAYS  *  * 

12010  NV=MEM— 500  :  NV=INT (NV/ (4* (NF+l ) ) ) -1 

12020  IF  MM  >  NV  THEN  PRINT  "INSUFFICIENT  SPACE  TO  READ  IN. 
INPUT  "PRESS  ENTER  TO  CONTINUE" ; A* 


12030  DIM  X ( N V , NF ) 

12040  MN=NV 
12050  RETURN 

13000  REM  *  *  *  RENAME  FILES  ROUTINE  *  * 

13010  CLS 

13020  PRINT  "Disk  Drive  To  Contain  Files  (0-3)?" 

13030  60SUB  240 

13040  SL=8240: GOSUB  1200  s  BOSUB  1400  :  PRINT  "" 

13050  FOR  J=1  TO  NF 

13060  PRINT  “Variable  4" ; J; "Currently  Uses  Filename 

— >";F1*(J) 

13070  INPUT  "New  Name";S*  :  GOSUB  11000 

13080  IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO 

13060 

13090  F*=S*  s  BOSUB  3500 

13100  IF  EX= 1  THEN  PRINT : PRINT  “A  File  Currently  Uses 

the  New  Fi lename. PRINT  "Da  You  Want  to  Use  it  Anyway 
( Y/N) ?"  ELSE  GOTO  13120 

13110  A*=INKEY*:IF  A«=“"  THEN  13110  ELSE  PRINT  A*:  IF  A*< > 

"Y”  THEN  13060 

13120  F1*(J)=F* 

13130  NEXT  J 
13140  RETURN 

14000  REM  *  *  *  NAME  FILES  ROUTINE  *  * 

14010  CLS 

14020  PRINT  "Disk  Drive  To  Contain  Files  (0-3) ?": GOSUB  240 
14030  SL=8330  :  GOSUB  1200  :  BOSUB  1400  :  PRINT  “" 

14040  PRINT 

14050  FOR  J=1  TO  NF 

14060  PRINT  "FILENAME  FOR  VARIABLE  4";J  :  INPUT  S*  :  GOSUB 
llOOO 

14070  IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO  14060 
14080  F*=S*  :  GOSUB  3500 

14090  IF  EX=1  THEN  PRINT  "FILENAME  IS  ALREADY  USED.  TYPE 
Y  TO  USE  NAME  ANYWAY  ELSE  N";:  INPUT  A*  :  IF  A*  <>  "Y"  THEN 
1 4060 

14100  F1*(J)=F* 

14110  NEXT  J 
14120  RETURN 

15000  REM  *  *  *  EXITING  ROUTINE  FOR  OPTIONS  1  S<  2  *  * 

15010  CLS 

15020  PRINT  "Do  You  Wish  To  Rename  Files  (Y/N)?" 

15030  ZI=INKEY$: IF  ZI=” "  GOTO  15030 
15040  IF  Z I = " Y "  THEN  GOSUB  13000 

15050  SL=8550  :  GOSUB  1200  :  GOSUB  1400  :  PRINT  ““ 

15060  PRINT: PRINT 
15070  FOR  K= 1  TO  NF 

15080  PRINT  "WRITING  TO  DISK  — >  ";F1*(K> 

15090  OPEN  "O", 1,F1*(K>+": “+DN* 

15100  PRINT  41, N  :  FOR  1=1  TO  N  :  PRINT  41,X(I,K);  : 

NEXT  I 
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15110  CLOSE  1  :  F*=F1*(K)  :  GOSUB  3500  :  IF  EX=1  THEN 

15130 

15120  FO=FO+1  :  F* <F0> =F1* <K>  :  GOSUB  3400 

15130  NEXT  K 
15140  RETURN 

19000  'ERROR  HANDLING  ROUTINE 

19015  IF  ERR=10  AND  ERL=8520  THEN  N*=,,“  :  RESUME  8530 

19020  IF  ERR  <>  106  THEN  19100 

19021  IF  ERL=5100  THEN  RESUME  5110 

19022  IF  ERL=1040  THEN  GOSUB  1100  :  RESUME  1050 
19025  IF  ERL  <>  1230  THEN  19100 

19030  IF  SL=1715  OR  SL=5010  OR  SL=5410  THEN  PRINT  "NO 
FILES  ON  THIS  DRIVE.  PRESS  ENTER  TO  CONTINUE":  INPUT  A*  : 
SL=0  :  RESUME  1300 

19045  IF  SLIsAS240  OR  SL=8330  OR  SL=8550  THEN  GOSUB  2300 
RESUME  1210 

19100  PRINT  " ERR= " ; ERR ; " ERL= " ; ERL  :  ON  ERROR  GOTD  O  : 
RESUME 


54 


PROGRAM  3 


CTABLE/BAS 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFSTR  Z 

40  DIM  F*(95) ,M*(11) ,N*(11) 

50  LB*=CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL*=CHR* ( 252 ) 
60  Z*="######" :  ZZ=“##.  ##C  I  C I " 

70  CLS 

SO  PRINT  TAB (20) ; "DATA  FILE  MANAGEMENT  PROGRAM" 

90  PRINT: PRINT: PRINT  "MENU:" 

100  PRINT  "(1)  Create  New  Data  Files" 

110  PRINT  "(2)  Edit  Existing  Files" 

120  PRINT  "(3)  Purge  Existing  Files" 

130  PRINT  "(4)  Print  Data  Files" 

140  PRINT  "(5)  Return  to  Master  Menu" 

150  KP=896: KS=5: GOSUB  190 

160  ON  10  GOSUB  1000,1100,1200,1300,180 

170  GOTO  70 

180  RUN  "MENU/BAS" 

190  REM  *  *  KEYBOARD  INPUT  SUBROUTINE  * 

200  PRINT  5>KP-*-5,  "OPTION  DESIRED  " ; LB* ; MB*; RB*; 

210  ZI= INKEY*: IF  ZI  =  ""  THEN  210  ELSE  IO=VAL(ZI) 

220  PRINT  3KP+21, ZI; : FOR  1=1  TO  100: NEXT 

230  IF  I0<1  OR  IO>KS  THEN  PRINT  S)KP,  BL*;  :  PRINT  3KP, "** 

NOT  A  VALID  OPTION  ***";:  FOR  1=1  TO  1000:  NEXT:  FRINT 
3KP,  BL*; : GOTO  200 
240  RETURN 

250  DN*= INKEY*: IF  DN*=""  THEN  250  ELSE  I=VAL(DN*> 

260  IF  ICO  OR  I >3  THEN  250  ELSE  RETURN 
lOOO  REM  *  *  *  NEW  FILE  ROUTINE  *  * 

1010  CLS  :  INPUT  "How  many  rows" ; RT 
1020  INPUT  "How  many  columns";CT 
1030  BT  = 1 

1040  GOSUB  14000  :  GOSUB  12000  :  N=i  :  CR=1  :  GOSUB  6040 

1050  CLS: GOSUB  15050: RUN 

1100  GOSUB  1500  :  CR=1  :  GOSUB  6040 

1110  GOSUB  15000  :  RUN 

1200  GOSUB  5000  :  RUN 

1300  GOSUB  1500  :  GOSUB  4000  :  RUN 

1500  REM  *  *  *  READ  DATA  ARRAY  *  * 

1510  DL*=" / " 

1520  CLS: PRINT  "Disk  Drive  Containing  Files  (0-3)?" 

1530  GOSUB  250 

1540  SL=1715: GOSUB  2400: IF  SL=0  THEN  GOTO  1520 
1550  GOSUB  2800  :  PRINT  "" 

1560  PRINT: IF  10  =  2  PRINT  "Enter  the  File  To  Be  Edited"; 
ELSE  PRINT  "Enter  the  File  To  Be  Printed"; 
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1570  INPUT  NV 

1580  IF  NV>F0  THEN  PRINT : PRINT “  *  t  *  INVALID  FILE  #  * 

* “ s  GOTO  1560 

1590  F1S=F*<NV> 

1600  OPEN  "I",1,F14 
1610  INPUT  #1 , RT 
1620  INPUT  #1,CT 
1630  INPUT  #1,BT 
1640  CLOSE  1 

1650  GOSUB  12000  :  PRINT 
1660  OPEN  “I", 1 ,  F 1  * 

1670  INPUT  #1,RT: INPUT  #1,CT: INPUT  #1,BT 
1680  FOR  1=1  TO  RT 
1690  FOR  J=1  TO  CT 
1700  INPUT  #1,X(I,J> 

1710  NEXT  J 
1720  NEXT  I 
1730  CLOSE  1 
1740  RETURN 

2000  REM  *  *  *  CHECK  FOR  FILE  DIRECTORY  *  * 

2010  IF  F0=0  THEN  RETURN 
2020  P2=0 

2030  FOR  Pl=l  TO  FO 

2040  OPEN  "I", 1,F*(P1-P2) 

2050  CLOSE  1 

2060  NEXT  PI 

2070  IF  P2=0  THEN  RETURN  ELSE  GOSUB  3400  :  RETURN 
2080  ’SUB  TO  CORRECT  DIRECTORY 
2090  IF  PI— P2=F0  THEN  2110 

2100  FOR  I=P1— P2  TO  FO-1  :  F*  ( I )  =F*  ( I  1 )  :  NEXT  I 

2110  F0=F0— 1  :  P2=P2+1  :  RETURN 

2400  REM  *  *  *  READ  DIRECTORY  OF  FILES  t  * 

2410  CLS 

2420  F*="CTABLE/DIR: "+DN* 

2430  OPEN  " I " , 1 , F* 

2440  INPUT  #1 , FO 
2450  IF  F0=0  THEN  2490 
2460  FOR  1=1  TO  FO 
2470  INPUT  #1,F*<I) 

2480  NEXT  I 

2490  CLOSE  1  :  GOSUB  2000 
2500  CLOSE  1 
2510  RETURN 

2800  REM  *  *  *  DISPLAY  DIRECTORY  *  * 

2810  CLS 

2820  IF  F0=0  THEN  2910 

2830  FOR  1=1  TO  FO  STEP  4 

2840  FOR  J=1  TO  4 

2850  IF  I+J-l >F0  THEN  2910 

2860  IF  J=1  AND  LEN <F* ( I+J-l ) ) < 12  AND  I>1  THEN 

PRINT 
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2870  PRINT  TAB(16*(J-1) ) ; " ; STR* (I+J-l > ; ") 

F*  < I+J— 1 ) ; 

2880  NEXT  J 

2890  NEXT  I 
2900  PRINT 
2910  RETURN 

3200  REM  *  *  *  CREATE  “CTABLE/DIR"  *  * 

3210  DN*="0“ 

3220  CLOSE  1 

3230  F*= "CTABLE/DIR: "+DN* 

3240  OPEN  "0“ ,  1 , F* 

3250  PRINT  #1,0 
3260  CLOSE  1 
3270  RETURN 

3400  REM  *  *  *  WRITE  NEW  FILES  IN  DIRECTORY  *  * 

3410  F*=*'CTABLE/DIR: "+DN* 

3420  OPEN  "0",1,F* 

3430  PRINT  #1 , FO 

3440  FOR  1=1  TO  FO 

3450  PRINT  #1,F*(I)  ;  ",  **; 

3460  NEXT  I 
3470  CLOSE  1 
3480  RETURN 

3600  REM  *  *  *  CHECK  IF  FILENAME  WAS  ALREADY  USED  *  * 
3610  EX=0 

3620  FOR  1=1  TO  FO 

3630  IF  F*=F*<I)  THEN  EX=1  :  RETURN 

3640  NEXT  I 
3650  RETURN 

4000  REM  *  *  *  PRINT  ROUTINE  *  * 

4010  CLS  :  INPUT  "Input  Printer  Width";PW 
4020  PF=0  :  ML=-1  :  DL*=" ! "  :  TL=-1 
4030  FOR  1=1  TO  N 
4040  FOR  J=1  TO  CT 

4050  IF  LEN<STR*(X<I, J) ) )  >  ML  THEN 

ML=LEN(STR*<X (I, J) ) ) 

4060  NEXT  J 

4070  NEXT  I 
4080  ML=ML+1 

4090  CLS: PRINT : PRINT  "Enter  Column  Headers.  Separating 
lines  with" 

4100  PRINT"the  delimiter  ’ ? ’  to  produce  a  two  column 
header " 

4110  FOR  1=1  TO  CT 

4120  PRINTzPRINT  "Enter  the  Column  Header  -for  Column";  I 
4130  INPUT  M*(I) 

4140  NEXT  I 
4150  MS=-1 
4160  FOR  1=1  TO  CT 
4170  P1=0 

4180  P2=  I NSTR  ( P 1  -t- 1 ,  M*  <  I )  ,  DL* ) 
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4190  L=P2-P1-1 

4200  IF  P2=0  THEN  L=LEN (MS < I > ) -PI 

4210  IF  L>MS  THEN  MS=L 

4220  IF  P2=0  THEN  4250 

4230  P1=P2 

4240  GOTO  4180 

4250  NEXT  I 

4260  IF  ML>MS  THEN  MX=ML  ELSE  MX=MS 

4270  FOR  I«1  TO  CT 

4280  P1=0  :  NL=0  :  TS="" 

4290  P2=INSTR(P1+1,MS(I> ,DLS) 

4300  L=P2-P1— 1 

4310  IF  P2=0  THEN  L=LEN (MS < I ) ) -PI 

4320  SF=INT< (MX-L)/2)  :  SB=MX-L-SF 

4330  TS=TS+STR I NGS ( SF  ,  " 

")+MIDS(MS(I)  ,P1+1,L)+STRINGS(SB,  H  ”> 

4340  NL=NL+1  :  IF  NL>TL  THEN  TL=NL 

4350  IF  P2=0  THEN  MS < I ) =TS  :  GOTO  4380 

4360  P1=P2 

4370  GOTO  4290 

4380  NEXT  I 

4390  PRINT : PRINT  "Would  you  like  to  number  the  rows  <Y/N)? 
» . 

4400  AS= INKEYS: IF  AS*”"  THEN  4400  ELSE  PRINT  AS 
4410  IF  AS="Y"  THEN  PI=5  ELSE  PI=0 
4420  MS* I NT (  (PW-MX  *  CT-PD/CT) 

4430  PRINT: PRINT'*The  maximum  #  of  spaces  that  can  be  used 
between  columns  is";MS;”." 

4440  PRINT: INPUT"Enter  the  number  of  spaces  between 
columns” ; SP 

4450  IF  SP<0  THEN  INPUT  ”T0  RENAME  COLUMNS  TYPE  Y  ELSE 

N " ; TS  :  IF  T*="Y"  THEN  4020  ELSE  RETURN 

4460  IF  SP<0  OR  SP>MS  THEN  RETURN 

4470  PRINT: INPUT  "Enter  Report  Title”;TS 

4480  P1=<CT-1)*(MX+SP)+PI+MX-1 

4490  P1  =  INT  <Pl/2)  :  P2=PW-P1  :  IF  LEN(TS)>PW  THEN 
T*=LEFT*(T*,PW) 

4500  IF  LEN  <TS) >P2  THEN  P1  =  INT C (PW-LEN (TS) ) /2) 

4510  LPRINT  TAB (PI ) ; TS 

4520  LPRINT  TAB (PI);  :  FOR  1=1  TO  LENCTS)  :  LPRINT 
; : NEXT  I  :  LPRINT 
4530  LPRINT  ""  :  LPRINT  "" 

4540  FOR  J=1  TO  TL 

4550  FOR  1  =  1  TO  CT 

4560  LPRINT 

TAB ( <1-1 ) * (MX+SP) +PI ) ; MID* (MS ( I ) , ( J— 1 ) X+l , MX) ; 

4570  NEXT  I 

4580  LPRINT  ” " 

4590  NEXT  J 

4600  AS*”"  :  FOR  1=1  TO  MX  :  AS=AS+"-"  :  NEXT  I 
4610  FOR  1=1  TO  CT 
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4620  LPRINT  TAB ( ( 1-1 > * (MX+SP) +PI > ; A*; 

4630  NEXT  I 

4640  LPRINT  ""  :  K=PI+MX-1-INT ( (MX— ML)  /2) 

4650  FOR  1=1  TO  RT 

4660  IF  PI=5  THEN  LPRINT  STR* ( I ) ; " : “ ; 

4670  FOR  J=1  TO  CT 

4680  LPRINT 

TAB ( ( J— 1 ) * (MX+SP) — LEN (STR* (X  ( I , J ) ))+K) ;X<I,J) ; 

4690  NEXT  J 

4700  LPRINT  "" 

4710  NEXT  I 
4720  RETURN 

5000  REM  *  *  *  PURGE  DATA  FILE  DIRECTORY  *  * 

5010  CLS 

5020  PRINT  "Disk  Drive  Containing  Files  to  Be  Purged 
(0-3)?" 

5030  GOSUB  250 

5040  SL=5010  :  GOSUB  2400  IF  SL=0  THEN  GOTO  5010 
5050  K=0 

5060  PRINT  TAB (20) ; "  ******  PURGE  DATA  FILES  »***" 

5070  FOR  1=1  TO  FO 

5080  A*="":PRINT  F*(I-K);"  KILL  IT?  (Y/N/Q)  " ; 

5090  ZI= INKEY*: IF  ZI=""  THEN  5090 
5100  IF  ZI="Q"  THEN  PRINT  ZI:GOTO  5130 

5110  IF  ZI="Y"  THEN  PRINT  ZI:GOSUB  5150  ELSE  PRINT  "N" 

5120  NEXT  I 
5130  GOSUB  3400 
5140  RETURN 
5150  KILL  F* ( I — K ) 

5160  IF  I— K=F0  THEN  5180 

5170  FOR  J=I-K  TO  FO-l  :  F* ( J ) =F* ( J+l )  :  NEXT  J : 

5180  F0=F0— 1  :  K=K+1 
5190  RETURN 

6000  REM  ***  EXAMINE  DATA  ** 

6010  INPUT  CR 

6020  IF  CR  <  1  THEN  CR=1 

6030  IF  CR  >  N  THEN  CR=N 

6040  E*=CHR*(31>  :  N*^'"  :P1=0:  CLS 

6050  C1=1:R1=1:B1=1 

6060  CR=R1:CC=C1 

6070  GOSUB  6760 

6080  GOSUB  6240  :  GOSUB  6290 

6090  P1=0 

6100  A*=  INKEY*  :  IF  A*=,,M  THEN  6100 

6110  IF  A*=CHR* (91 )  AND  P1=0  THEN  N*="":  GOSUB  6360  :  GOTO 
6100 

6120  IF  A*=CHR* (27)  AND  P1=0  THEN  N*=" ": CR=1 : GOSUB 
6380: GOTO  6100 

6130  IF  A*=CHR*  (10)  AND  P1=0  THEN  N*=,,,,:  GOSUB  6410  :  GOTO 
6100 

6140  IF  A*=CHR* (26)  AND  P1=0  THEN  N*= CR=RT : GOSUB 
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6440s GOTO  6100 

6150  IF  A*»CHR*<8>  AND  P1=0  THEN  NS*"":  SOSUB  6470  :  GOTO 
6100 

6160  IF  A*«CHR*<9>  AND  P1=0  THEN  N*=M,,s  GOSUB  6530  :  GOTO 
6100 

6170  IF  A*=CHR*<24>  AND  P1=0  THEN  N*= “ " : CC= 1 : GOSUB 
6500: GOTO  6100 

6180  IF  A*=CHR*<25)  AND  P1=0  THEN  N*=M CC=CT: GOSUB 
6560: GOTO  6100 

6190  IF  <A*>"/"  AND  A*<*':M)  OR  A*="EM  THEN  Pl  =  l  :  GOSUB 
6720  :  GOTO  6100 

6200  IF  A*~CHR*<13>  THEN  P1=0:  GOSUB  6590  :  GOTO  6100 
6210  IF  A*~CHR*(8>  AND  Pl=l  THEN  GOSUB  6740  :  GOTO  6100 
6220  IF  A*=M3>“  AND  P1=0  THEN  RETURN 
6225  IF  A*=CHR*<31>  AND  P1=0  THEN  RUN 
6230  GOTO  6100 

6240  REM  *  *  *  CLEAR  ENTRY  AREA  *  * 

6250  CC=1 
6260  PRINT  32, ” 

6270  PRINT  3  2,CHR*(95> 

6280  RETURN 

6290  REM  *  *  *  BRACKET  ENTRY  *  * 

6300  K=320+ 1 28*  < CR-R 1 ) + 1 1  *  <  CC-C 1 >  +5 
6310  TY=15360+K 

6320  POKE  TY, 183  :  POKE  TY+1 1,187 
6330  RETURN 

6340  POKE  TY, 128  :  POKE  TY+1 1,128 
6350  RETURN 

6360  REM  *  *  *  PROCESS  UP  ARROW  *  * 

6370  CR=CR— 1 : IF  CR< 1  THEN  CR=1:G0T0  6400 
6380  IF  CR<R1  THEN  R1=CR: RL=Rl+4: GOSUB  6760 
6390  GOSUB  6340: GOSUB  6290 
6400  RETURN 

6410  REM  *  *  *  PROCESS  DOWN  ARROW  *  * 

6420  CR=CR+ 1 

6430  IF  CR>RT  THEN  CR=RT:GOTO  6460 

6440  IF  CR >RL  THEN  RL=CR: Rl=RL-4: GOSUB  6340: GOSUB  6760 
6450  GOSUB  6340: GOSUB  6290 
6460  RETURN 

6470  REM  *  *  *  PROCESS  LEFT  ARROW  *  * 

6480  CC=CC-1 

6490  IF  CC<1  THEN  CC=1:G0T0  6520 

6500  IF  CC<C1  THEN  C1=CC: CL=C1 +4 : GOSUB  6760 

6510  GOSUB  6340: GOSUB  6290 

6520  RETURN 

6530  REM  *  *  *  PROCESS  RIGHT  ARROW  *  * 

6540  CC=CC+1 

6550  IF  CC  >  CT  THEN  GOSUB  7000 

6560  IF  CC>CL  THEN  CL=CC: Cl =CL-4: GOSUB  6760 

6570  GOSUB  6340: GOSUB  6290 

6580  RETURN 
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6590  REM  t  *  *  CHANGE  VALUE  IN  FIELD  *  * 

6600  GOSUB  6260 

6610  IF  N*=""  THEN  GOSUB  6530  :  RETURN 
6620  X(CR,CC)-VAL(N*) 

6630  N*=M" 

6640  K=«326+ 1 28*  < CR-R 1 )  + 11  *  ( CC-C 1 ) 

6650  IF  X(CR, CC> >999999  THEN  GOTO  6690 
6660  PRINT  3K+1 ,  CHR*  (  128)  ; 

6670  PRINT  USING  Z*;X(CR,CC); 

6680  RETURN 

6690  K=K— 2: PRINT  3K, CHR* ( 128) ; 

6700  PRINT  USING  ZZ;X(CR,CC); 

6710  RETURN 

6720  REM  *  *  *  PROCESS  NUMBER  ENTRY  *  * 

6730  N*=N*+A*  :  PRINT  32,  “  "  :  PRINT  3  2,N*;:PRINT 

CHR* (95)  :  RETURN 

6740  IF  LEN  (N*)  >0  THEN  N*=LEFT*  (N*,  LEN  (N*) -1 )  :  PRINT  3  2, '* 
PRINT  3  2,N* 

6750  RETURN 

6760  PRINT  3196, CHR* (143)  :  IF  CT<Cl+5  THEN  CL=CT  ELSE 
CL=Cl+4 

6770  IF  RT<Rl+4  THEN  RL  =  RT  ELSE  RL=Rl+4 
6780  FOR  I=C1-1  TO  CL-1 

6790  PRINT  3197+11*  ( I  +  l-Cl ),  STRING*  (4,CHR*(  143)  );  .-PRINT 
CHR* ( 1 28 ) ; : PR I NT  CHR* (65+1 );: PRINT  CHR* ( 128) PRINT 
STRING* (4, CHR* (143) ) 

6800  NEXT  I 

6810  PRINT  3195, CHR* ( 170) : PRINT  3259, CHR* ( 170) 

6820  FOR  I=R1— 1  TO  RL-1 

6830  IF  I  >  8  THEN  K1=0  ELSE  Kl=l 

6840  PRINT  3320+128* ( I  +  1  — R1 ) +K1 . STR* ( 1  +  1 ) ; : PRINT  CHR*(170> 
6850  FOR  J=C1-1  TO  CL-1 

6860  PRINT  3327+11* ( J+l-Cl > +128* (I+l-Rl) , CHR* (128) ; :PRINT 
USINGZ*; X(I+1,J+1) 

6870  NEXT  J 

6880  PRINT  3387+128* ( I+l-Rl ) ,CHR*<170) 

6890  NEXT  I 
6900  RETURN 

7000  REM  *  *  *  PROCESS  LAST  COLUMN  ENTRY  *  * 

7010  CC=1  :  C1=CC  :  CL=Cl+4 
7020  CR=CR+ 1  :  IF  CR  >  RT  THEN  CR=RT 
7030  IF  CR  >  RL  THEN  RL=CR  :  Rl=RL-4 
7040  GOSUB  6340 

7050  IF  CT  >  5  OR  CR  =  RL  THEN  GOSUB  6760 
7060  RETURN 

11000  REM  *  *  *  VERIFY  FILENAME  -  FW=1  IF  FAILS  *  * 

11010  FW=0 

11020  P 1 = I NSTR ( 1 , S* , " / " > 

11030  IF  PI  <>  O  THEN  IF  LEN(S*)-P1>3  THEN  FW=1  :  RETURN 
11040  IF  PI  >8  THEN  FW=1  :  RETURN 

11050  IF  P1=0  THEN  IF  LEN(S*)>8  THEN  FW=1  :  RETURN 
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11060  IF  PI  <>  O  THEN  IF  ASC  (MID*  (S*,  PH-1 , 1  >  >  <65  OR 
ASC(MID*(S*,P1+1, 1) > >90  THEN  FW*1  :  RETURN 
11070  IF  ASC (LEFT* (S*, 1 ) ) <65  OR  ASC (LEFT* (S*, 1 ) ) >90  THEN 
FW=1  :  RETURN 

11080  IF  P1=0  THEN  P2=LEN(S*)  ELSE  P2=P1-1 
11090  P3=l 

11100  FOR  I=P3  TO  P2 
11110  A*=M I D* ( S* ,1,1) 

11120  IF  (ASC (A*) >47  AND  ASC(A*)<58)  OR  (ASC (A*) >64 

AND  ASC ( A* ><91  )  THEN  11140 
11130  FW®1  :  RETURN 

11140  NEXT  I 

11150  IF  P3=l  AND  PI  <>  0  THEN  P3=P1+1  :  P2=LEN(S*>  :  GOTO 
11100 

11160  RETURN 

12000  REM  *  *  *  DIMENSION  ARRAYS  *  * 

12010  NV=MEM— 500  :  NV=INT <NV/ (4* (RT  *  CT  *  BT+1>>)-1 
12020  IF  MM  >  NV  THEN  PRINT  "INSUFFICIENT  SPACE  TO  READ  IN. 

INPUT  "PRESS  ENTER  TO  CONTINUE" ; A* 

12030  DIM  X (RT, CT) 

12040  MN=NV 
12050  RETURN 

13000  REM  *  *  *  RENAME  FILES  *  * 

13010  CLS 

13020  PRINT  "Disk  Drive  To  Contain  Files  (0-3)?" 

13030  GOSUB  250 

13040  SL=8240: GOSUB  2400  :  GOSUB  2800  :  PRINT  " " 

13050  INPUT  "New  Naroe";S*  :  GOSUB  HOOO 

13060  IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO 

13050 

13070  F*=S*  :  GOSUB  3600 

13080  IF  EX=1  THEN  PRINT : PRINT  "A  File  Currently  Uses 

the  New  Fi lenarae. ": PRINT  "Do  You  Want  to  Use  it  Anyway 
<Y/N>?"  ELSE  GOTO  13100 

13090  A*= INKEY*: IF  A*=n"  THEN  13090  ELSE  PRINT  At:  IF  A4< > 
"Y"  THEN  13050 
13100  F 1 *=F* 

13110  RETURN 

14000  REM  *  t  t  NAME  FILE  *  * 

14010  CLS 

14020  PRINT  "Disk  Drive  To  Contain  File  (0-3) ?": GOSUB  250 
14030  SL=8330  :  GOSUB  2400  :  GOSUB  2800  :  PRINT  “ " 

14040  PRINT 

14050  PRINT  "FILENAME  FOR  TABLE”;  :  INPUT  S*  :  GOSUB  HOOO 
14060  IF  FW=1  THEN  PRINT  "INVALID  FILENAME”  :  GOTO  14050 
14070  F*=S*  :  GOSUB  3600 

14080  IF  EX=1  THEN  PRINT  "FILENAME  IS  ALREADY  USED.  TYPE 

Y  TO  USE  NAME  ANYWAY  ELSE  N";:  INPUT  A*  :  IF  A*  <>  "Y"  THEN 

14050 

14090  F1*=F* 

14100  RETURN 
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15000  REM  *  *  *  EXITING  ROUTINE  FOR  OPTIONS  1  St  2  *  * 
15010  CLS 

15020  PRINT  “Do  You  Wish  To  Rename  File  <Y/N>?“ 

15030  ZI=INKEY*:IF  ZI=“”  GOTO  15030 
15040  IF  Z I  = " Y “  THEN  GOSUB  13000 

15050  SL=8550  :  GOSUB  2400  :  GOSUB  2BOO  :  PRINT  ,,n 
15060  PRINT: PRINT 

15070  PRINT  “WRITING  TO  DISK  — >  ";F1* 

15080  OPEN  ”0”, l,Fl*+“: “+DN* 

15090  PRINT  #1,RT: PRINT  #1 , CT: PRINT# 1 , BT 
15100  FOR  1=1  TO  RT 
15110  FOR  J=1  TO  CT 
15120  PRINT  #1,X(I,J); 

15130  NEXT  J 
15140  NEXT  I 

15150  CLOSE  1  :  F*=F1*  :  GOSUB  3600  :  IF  EX=1  THEN 

15170 

15160  F0=F0+1  :  F*<F0)=F1*  :  GOSUB  3400 

15170  RETURN 

19000  REM  *  *  *  ERROR  HANDLING  ROUTINE  *  * 

19010  IF  ERR=10  AND  ERL=6590  THEN  N*= “ "  :  RESUME  6630 

19020  IF  ERR  <>  106  THEN  19080 

19030  IF  ERL=5150  THEN  RESUME  5160 

19040  IF  ERL=2040  THEN  GOSUB  2080  :  RESUME  2050 

19050  IF  ERL  <>  2430  THEN  19080 

19060  IF  SL=1715  OR  SL=5010  OR  SL=5410  THEN  PRINT  "NO 
FILES  ON  THIS  DRIVE.  PRESS  ENTER  TO  CONTINUE":  INPUT  A* 
SL=0  :  RESUME  2500 

19070  IF  SL=8240  OR  SL=8330  OR  SL=8550  THEN  GOSUB  3200 
RESUME  2410 

19080  PRINT  " ERR= " ; ERR ; " ERL= " ; ERL  :  ON  ERROR  GOTO  0  : 
RESUME 


63 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFINT  PsDEFSTR  Z: DEFDBL  F,L 

40  DIM  S ( 1500, 1 )  , F* (95) , P (1500) ,  X  ( 1500) , B  <50 ) 

50  LB*=CHR*(123) : RB*=CHR* ( 125) : MB*=CHR* ( 176) : BL*=CHR* (252) 
60  REM  *  *  *  WILCOX SON  SIGNED  RANKS  TEST  *  * 

70  CLS:  PRINT  320,  "  WILCOX  SON  SIGNED  RANKS  TEST** 

75  PR I NT : PR I NT : PR I NT  "Disk  Drive  Containing  Data  Files 
(0-3) ?" 

80  DN*= INKEY*: IF  DN*="“  THEN  SO  ELSE  I=VAL(DN*) 

85  IF  ICO  OR  I >3  THEN  75 

90  SL= 1 00  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT 

100  PRINT: INPUT  "Enter  the  Two  Data  Files  (e.g.  3,5>";F1,F2 
110  IF  FI >F0  OR  F2>F0  OR  F1<1  OR  F2< 1  THEN  lOO 
120  OPEN  " I", 1,F*(F1)+": "+DN* 

130  INPUT  #1 ,  N 

140  FOR  1=1  TO  N  :  INPUT  #1,S(I,0)  :  NEXT  I 
150  CLOSE  1 

160  OPEN  " I " , 1 , F* (F2) +  "  : “+DN* 

170  INPUT  #1.M 

180  FOR  1=1  TO  M  :  INPUT  #1,S(I,1)  :  NEXT  I 

190  CLOSE  1 

200  IF  N=M  GOTO  250 

210  PRINT : PRINT  "The  Data  Sets  are  NOT  a-f  Equal  Size.  Do 
You  Want  to  Continue  using  the  Smaller  #  o-f  data  points. 
Enter  other  Data  Files,  or  Return  to  the  Master  Menu 
(C/E/R)?"; 

220  ZI=INKEY*:IF  ZI=""  THEN  220 

230  IF  ZI  =  *'C"  THEN  PRINT  ZI:GOTO  240:ELSE  IF  ZI  =  "E"  THEN 

PRINT  "E" : RUN  ELSE  IF  ZI"R"  THEN  PRINT  "R" : RUN"MENU/BAS" 

ELSE  GOTO  220 

240  IF  N  >  M  THEN  N=M 

250  J=i 

260  FOR  1=1  TO  N 

270  S(J, 1)=S(I,0)-S(I, 1) 

280  IF  S(J,1)=0  THEN  J=J-1 

290  J=J  +  1 

300  NEXT  I 
310  N=J-1 

320  FOR  1=1  TO  N  :  X ( I ) =ABS (S ( I , 1 ) )  :  NEXT  I 
330  M=1  :  N1=N  :  FOR  1=1  TO  N1  :  P(I)=I  :  NEXT  I  :  GOSUB 
10000 
340  JS=1 

350  FOR  K=2  TO  N 

360  R=ABS ( X (P ( JS) )/X(P(K) ) ) 
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370  IF  R  <  .999  OR  R  >  1.001  THEN  JE=K-1  :  GOSUB  14000 

:  JS=K 
380  NEXT  K 

390  JE=N  :  GOSUB  14000 
400  TP=0:TN=0 
410  FOR  1=1  TO  N 

420  IF  Stl„l)  >0  THEN  TP=TP+S (1,1)  ELSE 

TN=TN+ABS (S  < 1 , 1) ) 

430  NEXT  I 

440  IF  TP  >  TN  THEN  TS=TN  ELSE  TS=TP 

450  CLS: PRINT320, "WILCOXSON  SIGNED  RANKS  TEST" : PRINT: PRINT 
460  PRINT  "HYPOTHESIS  NULL 


ALTERNATIVE" 
470  PRINT  " 

(1) 

E  <  X  > 

E  <Y> 

E  ( X ) 

E  (Y)  " 

480  PRINT  " 

(2) 

E(X)  <= 

E  ( Y) 

E  (  X ) 

E  (Y)  " 

490  PRINT  " 

(3) 

E  ( X )  => 

E  (Y> 

E  (X) 

E  (Y)  " 

500  PRINT:  PRINT  "Where  X  and  Y  represent  Data  Sets  1  Sc  2 
Respectively" 

510  JP=640:KS=3: GOSUB  9000 

520  IF  10=1  V=TS  ELSE  IF  10=2  V=TP  ELSE  IF  10=3  V=TN 

530  CLS:PRINT5>20,  "WILCOXSON  SIGNED  RANKS  TEST" : PRINT: PRINT 

540  PRINT  "Test  Statistic  =  ”;V 

545  PRINT  "Sample  Size  =  ";N 

550  GOSUB  2000 

560  ON  ID  GOSUB  600,800,1000 

570  PRINT  5)960,  "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

530  IF  INKEY*=  ""  GOTO  5B0  ELSE  RUN  "MENU/BAS" 

600  REM  *  *  *  NULL  HYPOTHESIS  #1  *  * 

610  PRINT  "P-Value  =  PRINT  USING  "#.####"; 2  *  FOX 

'  ’.0  PRINT:  INPUT"Enter  Desired  Alpha  Level  " ;  AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  &  1 " : GOTO  490 

640  PRINT : PRINT  "For  Alpha  =";AL;",  "; 

650  IF  AL< 2  *  FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 
660  PRINT  "  the  hypothesis  that  ECO  =  E  (Y)  " 

670  PRINT  "—vs—  the  alternative  that  E(X)  =/=  E(Y)." 

680  RETURN 

300  REM  *  *  *  NULL  HYPOTHESIS  #2  *  * 

310  FRINT  "P-Value  =  ";:FRINT  USING  "#.####"; FOX 

320  PRINT: INPUT"Enter  Desired  Alpha  Level " ; AL 

830  IF  AL<'->  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  &  1 " : GOTO  620 

840  PRINT : PRINT  "For  Alpha  =";AL;",  " ; 

950  IF  AL<FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  E(X>  <=  E(Y>" 

870  PRINT  "-vs-  the  alternative  that  E(X)  >  E(Y>." 

980  RETURN 
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lOOO  REM  t  *  *  NULL  HYPOTHESIS  #3  *  * 

1010  PRINT  "P-Value  =  ";:PRINT  USING  ####' :  FOX 

lu20  PRINT: INPUT"Enter  Desired  Alpha  Level1 ;AL 

1030  IF  AL<0  OR  AL>1  THEN  PRINT  '•***  ERROR  ***  ENTER  A 

OLAUE  BETWEEN  0  S<  1 " : GOTO  620 

1040  PRINT: PRINT  "For  Alpha  =";AL;", 

1050  IF  AL<FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 
1060  PRINT  "  the  hypothesis  that  E(X)  =  >  E ( Y ) " 

1070  PRINT  "-vs-  the  alternative  that  E(X)  <  E(Y)." 

1080  RETURN 

2000  REM  *  *  *  WILCOXSON  DISTRIBUTION  *  * 

2020  PX=V: IN=N:FLAG=0 

2030  IF  V>IN*<IN*l)/2  OR  V<0  OR  IN<1  THEN  PRINT  "  *  *  ERROR 
*  < " : GOTO  lOO 

2040  IF  VMN*(IN+l)/4  THEN  V=IN*  (IN+1 )  /2-V-l :  FLAG=1 

2050  IF  0>27  OR  (IN>9  AND  V>IN>  THEN  GOSUB  2600  ELSE  GOSUB 

2070 

2060  RETURN 

2070  IF  I NT ( V ) =V  GOSUB  2110  ELSE  GOTO  2090 
2080  RETURN 

2090  V=INT<V) : GOSUB  21 10: TFOX=FOX: V=0+1 : GOSUB 
2110: FQX= (FOX+TFOX) /2:Q0X=1-F0X 
2100  RETURN 

2110  II (2) =0: II <3>=0: II <4>=0: II (51=0: II (6)=0 
2120  C=1 

2130  N=INT(.5*(~1+SQRU+8*V))  ) 

2140  FOR  1=1  TO  N 
2150  I1(1)=0 
2160  II (1)=I1 (1)+1 
2170  IF  I < 2  GOTO  2320 
2130  I 1 (2) =1 1 < 1 ) 

2190  II <2>=I1 (2)+l 
2200  IF  I <3  GOTO  2320 
2210  1 1 (3) =11 (2) 

2220  1 1 <3)=I 1 <3>+l 
2230  IF  I < 4  GOTO  2320 
2240  I 1(4) =11 (3) 

2250  II <4)=I1 <4)+l 
2260  IF  I <5  GOTO  2320 
2270  1 1 (5) =1 1 (4) 

2230  II <5)=I1 (5) +1 
2290  IF  I<6  GOTO  2320 
2300  1 1 (6) =1 1 (5) 

2310  II <6)=I1 <6)-H 

2320  S=I1 (1)+I1 <2)+Il (3) +11 (4) +1 1 (5) +11 (6) 

2340  IF  S>V  THEN  GOTO  2360 

2350  OOl  :0N  I  GOTO  2550,2540,2530,2510,2510,2500 

2360  IF  I>1  GOTO  2380 

2370  II < 1 ) =IN: GOTO  2550 

2380  FOR  IL=1  TO  1-1 

2390  S=0:J=I-IL 
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2400  IF  J-1<=0  THEN  GOTO  2440 
2410  FOR  K=1  TO  J-l 
2420  S=S+I 1 <K) 

2430  NEX'i  K 

2440  S=S+  < IL+ 1 )  1 <J)+(IL+l)*<IL+2>/2 
2450  IF  S>V  THEN  SOTO  2480 
2460  II (I+1-IL>=IN 

2470  ON  J  GOTO  2540,2530,2520,2510,2500 

2480  NEXT  IL 

2490  II <1>=IN:G0TQ  2550 

2500  IF  II (6) < IN  GOTO  2310 

2510  IF  II (5) < IN  GOTO  2280 

2520  IF  II (4) < IN  GOTO  2250 

2530  IF  II (3) < IN  GOTO  2220 

2540  IF  II (2) < IN  GOTO  2190 

2550  IF  I 1 ( 1 ) < IN  GOTO  2160 

2560  NEXT  I 

2570  F0X=C/2CIN:Q0X=1-F0X 

2575  IF  FLAG=1  THEN  FOX=QOX : Q0X=1 -FOX 

2580  RETURN 

2590  REM  *  *  *  WILCOX SON  SIGNED  RANK  APPROX  D I ST I BUT I ON  * 

2600  MU=IN* (IN+l) /4 

2610  VAR=(2N+1)  *  (IN+DN/24 

2620  NX=V+0.5 

2630  GOSUB  7010 

2640  L4= (3NC2+3N-1 ) / ( ION* ( IN+l  > * (2N+1 > ) 

2650  W= (NX— MU) /SQR(VAR) : F= (EXP (-O. 5*WC2) ) /SQR (6. 2831S) 

2660  F3=F* (WC3— 3*W) 

2670  F0X=F0X+L4*F3 

2680  F5— — F* (WC5— 10*WC3+15*W) 

2690  F7=— F * (WE  7-21*WC5+105*WC3-105*W) 

2700 

L6=4* (3NC4+6NC3-3N+1 ) / (35* (IN* ( IN+1) * (2N+1 > ) E2) 

2710  F0X=F0X+L6*F5+.5*F7*L4C2 
2720  IF  FLAG= 1  THEN  F0X=1-F0X 
2730  QOX=l-FaX 
274C  RETURN 

7000  REM  *  *  *  NORMAL  DISTRIBUTION  *  * 

7010  SD=SGR(VAR) 

7020  x;j=  (NX-MU)  /SD:  AX=ABS  ( XN) 

7030  X2=.39Q9423*EXP( (-AX*AX) /2> 

7040  NT=1/ ( 1+. 23 164 19* AX) 

7050 

QOX=XZ* ( ( . 3193815*NT) - ( . 3565638*NTC2) + ( 1 . 781478*NTC3) -(1.321 
256*NTC4)  +  ( 1 . 330274  *NT  C  5 ) ) 

7060  IF  XN<0  GOTO  7090 

7070  FOX= 1 — QOX 

7080  RETURN 

7090  FOX=QOX 

7100  Q0X=1-F0X 

7110  RETURN 
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8000  REM  *  *  *  OPTION  SELECTION  ROUTINE  *  * 

8010  PRINT  3 J P+5, "NUMBER  OF  OPTION  DESIRED 
";LB*;CHR*(176> ;RB* 

8020  ZI=INKEY*:IF  ZI=""  THEN  8020 
8030  PRINT  3JP+31, ZIj : IO=VAL(ZI> 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  S)JP,  "***  NOT  A  VALID 
OPTION  *  *  * “ ; : FOR  1  =  1  TO  1000: NEXT: PRINT  3JP,BL*:G0T0  8010 
lOOOO  REM  *  *  *  ARRAY  SORT  *  * 

10010  L=1  :  B  ( L )  =N  1  +  1 

10020  J=B(L)  ’SET  END  OF  ARRAY  SEGMENT 
10030  I=M— 1  ’SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230  ’HANDLE  1  OR  2  ELEMENTS  IN 
SPECIAL  CASE 

10050  M1=INT< (I+J>/2)  ’SET  COMPARE  ELEMENTS 

10060  REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 

10070  1=1+1 

10080  IF  I=J  THEN  10170 

10090  IF  X  <P < I )  )  <=  X(P(MD)  THEN  10070 

10100  REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

10110  J=J-1 

10120  IF  I=J  THEN  10170 

10130  IF  X(P(J>)  >=  X (P (Ml > )  THEN  10110 

10140  REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

10150  H=P  ( I )  :  Pm=P(J)  :  P(J)=H  :  GOTO  10070 

10160  REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

ELEMENT 

10170  IF  I  >  Ml  THEN  1=1-1 

10180  IF  J=M1  THEN  10210 

10190  H-P(I)  :  P (1) —P (Ml)  :  P(M1)=H 

10200  REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10210  L=L+1  :  B(L)=I 

10220  GOTO  10020 

10230  IF  J-M  <  2  THEN  10270 

10240  IF  X  ( P  ( M )  )  <=  X ( P ( M+ 1 ) )  THEN  10260 
10250  H=P(M)  :  P(M)=P(M+1)  :  P(M+1)=H 

10260  REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10270  M=B (L) +1  :  L=L-1 

10280  IF  L >0  THEN  GOTO  10020 

10290  RETURN 

11000  REM  *  t  *  READ  FILE  DIRECTORY  *  * 

11010  CLS: PRINT320, "WILCOXSON  SIGNED  RANKS  TEST" : PRINT 
11020  F*="EDITDATA/DIR: "+DN* 

11030  OPEN  " I", 1,F* 

11040  INPUT  #1„F0 
11050  FOR  1=1  TO  FO 
11060  INPUT  #1,F*<I) 

11070  NEXT  I 
11080  CLOSE  1 
1 1 090  RETURN 
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12000  REM  *  *  *  PRINT  FILE  DIRECTORY  *  * 

12010  CLS:  PRINT3)20, "WILCOXSON  SIGNED  RANKS  TEST “ : PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030  FOR  J=1  TO  4 

12040  IF  I+J  — 1 >F0  THEN  12090 

12050  IF  J=1  AND  LEN (F* ( I+J-l ) ) < 12  AND  I>1  THEN 

PRINT  "" 

12060  PRINT  TAB(  16*  (J-l>  );"  C^STR*  (I+J-l );  “) 

F*( I+J-l ) ; 

12070  NEXT  J 

12080  NEXT  I 
12090  RETURN 

14000  REM  *  *  *  AVERAGE  RANKS  *  * 

14010  IF  JS=JE  THEN  S (P ( JS> , 1 ) =SGN (S (P ( JS) , 1 ) )  *  JS  : 

RETURN 

14020  AV=0  :  FOR  I=JS  TO  JE  :  AV=AV+I  :  NEXT  I  : 

AV=AV/ ( JE—JS+1 ) 

14030  FOR  I=JS  TO  JE  :  S (P ( I ) , 1 ) =SGN (S (P ( I ) , 1 ) )  *  AV  :  NEXT 
I 

14040  RETURN 

19000  REM  *  *  *  ERROR  HANDLING  ROUTINE  *  * 

19010  IF  ERR= 1 06  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 
OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 
PRESS  ENTER  TO  CONT I NUE " ; A*  :  SL=0  :  RESUME  11080 
19100  PRINT  " ERR= " ; ERR ; " ERL  =  " ; ERL  :  ON  ERROR  GOTO  O 
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PROGRAM 


MANN/BAS 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFINT  P: DEFSTR  Z : DEFDBL  F, L 

40  DIM  F* (95) , P (3000) , X (3000) , B (50) 

50  LB*=CHR*  <123); RB»=CHR*  <125): MB*=CHR*  <176): BL  *=CHR* < 252 ) 
60  REM  *  *  *  MANN-WHITNEY  TEST  *  * 

70  CLS: PRINT  320, "MANN-WHITNEY  TEST- 

75  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0-3) ?- 

SO  DN*=INKEY*: IF  DN*=""  THEN  SO  ELSE  I=VAL(DN*) 

85  IF  I<0  OR  I >3  THEN  GOTO  75 

90  SL=1 00  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  " " 

10O  PRINT: INPUT  "Enter  the  Two  Data  Files  <e.g.  3,5>";F1,F 
110  IF  F 1 >F0  OR  F2>F0  OR  F1<1  OR  F2< 1  THEN  IOC 
120  OPEN  "I", 1,F*<F1) +“: "+DN* 

130  INPUT  #1 , N 

140  FOR  1=1  TO  N  :  INPUT  #1,X<I>  :  NEXT  I 

150  CLOSE  1 

160  OPEN  "I", 1,F*(F2)+": "+DN* 

170  INPUT  41, M 

ISO  FOR  I=N+1  TO  N+M  :  INPUT  #1,X<I)  :  NEXT  I 
190  CLOSE  1 
200  MT=M 

210  Nl=N-*-M:M=l  :  FOR  1=1  TO  N1  :  P(I)=I  :  NEXT  I  :  GOSUB 

lOOOO 

220  JS=1 

230  FOR  K=2  TO  N1 

240  IF  X <P<JS) )<>X <P<K) )  THEN  JE=K-1:  GOSUB  14000:  JS=K 
250  NEXT  K 

260  JE=N1  :  GOSUB  14000 
270  TS=0 

2S0  FOR  1=1  TO  N 
290  TS=TS+X  < I ) 

300  NEXT  I 

305  UX=TS-N* (N+l) /2  :  UY=NT-UX 

310  IF  UX  <  UY  THEN  UU=UX  ELSE  UU=UY 

320  CLS:PRINT320, "MANN-WHITNEY  TEST" : PRINT: PRINT 


330  PRINT  "HYPOTHESIS 
ALTERNATIVE" 

NULL 

340  PRINT 
ECY)  " 

II 

( 1 ) 

E  (X) 

=  E  (Y) 

ECO  =/ 

350  PR  IN'1- 
E  (Y)  " 

II 

(2) 

E  ( X  > 

<=  E  < Y) 

E  (  X )  > 

360  PRINT 
E  <  Y)  " 

II 

(3) 

E  <X) 

=>  E (Y) 

E  ( X  )  < 
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370  PRINTS  PRINT  "Where  X  and  Y  represent  Data  Sets  1  Sc  2 
Respecti vel y " 

380  JP=640: KS=3: GOSUB  8000 

385  IF  10=1  THEN  U=UU  ELSE  IF  10=2  THEN  U=UX  ELSE  IF  10=3 
THEN  U=UY 

390  CLS: PR I NT320, "MANN: -WHITNEY  TEST" : PRINT: PRINT 
400  PRINT  "Test  Statistic  =  ";U 

405  PRINT  "Sample  Sizes  are  “;N;: PRINT  "  and  MT 

410  GOSUB  2000 

420  ON  10  GOSUB  600,800,1000 

430  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

440  IF  INKEY*=  ""  GOTO  440  ELSE  RUN  "MENU/BAS" 

600  REM  *  *  *  NULL  HYPOTHESIS  #1  *  * 

610  PRINT  "P- Value  =  SPRINT  USING  "#.####"; 2  *  FOX 

620  PRINT: INPUT"Enter  Desired  Alpha  Level";AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  0  &  1":G0T0  620 

640  PRINT: PRINT  "For  Alpha  = " ; AL ; " ,  " ; 

650  IF  AL<2  *  FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 
660  PRINT  "  the  hypothesis  that  E(X)  =  E (Y) " 

670  PRINT  "—vs—  the  alternative  that  E(X)  =/=  E(Y).“ 

680  RETURN 

800  REM  *  *  *  NULL  HYPOTHESIS  #2  *  * 

810  PRINT  "P-Value  =  ";:PRINT  USING  "#.####"; FOX 

820  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

830  IF  AL<0  OR  AL>1  THEN  PF  NT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  &  1":G0T0  620 

840  PRINT: PRINT  "For  Alpha  =";AL;",  ”; 

850  IF  AL<F0X  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  E(X)  <=  E ( Y ) " 

370  PRINT  "—vs—  the  alternative  that  E(X)  >  E(Y>." 

880  RETURN 

lOOO  REM  *  *  *  NULL  HYPOTHESIS  #3  *  * 

1010  PRINT  "P— Value  =  PRINT  USING  "#.####"; FOX 

1020  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

1030  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A 

VALUE  BETWEEN  O  &  1":G0T0  620 

1040  PRINT: PRINT  "For  Alpha  =" ; AL; " ,  " ; 

1050  IF  AL<FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 
1060  PRINT  ”  the  hypothesis  that  E(X)  =>  E (Y) " 

1070  PRINT  "-vs-  the  alternative  that  E(X)  <  E ( Y ) . " 

1080  RETURN 

2000  REM  *  *  *  MANN-WHITNEY  DISTRIBUTION  *  * 

2010  M=MT 

2020  IF  M>N  THEN  M=N:N=MT 
2030  FLAG=0 

2040  IF  U>N/2— . 5  THEN  U=N— U— 1 : FLAG=1 
2050  IN=N+M:V=TS:I=N 

2060  IF  U>N  OR  U<0  OR  N< 1  OR  M< 1  THEN  PRINT  " 

*" :GOTO  100 


*  *  ERROR 


2070  IF  V>27  OR  M>9  THEN  GOSUB  2600  ELSE  GOSUB  2090 
2080  RETURN 

2090  IF  I NT  < V) =V  THEN  GOSUB  2130  ELSE  GOTO  2110 
2100  RETURN 

2110  V=INT<V) : GOSUB  2130: TFOX=FOX: V=V+1 : GOSUB 
2130: FOX* (FOX+TFOX ) /2: QOX= 1-FOX 
2120  RETURN 

2130  II (2) =0: II (3) *0: II <4)=0: II <5>=0: II (6>=0 

2140  C=0 

2150  I1(1>=0 

2160  II (1)*I1 (1)+1 

2170  IF  I<2  GOTO  2310 

2 ISO  II (2) =1 1  (1) 

2190  II <2>=I1 <2)+l 
2200  IF  I <3  GOTO  2320 
2210  II < 3 ) = 1 1 (2) 

2220  II (3) =11 (3) +1 
2230  IF  I <4  GOTO  2320 
2240  II (4) =1 1 (3) 

2250  II (4) =1 1 (4) +1 
2260  IF  I <5  GOTO  2320 
2270  II (5) =1 1 <4) 

2280  II <5>=*I1 (5)+l 
2290  IF  I <6  GOTO  2320 
2300  I 1 (6) *11 (5) 

2310  II  <6>=U  (6)+l 

2320  S=I1  <1)+I1  (2) +11  (3) +11  (4) +11  (5) +11  (6) 

2330  IF  S>V  THEN  GOTO  2350 

2340  C=C+1:0N  I  GOTO  2540,2530,2520,2510,2500,2490 

2350  IF  I>1  GOTO  2370 

2360  II (1)=IN:G0T0  2540 

2370  FOR  IL= 1  TO  1-1 

2380  S=0: J=I— IL 

2390  IF  J-I<=0  THEN  GOTO  2430 

2400  FOR  K*1  TO  J-l 

2410  S=S+I1(K) 

2420  NEXT  K 

2430  S=S+(IL+1) 1 <J)+<IL+l)*(IL+2)/2 
2440  IF  S>V  THEN  GOTO  2470 
2450  II (I+1-IL)=IN 

2460  ON  J  GOTO  2530,2520,2510,2500,2490 

2470  NEXT  IL 

2480  11(1) =IN: GOTO  2540 

2490  IF  II (6) < IN  GOTO  2310 

2500  IF  II (5) < IN  GOTO  2280 

2510  IF  II (4) < IN  GOTO  2250 

2520  IF  II (3) < IN  GOTO  2220 

2530  IF  IK2XIN  GOTO  2190 

2540  IF  II ( 1 ) < IN  GOTO  2160 

2550  S=1 

2560  FOR  K=0  TO  N-l : S=S* < IN-K) / (N-K) : NEXT  K 
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2570  F0X=C/S:Q0X=1-F0X 

2580  IF  FLAG®  1  THEN  F0X=Q0X:Q0X=1-F0X 

2590  RETURN 

2600  REM  *  *  *  MANN-WHITNEY  APPROXIMATE  DISTRIBUTION  *  * 

2610  MU=N*M/2:VAR=N*M* (N+M+l ) / 12: NX=U+. 5 
2620  GOSUB  7010 

2630  L4= (MC2+NC2+IN+M*N) / (20*M*N* < IN+1 ) ) 

2640 

L6= (2* (MC4+NC4) +4*M*N* (MC2+NC2) +6* (MC2) *  <NC2) +4* (MC3+NC3) +7 
*M*N*  (M+N)  +Mt2+NC2+2*M*N-IN>  /  (210*  (MC2)  *  (NC2)  *  <  ( IN+1 )  C2)  ) 
2650  W= (NX— MU) /SQR(VAR) : F= (EXP (-0. 5*WC2) ) /SQR (6. 28318) 

2660  F3=F*(WC3-3#W> 

2670  F5=-F* (WC5-10*WC3+15*W> 

2680  F7=-F* ( WC 7-21 *WC 5+1 05*WC 3—1 05*W) 

2690  F0X=F0X+L4  *  F3 

2700  F0X=FQX+L6  *  F5 

2710  FOX=FOX+. 5  *  F7  *  L4C2 

2720  IF  FLAG=1  THEN  FOX =1 -FOX 

2730  00X=1— FOX 

2740  RETURN 

7000  REM  *  *  *  NORMAL  DISTRIBUTION  *  * 

7010  SD=SQR(VAR) 

7020  XN=(NX-MU) /SD: AX=ABS(XN) 

7030  XZ=.3989423*EXP( (-AX*AX) /25 
7040  NT=1/ ( 1+. 23 164 19* AX) 

7050 

QOX=XZ* ( ( . 3193815*NT) - ( . 3565638*NTC2) + ( 1 . 781473*NTC3) - (1.821 
256*NTC4)+(1.330274*NTC5) ) 

7060  IF  XN<0  GOTO  7090 

7070  F0X=1-Q0X 

7080  RETURN 

7090  FOX=QOX 

7100  QOX=l— FOX 

7110  RETURN 

3000  REM  *  *  *  OPTION  SELECTION  ROUTINE  *  * 

8010  PRINT  3JP+5, "NUMBER  OF  OPTION  DESIRED 
" ; LB* ; CHR* (176); RB* 

8020  ZI= INKEY*: IF  ZI=""  THEN  8020 
8030  PRINT  S) JP+3 1 ,  Z  I ;  :  1 0=VAL  <  Z  I ) 

8040  IF  I0<1  OR  1 0 >KS  THEN  PRINT  3J P, ”***  NOT  A  VALID 
OPTION  *** "  ;  :  FOR  1  =  1  TO  lOOO:  NEXT:  PRINT  S)JP,  BL*:  GOTO  8010 
lOOOO  REM  *  *  *  ARRAY  SORT  *  * 

10010  L=1  :  B ( L ) =N 1+1 

10020  J=B(L)  *  SET  END  OF  ARRAY  SEGMENT 
10030  I=M— 1  'SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230  ’HANDLE  1  OR  2  ELEMENTS  IN 
SPECIAL  CASE 

10050  M1=INT ( ( I+J) /2)  ’SET  COMPARE  ELEMENTS 

10060  REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
10070  1=1+1 

10080  IF  I=J  THEN  10170 


10090  IF  X (P ( I ) )  <  =  X (P (Ml ) )  THEN  10070 

10100  REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

10110  J=J-1 

10120  IF  I=J  THEN  10170 

10130  IF  X<P(J))  >=  X (P (Ml ) )  THEN  10110 

10140  REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

10150  H=P < I )  :  P< I)=P(J>  :  P(J)=H  :  GOTO  10070 

10160  REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

ELEMENT 

10170  IF  I  >  Ml  THEN  1=1-1 

10180  IF  J=M1  THEN  10210 

10190  H=P(I)  :  P ( I ) =P (Ml )  :  P(M1)=H 

10200  REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10210  L=L+1  :  B(L)=I 

10220  GOTO  10020 

10230  IF  J-M  <  2  THEN  10270 

10240  IF  X (P  <M) )  <=  X  (P (M+l ) )  THEN  10260 
10250  H=P(M)  :  P(M)=P(M+1)  :  P(M+1)=H 

10260  REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10270  M=B<L)+1  :  L=L-1 

10280  IF  L>0  THEN  GOTO  10020 

10290  RETURN 

11000  REM  *  *  *  READ  FILE  DIRECTORY  *  * 

HOlO  CLS: PRINT320,  "MANN- WHITNEY  TEST”:  PR  I  NT 
11020  F*="EDITDATA/DIR: "+DN* 

11030  OPEN  "I",1,F* 

11040  INPUT  #1 , FO 
11050  FOR  1=1  TO  FO 
11060  INPUT  #1,F*<I> 

11070  NEXT  I 
11080  CLOSE  1 
11090  RETURN 

12000  REM  *  *  *  PRINT  FILE  DIRECTORY  *  * 

12010  CLS: PRINT320, "MANN— WHI TNEY  TEST" : PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030  FOR  J=1  TO  4 

12040  IF  I +J— 1 >F0  THEN  12090 

12050  IF  J=1  AND  LEN <F* < I+J-l ) ) < 12  AND  I>1  THEN 

PRINT  "" 

12060  PRINT  TAB  <  16*  <  J-l )  STR*  ( I+J-l )  ;  "  ) 

";F*< I+J-l) ; 

12070  NEXT  J 

12080  NEXT  I 
12090  RETURN 

14000  REM  *  *  t  AVERAGE  RANKS  *  * 

14010  IF  JS=JE  THEN  X(P(JS))=JS  :  RETURN 

14020  AV=0: FOR  I=JS  TO  JE:  AV=AV+I:NEXT  I:  AV=AV/ ( JE-JS+1 ) 
14030  FOR  I=JS  TO  JE: X <P < I ) ) =AV:  NEXT  I 
14040  RETURN 


74 


19000  ’ERROR  HANDLING 

19010  IF  ERR=10&  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 
OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 
PRESS  ENTER  TO  CONTINUE"; A*  :  SL=0  s  RESUME  11080 
19100  PRINT  “ ERR= " ; ERR ;  '* ERL=  “ ; ERL  :  ON  ERROR  GOTO  O 
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1 


ERQGRAM-^J _ SMiRNOy/BAS 


10  ON  ERROR  GOTO  19000 
20  DEFSTR  Z:  DEFINT  P 

30  DIM  F* (95) ,P(1500) , XI (1500) ,B<50) , X2(1500) , X (1500) 

40  LB*=CHR* (123) : RB*=CHR* ( 125) : MB*=CHR* ( 176) : BL*=CHR* (252) 
50  REM  *  *  *  SMIRNOV  TEST  *  * 

60  CLS: PRINT  320, "SMIRNOV  TEST" 

70  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0-3)?" 

SO  DN*= INKEY*:  IF  DN*="'*  THEN  80  ELSE  I=VAL(DN*) 

90  IF  I<0  OR  I >3  THEN  GOTO  70 

100  SL=1 10  :  GOSUB  11000  :  IF  SL=0  THEN  60  ELSE  GOSUB  13000 
:  PRINT  "" 

110  PRINT: INPUT  "Enter  the  Two  Data  Files  (e.g.  3,5)";F1,F2 
120  IF  FI >F0  OR  F2>FO  OR  F1<1  OR  F2< 1  THEN  110 
130  OPEN  "I", 1,F*(F1)+": "+DN* 

140  INPUT  #1 , N 

150  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 
160  CLOSE  1 

170  N 1 =N : M= 1 : FOR  1=1  TO  N: P ( I > =1 : NEXT  I : GOSUB  10000 
ISO  FOR  1=1  TO  N:X1 (I)=X(P(I) ) : NEXT  I 
190  OPEN  "I", 1,F*(F2)+": “+DN* 

200  INPUT  #1 , MM 

210  FOR  1=1  TO  MM  :  INPUT  #1,X(I)  :  NEXT  I 

220  CLOSE  1 

230  N 1 =MM : M= 1 : FOR  1=1  TO  MM: P ( I ) = I : NEXT  I : GOSUB  10000 
240  FOR  1=1  TO  MM: X2 ( I ) =X (P ( I ) ) : NEXT  I 
242  IF  MM=  >N  THEN  GOTO  250 

244  FOR  1=1  TO  N:X(I)=X1(I)  :  NEXT  I 

245  FOR  1=1  TO  MM  :  X1(I)=X2(I)  :  NEXT  I 

246  FOR  1  =  1  TO  N  :  X 2 ( I >  =  X  < I )  :  NEXT  I 
248  I =N  :  N=MM  :  MM=I 

250  1=1: J=l: XCDF=0: YCDF=0 
260  TP=0: TM=0 

270  IF  I >N  AND  J>MM  THEN  GOTO  380 

280  IF  X 1 ( I ) =X2 ( J )  THEN  1=1+1 : J=J+1 : GOTO  330 

290  IF  XI  (I)  <  X2(J)  THEN  GOTO  300  ELSE  GOTO  310 

300  IF  I  <=N  THEN  1  =  1  +  1: GOTO  330: ELSE  IF  J<  =MM  THEN 

J=J+1 : GOTO  330: ELSE  GOTO  380 

310  IF  XI  (I)  >  X2(J)  THEN  GOTO  320  ELSE  GOTO  330 

320  IF  J  <=MM  THEN  J=J+1  ELSE  IF  I  <=N  THEN  1=1+1  ELSE  GOTO 

380 

330  YCDF= ( I - 1 ) /N  :  XCDF= <3-1 ) /MM 

340  D=YCDF-XCDF 

350  IF  D<TM  THEN  TM=D 

360  IF  D >TP  THEN  TP=D 

370  GOTO  270 

380  CLS:  PR  I NT320,  "SMIRNOV  TEST"  :  PRINT:  PRINT 
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NULL 


390  PRINT  "HYPOTHESIS 
ALTERNATIVE" 


400  PRINT  " 

(1) 

F  ( X )  = 

G  (X ) 

F  ( X ) 

=  /= 

G  ( X )  " 

410  PRINT  “ 

(2) 

FCX)  <= 

G<X) 

F  ( X  > 

> 

G  ( X )  " 

420  PRINT  " 

(3) 

F  ( X )  => 

G  ( X  > 

F  ( X ) 

< 

G  <  X  >  " 

430  PRINT: PRINT  "Where  F(X>  and  G(X>  represent  Data  Sets  1 
Sc  2  Respectively" 

440  JP=640: KS=3s GOSUB  8000 

450  CLS: PR I NT320, "SMIRNOV  TEST" : PRINT: PRINT 
460  ON  10  GOTO  470,480,490 

470  IF  TP>ABS(TM)  THEN  D=TP  ELSE  D=ABS (TM) : GOTO  500 
480  D=TP: GOTO  500 
490  D=ABSCTM> 

500  PRINT  "Test  Statistic  =  ”;D 

505  PRINT  "Sample  Sizes  are  ";N;: PRINT  "  and  “ ; MM 
510  GOSUB  2000 

520  ON  10  GOSUB  600,800,1000 

530  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

540  IF  INKEY*=  ""  GOTO  540  ELSE  RUN  "MENU/BAS" 

600  REM  *  *  *  NULL  HYPOTHESIS  #1  *  * 

610  PRINT  "P-Value  =  PRINT  USING  "#.####"; 2  *  QOX 

620  PRINT:  INPUT"Enter  Desired  Alpha  Level'*;AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  Sc  1 "  :  GOTO  620 

640  PR I NT : PR I NT  "Far  Alpha  =  ";AL;  ",  " ; 

650  IF  AL<2*QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 
660  PRINT  "  the  hypothesis  that  F(X)  =  G(X)" 

670  PRINT  "-vs-  the  alternative  that  F(X)  ~/=  G(X)." 

680  RETURN 

800  REM  t  $  t  NULL  HYPOTHESIS  #2  *  * 

810  PRINT  "P-Value  =  "; SPRINT  USING  "#.####"; QOX 

820  PRINT: INPUT"Enter  Desired  Alpha  Level" ; AL 

830  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  Sc  1 "  :  GOTO  620 

840  PRINT:PRINT  "For  Alpha  =”;AL;”,  ”; 

850  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  F(X>  <=  S'X)" 

870  PRINT  "—vs—  the  alternative  that  F(X)  >  G(X)." 

880  RETURN 

lOOO  REM  tit  NULL  HYPOTHESIS  #3  *  * 

1010  PRINT  "P-Value  =  ";:PRINT  USING  "#.####";QOX 

1020  PRINT: INPUT"Enter  Desired  Alpha  Level " ; AL 

1030  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  t**  ENTER  A 

VLAUE  BETWEEN  O  Sc  1 " :  GOTO  620 

1040  PRINT : PRINT  "For  Alpha  =";AL;",  ”; 

1050  IF  AL<QOX  THEN  PRINT  "accept”;  ELSE  PRINT  "reject"; 
1060  PRINT  "  the  hypothesis  that  F<X>  =>  GiX)" 

1070  PRINT  ”— vs—  the  alternative  that  F(X)  <  G(X)." 
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1080  RETURN 

2000  REM  *  *  *  SMIRNOV  DISTRIBUTION  *  * 

2010  M=MM 

2020  IF  N<=. 1  THEN  S= (D-l / (2) ) * ( . 12+SQR <N) +  . 1 1 /SQR (N) ) 

ELSE  IF  (M/N)=INT(M/N>  THEN  S=D  *  SQR(N/(N+M>) 

+2/(3  *  SOR  <M) )  ELSE  S=D  *  SQR(N/(N+M)) 

+2/ (5  *  SQR (M) ) 

2030  QOX=(. 135) C <SC2) 

2040  IF  Q0X<0  THEN  QOX=0 
2050  IF  Q0X>1  THEN  Q0X=1 
2060  F0X=1— QOX 
2070  PX=D 
2080  RETURN 

8000  REM  *  *  *  OPTION  SELECTION  ROUTINE  *  * 

8010  PRINT  3>JP+5,  "NUMBER  OF  OPTION  DESIRED 
" ; LB$; CHR* (176) ; RB* 

8020  ZI=INKEY*: IF  21=""  THEN  8020 
8030  PRINT  5>JP+31,ZI;  :  IO=VAL<ZI> 

8040  IF  I0<1  OR  IO>KS  THEN  FRINT  5>JP,"***  NOT  A  VALID 
OPTION  ***"; :FOR  1=1  TO  lOOO: NEXT: PRINT  3JP,BL*:G0T0  8010 
lOOOO  REM  *  *  *  ARRAY  SORT  *  * 

10010  L=1  :  B(L)=N1+1 

10020  J=B(L)  ’SET  END  OF  ARRAY  SEGMENT 
10030  I=M— 1  'SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230  ’HANDLE  1  OR  2  ELEMENTS  IN 
SPECIAL  CASE 

10050  Ml  =  INTI  <I+J) /2)  ’SET  COMPARE  ELEMENTS 

10060  REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
10070  1=1+1 

10080  IF  I=J  THEN  10170 

10090  IF  X (P ( I ) )  <=  X  <P (Ml ) )  THEN  10070 

10100  REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

10110  J=J-1 

10120  IF  I=J  THEN  10170 

10130  IF  X(P(J))  >=  X ( P ( M 1 ) )  THEN  10110 

10140  REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

10150  H=P ( I )  :  P<I)=P(J>  :  P(J>=H  :  GOTO  10070 

10160  REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

ELEMENT 

10170  IF  I  >  Ml  THEN  1=1-1 

10130  IF  J=M1  THEN  10210 

10190  H=P ( I )  :  P ( I ) =P ( Ml )  :  P(M1)=H 

10200  REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10210  L=L+ 1  :  B (L) = I 

10220  GOTO  10020 

10230  IF  J-M  <  2  THEN  10270 

10240  IF  X (P (M) )  <=  X (P (M+l ) )  THEN  10260 
10250  H=P (M)  :  P<M)=P(M+1)  :  P(M+1)=H 

10260  REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10270  M=B(L)+1  :  L=L-1 
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10280  IF  L>0  THEN  GOTO  10020 
10290  RETURN 

HOOO  REM  *  *  *  READ  DATA  FILE  DIRECTORY  *  * 

11010  CLS:  PR  I  NTS>20,  "SMIRNOV  TEST"  sPRINT 
11020  F*="EDITDATA/DIR: “+DN* 

11030  OPEN  “I“,1,F* 

11040  INPUT  #1 ,  FO 
11050  FOR  1=1  TO  FO 
11060  INPUT  #1,F*<I) 

11070  NEXT  I 
11080  CLOSE  1 
11090  RETURN 

13000  REM  *  *  *  PRINT  DATA  FILE  DIRECTORY  *  * 

13010  CLS: PR I NT320, "SMIRNOV  TEST": PRINT 
13020  FOR  1=1  TO  FO  STEP  4 
13030  FOR  J=1  TO  4 

13040  IF  I+J-l >F0  THEN  13090 

13050  IF  J  =  l  AND  LEN <F* ( I+J-l >)< 12  AND  I>1  THEN 

PRINT  ,,M 

13060  PRINT  TAB (16* ( J  — 1 ) STR$ ( I+J-l ) ; " ) 

";F*< I+J-l) ; 

13070  NEXT  J 

13080  NEXT  I 

13090  RETURN 

19000  ' ERROR  HANDLING 

19010  IF  ERR=106  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  311 
OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 
PRESS  ENTER  TO  CONTINUE"; A*  :  SL=0  :  RESUME  11080 
19100  PRINT  " ERR= " ; ERR ; " ERL= " ; ERL  :  ON  ERROR  GDTD  0 


PROGRAM  7 


Ml _ X 


/BAS 


5  CLEAR  1000 

IO  ON  ERROR  GOTO  19000 

20  DEFSTR  Z:  DEFINT  P 

30  LB*=C.iR*( 123) : RB*=CHR* ( 125) : MB*=CHR* ( 176) : BL*=CHR* ( 252) 

40  DIM  F*(95) ,P( 1500) , XI ( 1500) ,B(50) , X (1500) 

50  LB* =CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL*=CHR* ( 252 ) 

60  REM  *  *  *  LILLIEFOR  TEST  FOR  NORMALITY  *  * 

70  CLS:  PRINT  5)20,  "N  ORMALITY  TEST" 

30  PRINT; PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0-3)? 

90  DN*= INKEY*: IF  DN4=""  THEN  90  ELSE  PRINT  DN*:  I=VAL(DN*) 
100  IF  I<0  OR  I >3  THEN  GOTO  80 

110  SL=1 10  :  GQSUB  11000  :  IF  SI  ~u  THEN  70  ELSE  GOSUB  12000 
:  PRINT  " " 

120  PRINT: INPUT  "Enter  the  Data  File  #”;F1 
130  IF  F 1 >F0  OR  Fl< 1  THEN  120 
140  OPEN  "I”, 1,F*(F1)+": "+DN* 

150  INPUT  #1,N 

160  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 

170  CLOSE  1 

ISO  CLS:  PRINT  35458,  ”***  ORDERING  DATA 

*  *" 

190  N1=N:M=1:F0R  1=1  TO  N; P ( I ) =1 : NEXT  I : GOSUB  10000 
200  SM=0 : SUM=0 : S2=0 
210  FOR  1=1  TO 

N: X 1 ( I ) =X (P ( I ) > :SM=SM+X1 ( 1) : DV=X1 ( I) -SUM: SUM=SUM+CV/ I : S2=S2+ 

DV* (XI (I) -SUM) : NEXT  I 

220  CLS: PRINT  "SELECT  OPTIGN  " 

230  F'RINT.-PRINT:  PRINT  "(1)  MEAN  and  VARIANCE  are  both 
unknown " 

240  PRINT:PRINT  "(2)  MEAN  is  known,  VARIANCE  is  unknown" 

250  PR I NT : PR I NT  "(3)  MEAN  and  VARIANCE  are  both  known" 

260  JP=896:KS=3: GOSUB  8000 
290  ON  10  GOTO  300,320,340 

300  MU=SM/N: VR=S2/ (N-l) :SD=SQR(VR) :XCDF=0:T=0 
310  GOTO  370 

320  CLS: INPUT  "Enter  the  MEAN”; MU 

330  VR=0: FOR  1=1  TO  N: VR=VR+ (XI ( I ) -MU) C2: NEXT 

I: SD=SQR(VR/ (N-l ) ) : GOTO  370 

340  CLS:  INPUT  "Enter  the  MEAN"; MU 

350  INPUT  "Enter  the  VARIANCE"; VR: SD=SGR(VR) 

360  IF  VR< 0  THEN  PRINT  "*  *  *  ERROR  t  %  *  VARIANCE  MUST 

BE  POSITIVE”: GOTO  350 

370  FOR  1=  1  TO  N 

330  X N= ( X 1 ( I) -MU) /SD: AX=ABS ( XN> 

390  XZ=. 3989423XP ( (~AX*AX)/2) 

400  NT  = 1 / ( 1  + . 23 16419tAX> 


30 


EBQiRAM_Zj _ NLILL/BAS 


5  CLEAR  1000 

10  ON  ERROR  GOTO  19000 

20  DEFSTR  Z:  DEFINT  P 

30  LB*=CHR* (123): RB*=CHR* (125): MB* =CHR* (176): BL*=CHR* ( 252 ) 

40  DIM  F* ( 95 ) ,P(1 500 ) , X 1 ( 1 500 )  ,  B  ( 50 )  ,  X  ( 1 500 ) 

50  LB*=CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL*=CHR* ( 252 ) 

60  REM  *  *  *  LILLIEFOR  TEST  FOR  NORMALITY  *  * 

70  CLS:  PRINT  5)20,  "N  ORMALITY  TEST" 

30  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0—3) ?  "; 

90  DN*= INKEY*: IF  DN*=“"  THEN  90  ELSE  PRINT  DN*:  I=VAL(DN*> 
100  IF  I<0  OR  I >3  THEN  GOTO  80 

110  SL= 1 1 O  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  "" 

120  PRINT: INPUT  "Enter  the  Data  File  #";F1 
130  IF  F 1 >F0  OR  F 1 < 1  THEN  120 
140  OPEN  "I", 1 ,F*(F1) +": "+DN* 

150  INPUT  #1,N 

160  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 
170  CLGSE  1 

130  CLS: PRINT  S45B,  "***  ORDERING  DATA 

*  *" 

190  N1=N. M=1 : FOR  1=1  TO  N: P ( I ) = I : NEXT  I : GOSUB  10000 
200  SM=0: SUM=0:S2=0 
210  FOR  1=1  TO 

N: X 1 ( I ) =X (P ( I ) ) : SM=SM+X 1(1): DV=X 1(1) -SUM: SUM=SUM+DV/ I : S2=S2+ 

DV * ( X 1 ( I > -SUM ) : NEXT  I 

220  CLS: PRINT  "SELECT  OPTION  " 

230  PR I NT : PR I NT : PR I NT  ” ( 1 )  MEAN  and  VARIANCE  are  both 
unknown " 

240  PRINT:PF;INT  "(2)  MEAN  is  known,  VARIANCE  is  unknown" 

250  PRINT : PR I NT  "(3)  MEAN  and  VARIANCE  are  both  known" 

260  JP=896:KS=3: GOSUB  8000 
290  ON  10  GOTO  300,320,340 

300  MIJ=SM/N:  VR=£0/  (N-l )  :  SD=SQR  <VR)  :  XCDF=0:  T=0 
310  GOTO  370 

320  CLS: INPUT  "Enter  the  MEAN" ; MU 

330  VR=0 :  FOR  1  =  1  TO  N:  VR=VR-i-  (X  1  (  I )  -MU)  C2:  NEXT 

I : SD=SQR ( VR/ ( N— 1 ) ) : GOTO  370 

340  CLS;  INPUT  "Enter  the  MEAN"; MU 

350  INPUT  "Enter  the  VAR I ANCE" ; VR: SD=SOR (VR) 

360  IF  VR< O  THEN  PRINT  ••  *  *  *  ERROR  *  *  *  VARIANCE  MUST 

BE  F'OS  I  T  I VE  "  :  GOT 0  350 

370  FOR  1=  1  TO  N 

330  XN= v  X 1 ( I ) -MU) /SD: AX=ABS (XN) 

390  XZ=. 3989423XP ( (-AX*AX) /2) 

400  NT=1/ ( 1 -t- .  23  1 64  19*AX  ) 


SC* 


410 

Q0X=XZ*  (  C.3193815*NT)-<.356563B*NTC2)+<1.781478*NTC3> -(1.821 
256*NTC4)+(1.33274*NTC5) ) 

420  IF  XN<0  THEN  Q0X=1-Q0X 

430  F0X=1-Q0X 

440  D=ABS ( XCDF— FOX ) 

450  IF  D>T  THEN  T=D 

460  XCDF=I/N 

470  D=ABS ( XCDF— FOX  > 

480  IF  D>T  THEN  T=D 
490  NEXT  I 

500  CLS: PRINT3)20. "N  ORMALITY  TEST" 

510  PR  I  NT :  PR  I  NT :  PR  I  NT  "Test  Statistic  =  '* ;  T 
520  PRINT  "Sample  Sice  =  ” ; M 
530  ON  IQ  GOSUB  2000,3000,4000 

580  PRINT  “P-Value  =  -.PRINT  USING  ####”;  QOX 

590  PRINT: INPUT"Enter  Desired  Alpha  Level" jAL 

600  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  8<  1"  .-GOTO  590 

610  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

620  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

630  PRINT  "  the  hypothesis  that  the  sample  is  -from" 

640  PRINT  "a  normal  population." 

800  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

S10  IF  INKEY*=  " "  GOTO  810  ELSE  RUN  "MENU/BAS" 

2000  REM  *  *  *  LILLIEFQR'S  DISTRIBUTION  -  MEAN  AND  VARIANCE 

UNKNOWN 

2020  DENOM=l 

2030  D=  T/DENOM 

2040  S=D* (— . Ol+SOR (N) +  . 85/SQR (N> ) 

2050  TDENOM=DENOM 
2060  GOSUB  2500 

2070  IF  FOX  <=  .85  THEN  DENOM  =  1-.24  *  FOX: GOTO  2200 
2080  IF  FOX  <  =  .9  THEN  DENOM  =  1.357  -  .66*  FOX  :  GOTO  2200 
2090  IF  FOX  <=  .95  THEN  DENOM  =  1.339  -  .64*  FOX  :  GOTO 
2200 

2100  IF  FOX  <=  .975  THEN  DENOM  =  1.795  -  1.12*  FOX  :  SOTO 
2200 

2110  DENOM  =  2.133  -  1.467*  FOX 

2200  IF  ABS (TDENOM  -  DENOM)  >  .001  THEN  GOTO  2030 

2210  RETURN 

2500  QOX  =  . 1351 ( S C 2 ) 

2510  IF  QOX  >  1  THEN  QOX  =  1 
2520  FOX  =  1  -  QOX 
2530  RETURN 


2540 

IF 

PR 

=  >  .70 

THEN 

XQ  =  0 . 743* PR 

+ 

0.2799 

GOTO 

2590 

2550 

IF 

PR 

=  >  .50 

THEN 

XQ  =  O. 546 *PR 

0.4187 

GOTO 

2590 

2560 

IF 

PR 

=  >  .30 

THEN 

XQ  =  0. 464*PR 

+ 

0.4588 

GOTO 

2590 

2570 

IF 

PR 

=  >  .20 

THEN 

XQ  =  0.485*PR 

4- 

0.4525 

GCTC 

2590 

2530 

IF 

PR 

=  >  .10 

THEN 

XQ  =  0. 587 *PR 

0. 4321 

GOTO 

2590 

2585 

XQ 

=  . 

846*PR 

+  0.  - 

4062 

81 


2590  XQ=XQ/DENQM 
2600  XG=XQ+.20/N 
2610  RETURN 

3000  REM  *  t  *  LILLIEFOR’S  DISTRIBUTION  -  MEAN  KNOWN  AND 
VARIANCE  UNKNOWN  :  *  * 

3010  DENOM  =  1  :  PX  =  D 

3020  DENOM  =  PX /DENOM 

3022  S=D* ( . 055+SQR (N) +  . 283/SQR (N) ) 

3026  TDENOM  =  DENOM 
3030  GOSUB  2500 

3040  IF  FOX  <=  .9  THEN  DENOM  =  1.83  -.8 *  FOX  :  GOTO  3200 
3050  IF  FOX  <=  .95  THEN  DENOM  =  1.47  -  .4*  FOX  :  GOTO  3200 
3060  IF  FOX  <=  .975  THEN  DENOM  =  1.85  -  .8*  FOX  :  GOTO  3200 
3070  DENOM  =  1.07 

3200  IF  ABS (TDENOM-DENGM)  >  .001  THEN  GOTO  3020 
3210  RETURN 

4000  REM  *  *  *  KOLMOGOROV  DISTRIBUTION  *  * 

4010  S=T*  (O.  12+SGF:  <N>  +  .  1 1/SQR  'N>  ) 

4020  QOX=. 1351 (SC2) 

4030  FOX= 1 — OCX 
4040  RETURN 

9000  REM  *  *  *  OPTION  SELECTION  ROUTINE  *  * 

8010  PRINT  3 JP +5, "NUMBER  OF  OPTION  DESIRED 
" ; LB4 ; CHRt  < 1 76 ) ; RB$ 

8020  21= INKEYt : IF  ZI  =  ""  THEN  8020 
8030  PRINT  .SJP+31,21;  :  IO=VAL(ZI) 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  3JP, ”***  NOT  A  VALID 
OPTION  **.*";  :FOR  1  =  1  TO  1000:  NEXT :  PRINT  S)JP,  BL$ :  GOTO  8010 
9050  RETURN 

10000  REM  *  *  *  ARRAY  SORT  *  I 
10010  L= 1  :  B(L)=N1+1 

10020  J=B(L)  ’ SET  END  OF  ARRAY  SEGMENT 
10030  I=M— 1  ” SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230  ’ HANDLE  1  OR  2  ELEMENTS  IN 

SPECIAL  CASE 

10050  Ml = I NT ( ( I +J ) /2)  ’ SET  COMPARE  ELEMENTS 

10060  REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
10070  1=1+1 


1 0080 
10090 
10100 
101  10 
10120 
10130 
10140 
10150 
10160 
ELEMENT 
10170 
10180 
1 0 1 90 
1  0200 


IF  I=J  THEN  10170 

IF  X • P ( I ) )  < =  X (P (Ml ) )  THEN  10070 

REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

J=J -1 

IF  I=J  THEN  10170 

IF  X(P(J))  >=  X ( P ( M 1 ) )  THEN  10110 

REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

H=P ( I )  :  P(I)=P(J>  :  P(J)=H  :  GOTO  10070 

REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

IF  I  >  Ml  THEN  1=1-1 
IF  J=M1  THEN  10210 

H=P ( I )  :  PCI) =P (Ml)  :  P(M1)=H 

REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 
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LARGE  ELEMENTS 

10210  L=L+1  :  B(L)=I 

10220  SOTO  10020 

10230  IF  J-M  <  2  THEN  10270 

10240  IF  X<P<M>)  <=  X (P (M+l ) )  THEN  10260 

10250  H=P(M)  :  P<M)=P(M+1>  :  P(M+1)=H 

10260  REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10270  M=B(L)+1  :  L=L-1 

102B0  IF  L>0  THEN  GOTO  10020 

10290  RETURN 

11000  REM  *  *  *  READ  DATA  FILE  DIRECTORY  *  * 

11010  CLS : PR I NT320 , " N  ORMALITY  TES  T  M : PR I NT 

11020  F*="EDITDATA/DIR: "+DN* 

11030  OPEN  "I", 1,F$ 

11040  INPUT  #1,F0 
11050  FOR  1=1  TO  FO 
11060  INPUT  #1,F*(I) 

11070  NEXT  I 
11080  CLOSE  1 
11090  RETURN 

12000  REM  *  *  *  LIST  DATA  FILE  DIRECTORY  *  * 

12010  CLS: PRINT320, "N  ORMALITY'  TES  T”: PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030  FOR  J=1  TO  4 

12040  IF  I+J-l >F0  THEN  12090 

12050  IF  J=1  AND  LEN <F* < I+J-l ) ) < 12  AND  I>1  THEN 

PRINT 

12060  PRINT  TAB ( 16* (J-l ; STR$ ( I+J-l ) ; “ ) 

” ; FS ( I+J-l ) ; 

12070  NEXT  J 

12080  NEXT  I 

12090  RETURN 

19000  ' ERROR  HANDLING 

19010  IF  £RR=  1 06  AND  (SL---20O5  OR  SL=2410  OR  SL=2710  OR  3110 
OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 
PRESS  ENTER  TO  CONTINUE"; A*  :  SL=0  :  RESUME  lioeo 
19100  PRINT  "ERR=";ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O 


S3 


PROGRAM  O 


/BAAS 


5  CLEAR  1000 

10  ON  ERROR  GOTO  19000 

20  DEFSTR  Z:  DEFINT  P 

30  LB*=CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL*=CHR* ( 252 ) 
40  DIM  F* (95) ,P(1500> , XI (1500) , B (50) , X (1500) 

50  LB*=CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL*=CHR* ( 252 ) 
60  REM  *  *  *  EXPONENTIAL  DISTRIBUTION  TEST  *  * 

70  CLS: PRINT  320, "E  XPONENT  IAL  TEST" 

80  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0-3)? 

90  DN*= INKEY*: IF  DN*=""  THEN  90  ELSE  PRINT  DN*:  I=VAL(DN*) 
lOO  IF  I<0  OR  I >3  THEN  GOTO  80 

110  SL=1 10  :  GOSUB  llOOO  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  "" 

120  PRINT: INPUT  "Enter  the  Data  File  #“;F1 
130  IF  F 1 >FQ  OR  F1<1  THEN  120 
140  OPEN  "I", 1,F*(F1)+": "+DN* 

150  INPUT  #1,N 

160  FOR  1=1  TO  N  :  INPUT  #1,X(I>  :  NEXT  I 

170  CLOSE  1 

180  CLS; PRINT  3458,  "*  *  *  ORDERING  DATA 

*  *" 

190  N1=N: M=1 : FOR  1=1  TO  N: P ( I ) =1 : NEXT  I : GOSUB  10000 
200  SM=Q 

210  FOR  1=1  TO  N: X 1 ( I ) =X (P ( I ) ) : SM=SM+X 1(1): NEXT  I : LA=N/SM 
220  CLS: PR I NT  "SELECT  OPTION  " 

230  PRINT: PRINT: PRINT  "(1)  LAMBDA  is  unknown” 

240  PRINT : PRINT  "(2)  LAMBDA  is  known" 

250  JP=896:KS=2: GOSUB  8000 

260  ON  10  GOTO  290,270 

270  CLS: INPUT  "Enter  LAMBDA" ; LA 

290  IF  LA  <=  0  THEN  PRINT  "*  *  *  ERROR  *  *  *  LAMBDA  MUST 

BE  POSITIVE" -.GOTO  270 

290  T=C: XCDF=C 

700  FCR  1=1  TO  N 

310  FOX  =  1-EXP (-LA1 ( I ) ) 

320  D=ABS ( XCDF-FQX ) 

330  IF  D  >  T  THEN  T=D 

340  XCDF  =  I/N 

350  D=ABS ( XCDF— FOX ) 

360  IF  D  >  T  THEN  T=D 
370  NEXT  I 

380  CLS: PRINT320, "E  XPONENT  IAL  TEST" 

390  PR  I NT : PR I NT : PR I NT  "Test  Statistic  =  ";T 
400  PRINT  "Sample  Size  =  " ; N 
410  ON  10  GOSUB  2000,3000 

420  PRINT  "P-Value  =  ";:PRINT  USING  "#.####“; QOX 


430  PRINT: INPUT"Enter  Desired  Alpha  Level  ";AL 

440  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  &  l":GOTO  430 

450  PRINTsPRINT  "For  Alpha  =";AL;",  "; 

460  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

470  PRINT  **  the  hypothesis  that  the  sample  is  from" 

400  PRINT  "a  normal  population.” 

490  PRINT  5)960,  "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

500  IF  INKEY*=  " “  GOTO  500  ELSE  RUN  "MENU/BAS" 

2000  REM  *  *  *  LILLIEFDR’S  DISTRIBUTION  -  MEAN  AND  VARIANCE 
UNKNOWN 

2020  DEN0M=SQR<N)+.26+.5/JQR(N) 

2040  IF  T  >  1.3442/DENOM  +  . 2/N  THEN  FOX  =  1  :  GOTO  2220 
2060  F0X=. 5: GOSUB  22S0 

2080  IF  T>XQ  THEN  F1=F0X: F2=l : GOTO  2120 
2100  F 1 =0  :  F2=F0X 
2120  FOX  =  <Fl+F2>/2 
2140  GOSUB  2280 

2160  IF  ABS (T— XQ) < . 0001  THEN  GOTO  2220 
2180  IF  T>XQ  THEN  F1=FQX  ELSE  F2=F0X 
2200  GOTO  2120 
2220  QOX= 1-FOX 


2230 

2280 

RETURN 
IF  FOX 

_  V 

.95 

THEN 

XQ 

5. 125* 

FOX 

_ 

3. 7808 

GDTO 

2460 

2300 

IF 

FOX 

_  -j. 

.90 

THEN 

XQ 

_ 

2.026* 

FOX 

0.8367 

GOTD 

2460 

2320 

IF 

FOX 

=  > 

.80 

THEN 

XQ 

= 

1. 124* 

FOX 

_ 

0.0249 

GOTO 

2460 

2340 

IF 

FOX 

—  > 

.70 

THEN 

XQ 

__ 

0.743* 

FOX 

+ 

0.2799 

GDTD 

2460 

2360 

IF 

FOX 

.  50 

THEN 

XQ 

= 

0.546* 

FOX 

+ 

0.4187 

GOTO 

2460 

2380 

IF 

FOX 

.  30 

THEN 

XQ 

0.464* 

FOX 

4* 

0 . 4588 

GOTO 

2460 

2400 

IF 

FOX 

.20 

THEN 

XQ 

_ 

0. 485* 

FOX 

+ 

0.4525 

GOTO 

2460 

2420 

IF 

FOX 

—  X 

.  10 

THEN 

XQ 

=r 

0. 587* 

FOX 

4- 

0.4321 

« 

GOTO 

2460 

2440 

XQ 

346* 

FOX 

+  0. 

4062 

2460  XQ=XQ/DENOM  +  .2/N 
2470  RETURN 

3000  REM  *  *  *  LILLIEFOR’S  DISTRIBUTION  -  LAMBA  KNOWN  :  * 

3010  S=T*<0. 12+SQR'N) +. ll/SQRCN) ) 

3020  QOX=. 135C (SC2) 

3030  IF  QOX >1  THEN  Q0X=1 
3040  F0X=1-Q0X 
3050  PX=T 
3060  RETURN 

3000  REM  *  *  *  OPTION  SELECTION  ROUTINE  *  * 
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8010  PRINT  3JP+5, "NUMBER  OF  OPTION  DESIRED 
";LB*;CHR*<176>  ;RB* 

8020  Z I  =  INKEY*:  IF  ZI  =  ,,“  THEN  8020 
8030  PRINT  3JP+31, ZI; : IO=VAL<ZI) 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  3JP, "***  NOT  A  VALID 
OPTION  ***";: FOR  1=1  TO  1000: NEXT: PRINT  8JP,BL*:G0T0  8010 
8050  RETURN 

lOOOO  REM  *  *  *  ARRAY  SORT  *  * 

10010  L=1  :  B(L)=N1+1 

10020  J=B(L>  ’SET  END  OF  ARRAY  SEGMENT 
10030  I=M— 1  ’SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230  ’HANDLE  1  OR  2  ELEMENTS  IN 
SPECIAL  CASE 

10050  Ml = INT ( ( I+J ) /2)  ’SET  COMPARE  ELEMENTS 

10060  REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
10070  1=1+1 

10080  IF  I=J  THEN  10170 

10090  IF  X  <P ( I ) )  <=  X(P!M1))  THEN  10070 

10100  REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

10110  J=J-1 

10120  IF  I=J  THEN  10170 

10130  IF  X<P<J))  >=  X (P (Ml ) )  THEN  1C110 

10140  REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

10150  H=P ( I )  :  P(I)=P(J>  :  P(J)=H  :  GOTO  10070 

10160  REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

ELEMENT 

10170  IF  I  >  Ml  THEN  1=1-1 

10180  IF  J=M1  THEN  10210 

10190  H=P ( I )  :  P ( I ) =P (Ml  >  :  P(M1)=H 

10200  REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10210  L=L+1  :  B«L)=I 

10220  GOTO  10020 

10230  IF  J-M  <  2  THEN  10270 

10240  IF  X (P (M) )  <=  X (P (M+l ) )  THEN  10260 
10250  H=P(M)  :  P(M)=P<M+1>  :  P(M+1)=H 

10260  REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT  OF 

LARGE  ELEMENTS 

10270  M=B(L)+1  :  L=L-1 

10280  IF  L>0  THEN  GOTO  10020 

10290  RETURN 

11000  REM  *  *  *  READ  DATA  FILE  DIRECTORY  *  * 

11010  CLS: PRINT320, "E  XPONENTIAL  TES  T": PRINT 
11020  F*="EDITDATA/DIR: "+DN* 

11030  OPEN  “  I  ,  1 ,  F* 

11040  INPUT  #l,FO 
11050  FOR  1=1  TO  FO 
11060  INPUT  #1 , F*  ( I ) 

11070  NEXT  I 
11080  CLOSE  1 
11090  RETURN 

12000  REM  *  *  *  LIST  DATA  FILE  DIRECTORY  *  * 


12010  CLS: PRINT320, "E  XPONENTIAL  T  E  S  T": PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030  FOR  J=1  TO  4 

12040  IF  I+J— 1 >F0  THEN  12090 

12050  IF  J=1  AND  LEN (F* ( I+J-l ) ) < 12  AND  I>1  THEN 

PRINT 

12060  PRINT  TAB ( 16* ( J-l > ) ; ” ; STR* < I+J-l);”) 

’*;F*(I+J-1> ; 

12070  NEXT  J 

12080  NEXT  I 

12090  RETURN 

19000  ' ERROR  HANDLING 

19010  IF  ERR= 106  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 
OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 
PRESS  ENTER  TO  CONTINUE" ; :  SL=0  :  RESUME  11080 
19100  PRINT  *'ERR=”  ;  ERR;  ',ERL='’ ;  ERL  :  ON  ERROR  GOTO  O 


87 


_ 2iGyisQy/gAs 


5  CLEAR  1000 

10  ON  ERROR  GOTO  19000 

20  DEFSTR  Z:  DEFINT  P 

30  LB*=CHR* (123): RB*=CHR* (125): MB*=CHR* (176): BL$=CHR* ( 252 ) 
40  DIM  F*(95) 

50  LB*=CHR*(123) : RB$=CHR$ ( 125) : MBS=CHR* ( 176) : 3L$=CHR* (252) 
60  REM  *  *  *  CHI -SQUARE  TEST  *  * 

70  CLS: PRINT  320, "C  HI-SOUARE  TEST" 

30  PR I NT : PR I NT : PR I NT  "Disk  Drive  Containing  Data  Files 
(0-3)? 

90  DN*=INKEY$: IF  DN*=M"  THEN  90  ELSE  PRINT  DNS:  I=OAL(DNS) 
100  IF  ICO  OR  I >3  THEN  GOTO  SO 

110  SL=1 10  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  " " 

120  PRINT: INPUT  "Enter  the  Data  File  #" ;F1 

130  IF  F 1 >F0  OR  F1C1  THEN  120 

140  OPEN  " I " , 1 , FS  <F1 ) +" : "+DNS 

150  T NPUT  #1,RT: INPUT  #1,CT: INPUT  #1,BT 

160  DIM  X (RT,CT) ,R (RT) ,C(CT) 

170  FOR  Rl=l  TO  RT 
ISO  FOR  Cl  =  l  TO  CT 
190  FOR  Bl  =  l  TO  BT 
200  INPUT  #1,X(R1,C1) 

210  NEXT  B1 
220  NEXT  Cl 
230  NEXT  R1 

240  N=0:DF=(RT-1> * (CT-1 ) 

250  FOR  1=1  TG  RT 
260  R ( I ) =0 
270  FOR  J=1  TO  CT 
280  R ( I ) =R ( I ) +X  ( I , J ) 

290  NEXT  J 
300  N=N+R(I) 

310  NEXT  I 

320  FOR  J=1  TO  CT 

330  C(J>=0 

340  FOR  1=1  TO  RT 

350  C(J)=C(J)+X(I, J) 

360  NEXT  I 

370  NEXT  J 

330  T =0 : CZ  =0 : CF =0 

390  FOR  1=1  TO  RT 

400  FOR  J=1  TO  CT 

410  E=R(I>*  C!J)/N 

420  T=T+( (X(I, J)-E) C2)/E 

430  IF  E  <  5  THEN  CF=CF+1 

440  IF  E  <  1  THEN  CZ=CZ+1 
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450  NEXT  J 
460  NEXT  I 

465  IF  RT =2  AND  CT=2  THEN  T=  < (ABS (X ( 1 ,  1 )  *  X (2,2) -X (1,2) 

X (2, 1) )-0.5*N) C2> *  N/(R(1)«  R<2>*  C(l>*  C(2)) 

470  CLS: PRINT  320, “C  HI-SQUARE  TEST" 

480  PRINT: PRINT: PRINT  "Test  Statistic  -  " ; T 
490  PRINT  "Degrees  o-f  Freedom  =  ";DF 

500  PRINT  "Cells  with  expected  -frequency  <  5  is  ";100 
CF/(RTt  CT);:  PRINT  "7." 

510  PRINT  "Number  o-f  cells  with  expected  -frequency  <  1  is 
" ;  CZ 

520  GOSUB  2000 

530  PRINT  "P-Value  =  PRINT  USING  "#.####"; QOX 

540  PRINT: INPUT"Enter  Desired  Alpha  Level " ; AL 

550  IF  AL<0  OR  AL>1  THEN  PRINT  ”***  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  0  &  1 "  :  GOTO  540 

560  PRINT : PRINT  "For  Alpha  =";AL;",  "; 

570  IF  ALCQOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

580  PRINT  "  the  hypothesis  that  the  column  cells" 

590  PRINT  "have  equal  probabilities." 

600  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 
main  menu"; 

610  IF  INKEY*=  "“  GOTO  610  ELSE  RUN  "MENU/BAS" 

2000  REM  *  *  *  CHI-SQUARE  DISTRIBUTION  *  * 

2020  CX=T 

2040  IF  DF>30  GOTO  2460 
2060  DP=1 

2080  FOR  I=INT (DF)  TO  2  STEP  -2 
2100  DP=DP 
2120  NEXT  I 

2140  NP =CXC ( I NT ( (DF+1 > /2) ) *EXP(-CX/2) /DP 
2160  IF  I NT (DF/2) =DF/2  GOTO  2220 
21  SO  L1=SQR (2/CX/3. 141593) 

2200  GOTO  2240 

2220  L 1 = 1 

2240  LC=1 

2260  MC= 1 

2280  D=DF 

2300  D=D-»-2 

2320  MC=MC*CX/D 

2340  IF  MC< . OOOOO 1  GOTO  2400 

2360  LC=LC+MC 

2380  GOTO  2300 

2400  F0X=L1*NP*  LC 

2420  QOX= 1 -FOX 

2440  RETURN 

2460  NX=( (CX/DF) C (1/3)  -  (l-<2/(9*  DF ) ) ) ) /SQR ( 2/ ( 9*  DF) ) 
2480  MU=0  :  VAR=1  :  12=1 

3000  REM  $  *  t  NORMAL  DISTRIBUTION 

«  * 

3010  SD=SQR (VAR ) 

3020  XN= (NX-MU) /SD: AX=ABS (XN) 
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3030  XZ=. 3989423*  EXP ( <-AX*AX> /2> 

3040  NT=1/ ( 1+. 2316419*  AX) 

3050 

Q0X=XZ* ( ( . 3193815*NT) - ( . 3565638*NTC2) -Ml. 781478*NTC3) - (1 . 821 
256*NTE4) Ml . 330274*NTL5) ) 

3060  IF  XN<0  GOTO  3090 
•070  F0X  =  1— QOX 
3080  RETURN 
3090  FOX=QOX 
3100  00X=1— FOX 
3110  RETURN 

11000  REM  *  *  *  READ  DATA  FILE  DIRECTORY  *  * 

11010  CLS:  PRINT  5)20,  "C  HI-SQUARE  TES  T":  PRINT 

11020  Ft="CTABLE/DIR: "+DNt 

11030  OPEN  ” I " , 1 , Ft 

11040  INPUT  #1,F0 

11050  FOR  1=1  TO  FO 

11060  INPUT  #l,Ft(I) 

11070  NEXT  I 
11080  CLOSE  1 
11090  RETURN 

12000  REM  *  *  *  PRINT  OUT  DIRECTORY  CF  DATA  FILES  *  * 

12010  CLS:  PRINT  5)20,  ”C  HI-SQUARE  TES  T" -.PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030  FOR  J=1  TO  4 

12040  IF  I+J-1'FO  THEN  12090 

12050  IF  J  =  1  AND  LEN  (Ft  ( I+J-J. )  )  <  12  AND  I>1  THEN 

PRINT  ”” 

12060  PRINT  TAB ( 16* ( J— 1 ));"<"; STRt ( I+J-l )  ;  " ) 

" ;  Ft  ( I+J-l ) ; 

12070  NEXT  J 

12080  NEXT  I 

12090  RETURN 

19000  ?  ERROR  HANDLING 

19010  IF  ERR=106  AND  <SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 
OR  SL=3510  CR  SL=3910'  THEN  INPUT  ”N0  FILES  DN  THIS  DRIVE. 
PRESS  ENTER  TO  CONTINUE"; At  ;  3L=0  :  RESUME  11080 
19100  PRINT  "ERF:="  ;  ERR;  "ERL= :  ERL  :  ON  ERFlOR  GOTO  0 


R-FSOGR^M  1  O 


N  F*  S  T  #=*  “T  /  J3  «=*  S 


10  CLEAR  350 

20  DEFINT  I-K  :  DEFSTR  Z 
30  DEFDBL  F, L 

40  LB*=CHR*<123)  :  RB*=CHR* ( 1 25)  :  MB*=CHR* < 1 76)  : 

BL*=CHR* (252) 

50  CLS  i  PRINT  TABE20)  ;  *****  MENU  ***** 

60  PRINT" (1)  Probability  Distributions 
70  PRINT" (2)  Inverse  Probability  Distributions 
30  PRINT  "(3)  Return  to  Program  Selection  Menu" 

90  KP=704  :  KS=3  :  GQSUB  130 
100  ON  10  GOSUB  500,1000,120 
110  GOTO  50 
120  RUN  "MENU /BAS'* 

130  REM 

140  PRINT  DKP+5, "OPTION  DESIRED  " ; LB$ ; MB* ; RB$ ; 

150  Z I  =  I NKEY t  :  IF  ZI  =  ""  THEN  150ELSE  I0=VAL(ZI) 

160  PRINT  S/KP+21,  ZI;:  FOR  1  =  1  TO  100  :  NEXT 

170  IF  I0<1  OR  I  OKS  THEN  PRINT  SKP,  BL*;  :  PRINT  3)KP ,  "  *  * 

NOT  A  VALID  OPTION  **#**;:  FOR  1=1  TO  lOOO  :  NEXT  :  PRINT 
5)KP,  BL*;  :  GOTO  140 
180  RETURN 

190  PRINT  DKP+5, LB4; “  ";L*;"  or  ";R$; "  ";RB$;" 
";LB*;MB*;RB*; 

200  ZQ-INKEY*  :  IF  ZQ=“"  THEN  200ELSE  PRINT  3KF+13,  ZQ;  :  FCR 
1=1  TO  100  :  NEXT 

210  IF  ZQOL*  AND  ZQOCHR*  <  ASC  <L*>) -32)  AND  ZQ<  >R$  AND 
Z3< >CHR* < ASC  < R* ) -32 )  THEN  1 90 
220  RETURN 


500 

510 

520 

530 

540 

550 

560 

5/0 

580 

590 

600 

Dist 


REM 

CLS  : 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PR  I  NT 

PRINT 

PR  I  NT 

PRINT 

ribut 


PR 
:  PR 
"ME 
"  ( 1 
"  (2 
"  ( 
( 


INT  TAB( 13) "PROBABILITY  DISTRIBUTIONS 
I NT:  12=0 
NU:  " 

>  Normal  Distribution 

)  Wilcoxson  Signed-Ranks  Distribution'' 

3)  Mann-Whitney  Distribution" 

4)  Kologmorov  Distribution" 

5)  Smirnov  Distribution" 

6)  Lillifor’s  Normality  Test  Distribution" 

7)  Li  1 1  le-for7  s  E;: ponent i  al  i  ty  Test 


l  on 


610  FRINT  " (S)  Chi-Square  Distribution " 

620  KP=896  :  kS=S  :  GOSUB  130 

630  ON  10  GOSUB  2000 , 3000 , 5000 , 6000 , 7 000 , 3000 , 9000 , 1 0000 
640  IF  FLAG= 1  THEN  FLAG=0: F0X= 1 -FGX : Q0X= 1 -OCX 
650  PRINT  :  PRINT  :  PRINT  "F'rob.  i  X  <=";PX;">  =  :  PRINT 

USING"#. 4###"; FOX 

660  PRINT  "Frcb.  <  X  >  *\;PX;">  =  :  PRINT 


USING"#. ####" j  DQX 

670  PRINT  3960,"  Run  Again  or  Quit?'*;:  KP=97S  :  L*="r"  : 
R*="q"  :  GOSUB  190 

680  CLS  :  12=0  :  IF  ZQ=“r"  OR  ZQ='*R**  GOTO  630ELSE  RETURN 
690  GOTO  50 

1000  CLS  :  PRINT  TAB < 15) “ INVERSE  PROBABILITY  DISTRIBUTIONS 
1010  PRINT  3192, "MENU: PRINT" (1)  Inverse  Normal 
1020  PRINT"  (2)  Inverse  Wi  1  coxson— Si gned  Ranks'* 

1030  PRINT  "(3)  Inverse  Mann— Whitney  Distribution" 

1040  PRINT  "(4)  Inverse  Kolmogorov  Distribution" 

1050  PRINT  "(5)  Inverse  Smirnov  Distribution" 

1060  PRINT  "(6)  Inverse  Lillie-for  Normality  Test 
Di stribution" 

1070  PRINT  "(7)  Inverse  Lillie-for  Exponential  Test 
Di stribution " 

1080  PRINT  "(8)  Inverse  Chi-Square  Di str i buti on " 

1090  KP=896:KS=B: GOSUB  130 
1100  12=0 
1110  ON  10  GOSUB 

13000, 14000, 15000, 16000, 17000, 18000, 19000, 19220 
1120  PRINT  :  PRINT  "For  F(>s)  =  Prob.  <  X  <=  )  =  **  ;PR 

1130  IF  XQ=9E36  THEN  PRINT  =  infinity"  ELSE  IF  XQ=9E-36 
PRINT  ">!  =  -  infinity" 

1140  IF  XQ=9E36  OR  XQ=9E-36  THEN  GOTO  1160 

1150  PRINT  "x  =  '*;  :  IF  ABS(X0)>1000  THEN  PRINT  IMT(XQ) 

ELSE  IF  ABS ( XQ)  >100  PRINT  USING  ’####.#"; XO  ELSE  IF 
ABS  ( XQ)  >10  PRINT  USING**###.  ##";  XQ  ELSE  PRINT 
USING"##.####"; XQ 

1160  PRINT  3960,"  Run  Again  or  Quit”::  LS="r"  :  F.$="q"  : 
KP=976  :  GOSUB  190 

1170  IF  ZQ=**R"  OR  ZQ=  "r "  GOTO  1 100EL3E  RETURN 
2000  REM  *  *  *  NORMAL  DISTRIBUTION  *  * 

2010  IF  12=1  GOTO  7r'7Q 

2020  CLS: PRINT  32'  NORMAL  DISTRIBUTION" 

2030  PRINT : PRINT : .  <FUT"Enter  Mean";MU 
2040  I NPUT " En t er  Var lanes"; VAR 

2050  IF  VAR  1=0  PRINT’’:*:*  ERROR  **  ":VAR"IS  NOT  A  VALID 
VARIANCE" : GOTO  2040 

2060  INPUT"Enter  Normal  Statistic  " ;  NX  :  IF  12  =0  THEN 
PX=NX 

2070  sd=sqr:var> 

2080  XN= (NX-MU) /SD: AX=ABS CXN> 

2090  X  Z= . 3989423XP ( <-AX*AX>/2> 

2100  NT  = 1 / ( 1  + . 23 16419*AX) 

2110 

QOX  =  XZ*  (  i  ,3193815*NT)  -  ( .  3565638*NTC2>  +  <  1 . 731473*NTI3)  -  ■;  1  .  £21 
256*NTC4 ) + ( 1 . 330274CNTC5) ) 

2120  IF  X N< O  GOTO  2150 
2130  FOX= 1 -QOX 
2140  RETURN 
2150  FOX=QOX 
2160  QOX= 1 -FOX 


a— 


2170  RETURN 

3000  REM  *  *  *  WILCOXSON  SIGNED-RANKS  DISTRIBUTION  *  * 
3010  CLS: PR I NTS)  10, "WILCOXSON  SIGNED-RANKS  DISTRIBUTION" 

3020  PR I NT: PR I NT 

3030  INPUT  "ENTER  Wilcox son  Stat i sti c " ; V: PX=V 

3040  IF  V  <  O  THEN  PRINT  "*  *  *  ERROR  *  *  *  STATISTIC  MUST 

BE  POSITIVE": GOTO  3030 

3050  INPUT  "ENTER  Sample  Size"; IN 

3060  IF  IN  <  1  THEN  PRINT  "*  *  *  ERROR  *  *  *  SAMPLE  SIZE 
MUST  BE  >  1 " : GOTO  3050 

3070  IF  V>IN* (IN+1) /2  THEN  PRINT  "  *  *  ERROR  *  *  ” : GOTO 

oOoO 

3080  IF  V>IN*(IN+l>/4  THEN  V=IN* ( IN+1 ) /2— V— 1 : FLAG=1 

3090  IF  V> 27  OP.  (IN>9  AND  V>IN)  THEN  GOSUB  3620  ELSE  GOSUB 

31 10 

3100  RETURN 

3110  IF  I NT ( V ) =V  GCSUB  3150  ELSE  GOTO  3130 
3120  RETURN 

3 1 30  V=INT (V) : GOSUB  3150: TFOX=FQX : V=V+1 : GOSUB 
3150: FOX  = (FOX+TFOX ) /2: QOX= 1 -FOX 
3140  RETURN 

3150  1 1 (2) =0: 11(3) =0: 11(4) =0:  1 1 (5) =0: 11(6) =0 
3160  C=1 

3170  N= I NT ( . 5* (-1+SQR  < 1+S#V> )  ) 

3180  FOR  1=1  TO  N 
3190  U(1)=0 
3200  I1(1)=I1(1)+1 
3210  IF  I < 2  GOTC  3360 
3220  1 1 <  2) = 1 1 <  1  > 

3230  1 1 (2)  =  11 (2) +1 
3240  IF  I <3  GOTO  3360 
3250  I  1 ( 3 ) = 1 1 ( 2 ) 

3260  II  (3)  =1  1  (3)  -t-1 
3270  IF  1-4  GOTO  3360 
3280  1 1 (4) =1 1  (3) 

3290  I  1 ( 4 ) = I  1 ( 4 ) + 1 
3300  IF  I  5  GOTO  3360 
3310  II <5>=I l (4) 

3320  I  1  <5>=I 1  :5>  +1 
3330  IF  16  GOTO  3360 
3340  II '6) =11 (5) 

3330  1 1 ( 6 ) = I  1  C  6  > - 1 

3360  5=  II  ( 1 )  + 1  1  <  2  >  + 1  1  <  3 )  +  I  1  (  4 )  +  I  1  ( 5 )  ■*■  I  1  6) 

3370  IF  S  V  THEN  GOTO  3390 

33 SO  C=C+ 1 : ON  I  GOTO  3580, 3570, 3360. 3540, 3540, 3330 

3390  IF  II  GOTO  3410 

3400  I  1  •  1 ) = IN: GOTO  3580 

3410  FOR  I L= 1  TO  1-1 

3420  3=0: J= I - I L 

3430  IF  J-l  =0  THEN  GOTO  3470 

3440  FOR  K= 1  TO  J-l 

3450  S=s+Il-:K) 
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3460  NEXT  K 

3470  S-S+(IL+1> 1 (J)+(IL+l)*(IL+2)/2 
3480  IF  S>V  THEN  GOTO  3510 
3490  II  (I+l-ID-IN 

3500  ON  J  GOTO  3570,3560,3550,3540,3530 

3510  NEXT  IL 

3520  11(1 ) =IN: GOTO  3580 

3530  IF  I 1 (6) < IN  GOTO  3350 

3540  IF  II (5) < IN  GOTO  3320 

3550  IF  I 1 (4) < IN  GOTO  3290 

3560  IF  II (3) < IN  GOTO  3260 

3570  IF  II (2) < IN  GOTO  3230 

3580  IF  I 1 ( 1 ) < IN  GOTO  3200 

3590  NEXT  I 

3600  F0X=C/2CIN:QDX=1— FOX 
3610  RETURN 

3620  REM  *  *  *  WILCOXSON  SIGNED  RANK  APPROX  DISTIBUTION  * 

3630  MU=IN*(IN+1> /4 

3640  VAR=(2N+l)*(IN+l)N/24 

3650  NX-V+0.5 

3660  GOSUB  2070 

3670  L4= (3NC2+3N-1 ) / ( ION* ( IN+1 ) * (2N+1 ) ) 

3680  W= (NX— MU) /SQR(VAR) : F- (EXP (-0. 5*WC2) ) /SCR (6. 28318) 
3690  F3=F* (WC3-3*W) 

3700  FOX=FOX+L4*  F3 

3710  F5— F*(WC5-10*WC3+15*W> 

3720  F7=-F* (WC7— 21*WC5+105*WC3— 105*W> 

3730 

L6—4*  (3NC4+6NE3— 3N+1  >  /  <35*  (IN*  ( IN+1 )  *  (2N+1 )  )  C2> 

3740  FOX—FOX+L6*  F5+.5*F7*  L4E2 
3750  Q0X=1— FOX 
3760  RETURN 

5000  REM  *  *  *  MANN- WHITNEY  DISTRIBUTION  *  * 

5010  CLS: PRINT310,  "MANN- WHITNEY  DISTRIBUTION" 

5020  PRINT : PRINT 

5030  INPUT  "ENTER  Mann-Whitney  Statistic“;U:PX=U 

5040  IF  U  <  O  THEN  PRINT  "*  *  *  ERROR  *  *  *  STATISTIC 

MUST  BE  POSITIVE  "sGOTO  5030 

5050  INPUT  "ENTER  Smaller  Sample  5ize";N 

5060  IF  N  <  1  THEN  PRINT  "*  *  *  ERROR  *  *  *  SAMPLE  SIZE 
MUST  BE  >  1 GOTO  5050 

5070  INPUT  "ENTER  Larger  Sample  Size";M 

5080  IF  M  <  N  THEN  PRINT  "  *  *  *  ERROR  *  *  *  YOU  ENTERED 
THE  SMALLER  SAMPLE  SIZE  LAST  ":GOTO  5050 
5090  FLAG-0 

5100  IF  U>N*M/2— . 5  THEN  U=N*M-U- 1 : FLAG- 1 
5110  IN— N+Ms V— U+N* (N+l ) /2s I-N 

5120  IF  U>N*M  THEN  PRINT  "  *  *  ERROR  *  *  THE  STATISTIC 
MUST  BE  >  THE  PRODUCT  OF  THE  SAMPLE  SIZES" sGOTO  5030 
5130  IF  V>27  OR  M>9  THEN  GOSUB  5650  ELSE  GOSUB  5150 
5140  RETURN 

5150  IF  I NT ( V ) — V  THEN  GOSUB  5190  ELSE  GOTO  5170 
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5160  RETURN 

5170  V*INT(V> :G0SUB  5190: TFOX=FOX : V=V+ 1 : GOSUB 
5190: FOX* (FOX+TFOX) /2:Q0X=1-F0X 
5180  RETURN 

5190  II  (2)  *0:  II  (3)  =0:  1 1  (4)  =0: 1 1  (5> *0:  1 1  (6)  =0 

5200  C=0 

5210  I1(1>*0 

5220  II (1>=I1 <1>+1 

5230  IF  I<2  GOTO  5370 

5240  I 1 <2) =I1<1) 

5250  II (2) *11 (2) +1 
5260  IF  I <3  GOTO  5380 
5270  II (3) *11 (2) 

5280  II (3) *11 <3)+l 
5290  IF  I <4  GOTO  5380 
5300  I 1(4) =11 (3) 

5310  II <4)=I1 <4)+l 
5320  IF  I <5  GOTO  5380 
5330  II (5) =11 (4) 

5340  II (5) =11 (5) +1 
5350  IF  I<6  GOTO  5380 
5360  II (6) =11 (5) 

5370  II <6>=I1 <6>+l 

5380  S=I1 (1)+I1 (2) +11 (3) +11 (4) +11 (5) +11 (6) 

5390  IF  S>V  THEN  GOTO  5410 

5400  C=C+1 : ON  I  GOTO  5600,5590,5580,5570,5550,5550 

5410  IF  I>1  GOTO  5430 

5420  II (1)=IN:G0T0  5600 

5430  FOR  IL=1  TO  l-l 

5440  S=0: J=I-IL 

5450  IF  J-1<=0  THEN  GOTO  5490 

5460  FOR  K=1  TO  J-l 

5470  S=S+I 1 (K) 

5480  NEXT  K 

5490  S=S+ ( IL+1 ) 1 (J)+(IL+1) * (IL+2) /2 
5500  IF  S>V  THEN  GOTO  5530 
5510  II (I+1-IL)=IN 

5520  ON  J  GOTO  5590,5580,5570,5560,5550 

5530  NEXT  IL 

5540  II (1)=IN:G0T0  5600 

5550  IF  II (6) < IN  GOTO  5370 

5560  IF  II <5) < IN  GOTO  5340 

5570  IF  I 1 (4) < IN  GOTO  5310 

5580  IF  II (3) < IN  GOTO  5280 

5590  IF  I 1 (2) < IN  GOTO  5250 

5600  IF  1 1  ( 1 X  IN  GOTO  5220 

5610  S-l 

5620  FOR  K*0  TO  N-l : S-S* < IN-K> / (N-K) : NEXT  K 
5630  FOX=C/S: Q0X=1— FOX 
5640  RETURN 

5650  REM  *  *  *  MANN-WHITNEY  APPROXIMATE  DISTRIBUTION  *  * 
5660  MU=N/2: VAR*Nt (N+M+l 5 /12:NX=U+.5 


5670  GOSUB  2070 

5680  L4=(MC2+NC2+IN+M*N)/(20*N*(IN+1) ) 

5690 

L6- (2* (MC4+NC4) +4*N* (MC2+NC2) +6* (MC2) * (NC2) +4* (MC3+NC3) +7 
*Nt (M+N) +MC2+NC2+2*N-IN> / (210* (MC2) * (NE2) * ( ( IN+1 ) E2) ) 

5700  W“ (NX— MU) /SQR(VAR) : F* (EXP (-0. 5C2) ) /SQR (6. 28318) 

5710  F3»F*(WE3-3> 

5720  F5®>— F*  (WES— 10E3+15) 

5730  F7«— F* (WE7— 21 C5+105C3— 105) 

5740  F0X=F0X+L4*  F3 
5750  F0X*F0X+L6*  F5 
5760  F0X=F0X+.5*  F7*  L4E2 
5770  RETURN 

6000  REM  *  *  *  KOLMOGROV  DISTRIBUTION  *  * 

6010  CLS:  PRINTS)  10,  "KOLMOGOROV  DISTRIBUTION" 

6020  PRINT: PRINT 

6030  INPUT  "ENTER  SKolmogorov  Statistic";D 

6040  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  *":GOTO  6030 

6050  INPUT  "ENTER  Sample  Size";N 

6060  IF  N<1  THEN  PRINT  "*  *  ERROR  *  *":GOTO  6050 
6070  S=D* (O. 12+SQR (N) +0. 1 1/SQR (N) ) 

6080  GOSUB  6110 
6090  PX=D 
6100  RETURN 
6110  QOX=(. 135) C (SE2) 

6120  IF  QOX<0  THEN  Q0X=0 
6130  IF  Q0X>1  THEN  Q0X=1 
6140  F0X=1— QOX 
6150  RETURN 

7000  REM  *  *  *  SMIRNOV  DISTRIBUTION  *  * 

7010  CLS: PRINTS)  10,  "SMIRNOV  DISTRIBUTION" 

7020  PR I NT: PR I NT 

7030  INPUT  "ENTER  Smirnov  Statist: c";D 

7040  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  *":GOTO  7030 
7050  INPUT  "ENTER  Smaller  Sample  5ize";N 
7060  IF  N< 1  THEN  PRINT  "*  *  ERROR  *  *  SAMPLE  SIZE  MUST 
BE  >  1 " : GOTO  7050 

7070  INPUT  "ENTER  Larger  Sample  Si2e"?M 

7080  IF  M<N  THEN  PRINT  "*  *  ERROR  *  *  YOU  ENTERED  THE 

SMALLER  SAMPLE  SIZE  LAST”: GOTO  7050 

7090  IF  N<“. 1  THEN  S» <D-1 / (2) )*(. 12+SQR (N) +. 1 1/SQR (M) ) 
ELSE  IF  (M/N)=*INT(M/N>  THEN 

S=D*  SQR (N/ (N+M) ) +2/ (3*  SQR (M) ) 

ELSE  S*D*  SQR (N/ (N+M) ) +2/ (5*  SQR  CM)) 

7100  GOSUB  6110 
7110  PX=D 
7120  RETURN 

8000  REM  *  *  *  L ILL IFOR'S  TEST  FOR  NORMALITY  *  * 

8010  J»IO 

8020  CLS: PRINT  "SELECT  OPTION" : PRINT: PRINT 

8030  PRINT  "(1)  MEAN  and  VARIANCE  are  unknown" : PRINT 

8040  PRINT  "(2)  MEAN  is  known,  VARIANCE  is  unknown" 
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8050  KP=896: KS=2: GOSUB  140 

8060  CLS: PRINT310,  "LILLIFOR’S  TEST  FOR  NORMALITY" 

8070  PRINT: PRINT 

8080  INPUT  “ENTER  SKolmogorov  Statistic" |D 

8090  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  *“:QOTO  8080 

8100  INPUT  "ENTER  Sample  Size“$N 

8110  IF  N<1  THEN  PRINT  "*  *  ERROR  »  *  SAMPLE  SIZE  MUST  BE 
POSITIVE": SOTO  8100 
8120  DENOM* 1.0:PX=D 
8130  D=PX/DEN0M 

8140  IF  10=1  THEN  S=D* (-. 01+SQR (N) +0. 85/SBR <N> )  FI  RF 
S=D* <0. 055+SQR (N) +. 283/SQR (N) > 

8150  TD£NQM=DENOM 

8160  60SUB  6110 

8170  IF  ID>1  THEN  GOTO  8260 

8180  IF  FOX  <=  .85  THEN  DEN0M=l-.24*  FOX: GOTO  8230 

8190  IF  FOX  <=  .90  THEN  DEN0M=1.357-  .66*  FOX: GOTO  8230 

8200  IF  FOX  <=  .95  THEN  DENOM  =  1.339  -  .64*  FOX  :  GOTO 

8230 

8210  IF  FOX  <=  .975  THEN  DENOM  =  1.795  -  1.12*  FOX  :  GOTO 
8230 

8220  DENOM  =  2.133-1.467*  FOX 

8230  IF  ABS (TDENOM— DENOM)  >  .001  THEN  GOTO  8130 
8240  IO=J 
8250  RETURN 

8260  IF  FOX  <=  .9  THEN  DEN0M=1 . 83-. 8*  FOX  : GOTO  8230 

8270  IF  FOX  <=  .95  THEN  DENOM  =1.47  -  .4*  FOX  :  GOTO  8230 

8280  IF  FOX  <=  .975  THEN  DENOM  =  1.85  -  .8*  FOX  :  GOTO  8230 

8290  DENOM  =  1.07  :  GOTO  8230 

9000  REM  *  *  *  LILLIFOR’S  TEST  FOR  EXPONENT I AL I TY  *  * 
9010  CLS: PRINTS 10,  "LILLIFOR’S  TEST  FOR  EXPONENT I AL I TY" 

9020  PR  I  NT:  PR  I  NT 

9030  INPUT  "ENTER  SKolmogorov  Statistic" ;D 

9040  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  *":GOTQ  8080 

9050  INPUT  "ENTER  Sample  Size";N 

9060  IF  N< 1  THEN  PRINT  "*  *  ERROR  *  *  SAMPLE  SIZE  MUST  BE 
POSITIVE": GOTO  8100 

9070  IF  D  >  1.3442/ (SQR(N) +.26+. 5/SQR(N) )+.2/N  THEN  PR=1  : 
GOTO  9160 

9080  PR= .  5 :  GOSUB  19090 

9090  IF  D>XQ  THEN  P 1 =PR : P2= 1 : GOTO  9110 
9100  P1=0:P2=PR 
9110  PR=(Pl+P2)/2 
9120  GOSUB  19090 

9130  IF  ABS  (D-XQX.  0001  THEN  GOTO  9160 

9140  IF  D>XQ  THEN  P1=PR  ELSE  P2=PR 

9150  GOTO  9110 

9160  FOX=PR 

9170  Q0X=1— FOX 

9180  PX  »  D 

9190  RETURN 
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10010  IF  12*1  SOTO  10080 

10020  CLSs  PRINT  820, "CHI-SQUARE  DISTRIBUTION" 

10030  PRINT  8192, " " j : INPUT"Enter  Chi-Square  Statistic  x 
-}CX 

10040  IF  CX<=0  PRINT”**  ERROR  **  CHI-SQUARE  STATISTIC 
MUST  BE  >  O" : SOTO  10030 
10050  IF  12*0  THEN  PX*CX 

10060  INPUT"Enter  Degrees  erf  Freedom" ;  DF 

10070  IF  DFOO  OR  DFOINT(DF)  THEN  PRINT"**  ERROR  ** 

ENTER  AN  INTESER  >  0":G0T0  10060 
10080  IF  DF>30  SOTO  10290 
10090  DP*1 

10100  FOR  I* I NT (DF)  TO  2  STEP  -2 
10110  DP*DP 
10120  NEXT  I 

10130  NP=CXC ( INT ( (DF+1 ) /2) ) *EXP (-CX/2) /DP 
10140  IF  INT (DF/2)=DF/2  GOTO  10170 
10150  PC=SQR(2/CX/3. 141593) 

10160  SOTO  10180 

10170  PC=1 

10180  LC=1 

10190  MC=1 

10200  D=DF 

10210  D=D+2 

10220  MC=MC*  CX/D 

10230  IF  MC< . OOOOOl  GOTO  10260 

10240  LC=LC+MC 

10250  GOTO  10210 

10260  FOX=PC*NP*  LC 

10270  QOX=l— FOX 

10280  RETURN 

10290  NX=< (CX/DF) C (1/3)  -  (1- (2/ (9*DF) ) ) ) /EQR (2/ (9*DF) ) 

10300  MU=0  s  VAR= 1  :  12=1 
10310  GOSUB  2070 
10320  RETURN 

13000  REM  *  *  *  INVERSE  NORMAL  DISTRIBUTION  *  * 

13010  IF  12=1  GOTO  13100 

13020  CLS  :  PRINT  820,  "INVERSE  NORMAL” 

13030  PRINT  :  PRINT  :  INPUT"Enter  Mean”; MO 
13040  INPUT  "Enter  Var i ance" ; VR 

13050  IF  VR< =0  PRINT"**  ERROR  **  ";VR;"IS  NOT  A  VALID 
VAR I ANCE": GOTO  13040 

13060  INPUT"Enter  Probability  (i.e.  F(x)  =  Prob.  (  X  <=  x 
)  =  " ; PR 

13070  IF  PR*0  THEN  XQ=9E~36  ELSE  IF  PR=1  XQ=9E+36 
13080  IF  PR*0  OR  PR*1  RETURN 

13090  IF  PR<0  OR  PR>1  THEN  PRINT"**  ERROR  **  ";PR;"IS  NOT 
A  VALID  PROBABILITY": GOTO  13060 

13100  IF  PR>. 5  THEN  RL=SQR (LOG ( 1 / ( ( 1-PR) C2) ) )  ELSE 
RL=SQR  < LOG  <  1/(PRC2)>) 

13110 

XQ-RL- < (2. 515517+ ( . 802853*RL>  +  ( . 01032S*RLC2) ) / ( 1+ ( 1 * 432788*R 
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L)+(. 189269*RLC2)+(.001308*RLC3) > ) 

13120  IF  PR>«.5  THEN  XQ-M0+ (SDR <VR) *XQ)  ELSE 
XQ-M0- (SQR ( VR) *XQ) 

13130  RETURN 

14000  REM  *  *  *  INVERSE  WILCOXSON  SIGNED-RANKS  *  * 

14010  CLSs PRINT  "INVERSE  WILCOXSON  SIGNED-RANKS 
DISTRIBUTION" 

14020  PRINT! PRINT 

14030  INPUT  "ENTER  Sample  Size"; IN 

14040  IF  IN  <  1  THEN  PRINT  "*  *  *  ERROR  *  *  *  SAMPLE  SIZE 
MUST  BE  POSITIVE"! GOTO  14030 
14050  M0»IN*(IN+l>/4-.3 
14060  VR» (2N+1 ) t ( IN+1 ) N/24 

14070  GOSUB  13060 

14071  V-INT(XQ) 

14072  GOSUB  3090 

14073  TFOX=FOX 

14074  V=V+1 

14075  GOSUB  3090 

14076  XQ*INT (XQ)  +  <CDBL (PR) —TFOX) / (FOX-CDBL (TFOX) ) 

1 4080  MO-XQ-INT(XQ) 

14090  IF  MO< . 26  THEN  XQ=INT ( XQ) : RETURN 

14100  IF  MO< . 75  THEN  XQ=INT ( XQ) +. 5  ELSE  XQ=INT(XQ)+1 

14110  RETURN 

15000  REM  *  *  *  INVERSE  MANN- WHITNEY  *  * 

15010  CLS: PRINT  "INVERSE  MANN- WHITNEY  DISTRIBUTION” 

15020  PRINT: PRINT 

15030  INPUT  "ENTER  Smaller  Sample  Size'*;N 

15040  IF  N< 1  THEN  PRINT  "t  *  *  ERROR  *  *  *  SAMPLE  SIZE  MUST 
BE  POSITIVE”: GOTO  15030 

15050  INPUT  "Enter  Larger  Sample  Size";M 

15060  IF  MCI  THEN  PRINT  "*  *  *  ERROR  *  *  *  SAMPLE  SIZE  MUST 

*JE  POSITIVE":  GOTO  15050 

15070  MO=N*M/2 

15080  VR=N*M*(N+M+1> /12 

15090  GOSUB  13060 

15091  U=INT(XQ) 

15092  GOSUB  5130 

15093  TFOX=FOX 

15094  U=U+1 

15096  GOSUB  5130 

15097  XQ=INT ( XQ) + (CDBL (PR) — CDBL ( TFOX ) ) / (FOX-CDBL (TFOX ) ) 
15100  MO=XQ-INT(XQ> 

15110  IF  M0< . 26  THEN  XQ=INT (XQ) : RETURN 

15120  IF  M0< . 75  THEN  XQ=INT ( XQ) +. 5  ELSE  XQ=INT(XQ)+1 

15130  RETURN 

16000  REM  *  *  *  INVERSE  KOLMOGOROV  *  * 

16010  CLS  :  PRINT  320,  "INVERSE  KOLMOGOROV" 

16020  PRINT: PRINT 

16030  INPUT "ENTER  Sample  Size";N 

16040  IF  N<1  THEN  PRINT  "*  It  ERROR  *  fiGOTO  16030 

16050  INPUT"Enter  Probability  (i.e.  F(x)  *  Prob-  (  X  <=  x 
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)  -  " ;  PR 

16060  IF  PR-1  THEN  XQ-  0.0  ELSE  IF  PR-0  XQ-9E+36 
16070  IF  PR-0  OR  PR-1  RETURN 

16080  IF  PR<0  OR  PR>1  THEN  PRINT”**  ERROR  **  "jPR;"IS  NOT 
A  VALID  PROBABILITY": GOTO  16050 
16090  S— SQR (LOG ( 1— PR) /LOG  < . 135) ) 

16100  D— S/ (. 12+SQR <N) +. 1 1/SQR (N) > 

16110  XQ-D 
16120  RETURN 

17000  REM  *  *  *  INVERSE  SMIRNOV  *  * 

17010  CLS  :  PRINT  320,  "INVERSE  SMIRNOV" 

17020  PRINT: PRINT 

17030  INPUT  "ENTER  Smaller  Sample  Size‘*;N 

17040  IF  N<  1  THEN  PRINT  "*  *  ERROR  *  ***:GOTO  1 7030 

17050  INPUT  "ENTER  Larger  Sample  Size‘*;M 

17060  IF  M<N  THEN  PRINT  "*  *  ERROR  *  *":GOTO  17030 

17070  INPUT’*Enter  Probability  (i.e.  F(x)  -  Prob.  (  X  <=  x 

)  -  “ ; PR 

17080  IF  PR-1  THEN  XQ-  0.0  ELSE  IF  PR-0  THEN  XQ-9E+36 
17090  IF  PR-O  OR  PR-1  RETURN 

17100  IF  PR<0  OR  PR>1  THEN  PRINT"**  ERROR  **  "jPR;"IS  NOT 
A  VALID  PROBABILITY": GOTO  17060 
17110  S— SQR (LOG ( 1— PR) /LOG ( . 135) ) 

17120  IF  N<— O. 1 *M  THEN  D-S/ (. 12+SQR (N> +. 1 1/SQR (N) > +1 / (2*M) 
ELSE  IF  (M/N)— INT (M/N)  THEN 

D— (S— 2/ (3*  SQR (M) ) ) /SQR  <N*M/ (N+M) ) 

ELSE  D- (S-2/ (5*  SQR (M) ) ) /SQR (N*M/ (N+M) ) 

17130  XQ-D 
17140  RETURN 

18000  REM  *  *  *  INVERSE  LILLIEFOR’S  NORMALAL I TY  *  * 

18010  J-IO: CLS; PRINT  "SELECT  OPTION" : PRINT: PRINT 
18020  PRINT  "(1)  MEAN  and  VARIANCE  are  unknown" : PRINT 
18030  PRINT  "(2)  MEAN  is  known,  VARIANCE  is  unknown" 

18040  KP— 896: KS— 2: GOSUB  140 

18050  CLS: PRINT  813,  "INVERSE  LILLIEFOR'S  NORMALITY  TEST 
DISTRIBUTION" 

18060  PR I NT: PR I NT 

18070  INPUT "ENTER  Sample  Size";N 

18080  IF  N< 1  THEN  PRINT  "*  *  ERROR  *  *":GOTO  16030 

18090  INPUT"Enter  Probability  (i.e.  F(x>  -  Prob.  (  X  <*  x 

)  -  " ; PR 

18100  IF  PR-0  THEN  XQ-9E-36  ELSE  IF  PR-1  XQ-9E+36 
18110  IF  PR-0  OR  PR-1  THEN  I 0-J: RETURN 

18120  IF  PR<0  OR  PR>1  THEN  PRINT”**  ERROR  **  ";PR?"IS  NOT 
A  VALID  PROBABILITY": GOTO  13060 
18130  S— SQR (LOG ( ( 1— PR) ) /LOG (. 135) ) 

18140  IF  10-1  THEN  D-S/ (-. 01+SQR (N) +. 85/SQR (N) )  ELSE 
D-S/ ( . 055+. 2B3/SQR (N) +SQR (N) ) 

18150  IF  I0>1  THEN  GOTO  18240 

18160  IF  PR  <»  .85  THEN  DEN0M-1-.24*  FOX: GOTO  18210 

18170  IF  PR  <-  .90  THEN  DENOM-1.357-  .66*PR:G0T0  18210 

18130  IF  PR  <*  .95  THEN  DENOM  «  1.339  -  0.64*PR  :  GOTO 
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1.12*PR  :  GOTO 


18210 

lfc.190  IF  PR  <  =  .975  THEN  DENOM  =  1.795  - 
18210 

18200  DENOM  *  2. 133-1 . 467*PR 
18210  XQ=D* DENOM 
18220  10  =  J 
18230  RETURN 

18240  IF  PR  <  =  .9  THEN  DENOM  =  1.83  -  .8*PR  s  GOTO  18210 
18250  IF  PR  <  =  .95  THEN  DENOM  =  1.47  -  .4*PR  s  GOTO  18210 
18260  IF  PR  <=  .975  THEN  DENOM  =  1.85  -  .8*PR  :  GOTO  18210 
18270  DENOM  =  1.07: GOTO  18210 

19000  REM  *  *  *  INVERSE  LILLIEFOR  EXPONENTIAL  TEST  *  * 
19010  CLS: PRINT  813,  "INVERSE  LILLIEFOR' S  EXPONENTIAL  TEST 
DISTRIBUTION" 

19020  PRINT: PRINT 

19030  INPUT "ENTER  Sample  Size";N 

19040  IF  N< 1  THEN  PRINT  "*  *  ERROR  *  *":GOTO  16030 

19050  INPUT“Enter  Probability  (i.e.  F(>:>  =  Prob.  (  X  <=  x 

)  =  " ;PR 

19060  IF  PR=0  THEN  XGN9E-36  ELSE  IF  PR=1  XQ=9E+36 
19070  IF  PR=0  OR  PR=1  RETURN 

19080  IF  FR<0  OR  PR>1  THEN  PRINT"**  ERROR  **  “;PR;"IS  NOT 
A  VALID  PROBABILITY": GOTO  19050 
19090  DENOM=SQR (N) +. 26+. 5/SQR  <N> 

19100  IF  PR  >  .95  THEN  XQ  =  5. 125*PR  -  3.7808  :  GOTO  19190 

19110  IF  PR  =  >  .90  THEN  XQ  =  2.026*PR  -  .8367  :  GOTO  19190 

19120  IF  PR  »>  .80  THEN  XQ  =  1.124*PR  -  .0249  :  GOTO  19190 

19130  IF  PR  =>  .7  THEN  XQ  =  . 743*PR  +  .2799  :  GOTO  19190 

19140  IF  PR  =  >  .5  THEN  XQ  *  .546*PR  +  .41B7  :  GOTO  19190 

19150  IF  PR  =>  .3  THEN  XQ  -  .464*PR  +  .4588  :  GOTO  19190 

19160  IF  PR  =>  .2  THEN  XQ  =  . 485*PR  +  .4525  :  GOTO  19190 

19170  IF  PR  =»>  .1  THEN  XQ  =  .  587*PR  +  .4321  :  GOTO  19190 

19180  XQ  =  . 846*PR  +  .4062 
19190  XQ=XQ/DENOM 
19200  XQ  =  XQ  +  . 2/N 
19210  RETURN 
19220  REM 

19230  IF  12=1  GOTO  20070 

20000  CLS  :  PRINT  820, "INVERSE  CHI-SQUARE" 

20010  PRINT  :  INPUT"Enter  Degrees  o-f  Freedom"  ;DF 
20020  IF  DF <  =0  OR  DFOINT(DF)  THEN  PRINT"**  ERROR  ** 

ENTER  AN  INTEGER  >  0":GOTO  20010 

20030  INPUT  "Enter  Probability  i.e.  F<x>  =  Prob.  (  X  <=  x 
)  »  " ; PR 

20040  IF  PR=0  THEN  XQ=0  :  RETURN 
20050  IF  PR=1  THEN  XQ=9E36  :  RETURN 

20060  IF  PR<0  OR  PR>1  THEN  PRINT"**  ERROR  **  ";PR;"IS  NOT 

A  VALID  PROBABILITY": GOTO  20030 

20070  MO-O  :  VR=1  :  12=1  :  PZ=PR 

20080  IF  DF>1  GOTO  20130 

20090  PR=. 5*  < 1-PZ) 

20100  GOSUB  13100 
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20110  XQ»XQC2  s  PR=PZ 
20120  RETURN 

20130  IF  DF>2  GOTO  20160 
20140  XQ— 2*  LOG  (1-PZ) 

20150  RETURN 
20160  G0SU6  13100 
20170  DQ-XQ 

20180  IF  DF> (2+INT (4tABS (DQ) ) )  GOTO  20280 
20190  X9»l. 0000886  -  .2237368/DF  -  .01513904/DFC2 
20200  X9=X9+ (DFC— . 5  *  DQ  *  (.4713941  -02607083/DF  - 

. 008986007 /DF  C  2 ) > 

20210  X9=X9  +  ( 1/DF  *  DQC2  *  (.0001348028  +  .01128186/DF  ♦ 
.02277679/ DFC 2) ) 

20220  X9=X9  +  (DFC-1.5  *  DQC3  *  (-.008553069  -  .01153761/DF 
-  . 01323293/DFC2) ) 

20230  X9=X9  +  (DFC -2  *  DQC4  *  (.00312558  +  . 005169654/DF  - 
. 006950356 / DF  C2 ) ) 

20240  X9=X9  +  (DFC-2.5  *  DQC5  *  (-.0008426812  + 

. 00253001 /DF  +  . 001060438/DFC2) > 

20250  X9=X9  +  (DFC-3  *  DQC6  *  (.00009780499  -  . 0014501 17/DF 
+  . OO 1 565326 /DF C  2 )  ) 

20260  XQ=DF*X9C3 
20270  RETURN 

20280  X9=l  -  2/(9*DF)  +  (4*DQC4  +  16*DQC2  -28) / ( 1215»DFC2> 
20290  X9=X9  +  (8*DQC6  +  7204DQC4  +  3126*DQC2  + 

2904) / ( 229635 *  DF  C  3 ) 

20300  X8=DQ/3  +  (-DQC3  +  3*DQ> / ( 162*DF> 

20310  X8— X8  -  (3*DQC5  +  40*DQC3  +  45*DQ) / (5832*DFC2) 

20320  X8=XB  +  (30l*DQC7  -  1517tDQC5  -32769*DQC3  - 
79349*DQ) / (7873200IDFC3) 

20330  X8=  SQR (2/DF)  *  X8 
20340  XQ=DF  *  <(X8+X9)C3) 

20350  RETURN 
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